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PREFACE 


This  report  is  the  second  issue  of  the  Air  Force  Academy  Aeronautics  Digest.  The 


first  issue,  Spring  1978.  can  be  ordered  from  the  Defense  Documentation  Center  (DDC) 


using  the  AD  Number:  ADA060207.  The  DDC  address  is  Cameron  Station,  Alexandria,  VA  22314 


Contributions  for  this  digest  represent  recent  scholarly  work  by  students  and  faculty 


of  the  Department  of  Aeronautics,  members  of  other  departments  of  the  Academy  and 


The  F.  J.  Seiler  Research  Laboratory,  researchers  directly  or  indirectly  involved  with 


USAFA-sponsored  projects  and  authors  in  fields  of  Interest  to  the  USAFA.  In  addition 


to  complete  papers,  the  Digest  also  Includes,  when  appropriate,  abstracts  of  lengthier 


reports  and  articles  published  in  other  formats.  The  editors  will  consider  for  publi 


cation  contributions  in  the  general  field  of  Aeronautics,  including 


-Aeronautical  Engineering 
— Flight  Mechanics 
--Propulsion 
— Structures 
■Fluid  Mechanics 

-Thermodynamics  and  Heat  Transfer 
-Engineering  Education 


Papers  on  other  topics  will  be  considered  on  an  individual  basis.  Contributions  should 


be  sent  to 


Editor,  Aeronautics  Digest 
DFAN 

US  Air  Force  Academy,  CO  80840 


The  Aeronautics  Digest  is  presently  edited  by  Capt  E.  J.  Jumper,  PhD,  and  Capt  M.  M 


Tower,  PhD,  with  the  cooperation  of  the  Department  of  English.  In  particular,  the  Digest 


would  like  to  thank  MaJ  J.  P.  Eaton  who  served  at  liaison  officer  from  DFENG,  gave  many 


helpful  suggestions  and  provided  the  final  editorial  review. 
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A SMOKE  TUNNEL  INVESTIGATION  OF  DYNAMIC  SEPARATION 

★ 

Arthur  C.  Deekens  and 
William  R.  Kuebler  Jr. 


Abstract 


An  NACA  0015  wing  section  was  tested  in  a smoke  tunnel  at  the  United  States  Air 
Force  Academy  Aeronautics  Laboratory.  These  tests  revealed  that  rotating  the  wing 
section  through  the  static  stall  angle  of  attack  at  high  rates  caused  separation  to  be 
delayed  until  an  incidence  of  as  much  as  48  degrees.  Since  separation  is  associated 
with  lift,  the  aerodynamic  performance  of  a wing  could, thus,  be  momentarily  increased. 
By  plotting  separation  incidence  vs.  non-dlmenslonalized  constant  angular  rate,  a pre- 
dictable relationship  between  these  parameters  resulted.  This  report  furnishes  the 
analysis  and  data  base  for  this  result. 


I.  Introduction 


Aerodynamlclsts  have  continuously  studied  ways  to  increase  the  performance  of  an 
airfoil  at  high  angles  of  attack.  Although  not  completely  understood,  dynamic  stall 
seems  to  hold  some  promise  of  increasing  performance  for  short  periods  of  time.  This 
stall  phenomenon  is  ultimately  connected  with  a difference  in  the  separation  angle  of 
attack  between  a fixed  angle  stall  (static  separation)  and  the  angle  of  stall  when  the 
angle  of  attack  is  changing  (dynamic  separation) . 

II.  Theory 

Separation  effects  can  be  studied  by  placing  a wing  section  In  a smoke  tunnel  and 
observing  the  changing  streamline  patterns  as  this  wing  section's  angle  of  attack  changes. 
When  the  airfoil  la  level,  at  an  angle  of  attack  of  zero,  the  streamlines  flow  smoothly 
around  the  airfoil.  Aa  the  angle  of  attack  Increases,  the  streamlines  constrict  along 
the  upper  surface  Indicating  that  a pressure  differential,  irfilch  gives  rise  to  lift, 
has  occurred.  At  low  angles,  the  flow  remains  attached  to  the  wing.  At  some  critical 
angle  of  attack,  however,  the  streamlines  separate  from  the  upper  surface  of  the  wing, 
giving  rise  to  a region  of  reverse  and  turbulent  flow.  The  angle  of  attack  at  which 
separation  occurs,  org  , is  of  Interest  to  the  experiment  described  here. 
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As  previously  mentioned,  there  ere  two  types  of  separation,  static  and  dynamic, 
associated  with  static  stall  and  dynamic  stall  respectively.  The  static  separation 
angle  of  attack  is  the  lowest  angle  of  attack  at  which  the  airfoil,  when  maintained  at 
this  angle  for  an  Indefinite  period  of  time,  will  experience  separation.  This  angle  is 
essentially  fixed  for  a fixed  airfoil  geometry. 

The  dynamic  separation  angle  of  attack, at  least  as  defined  for  this  experiment,  is 
the  angle  at  which  the  airfoil  first  experiences  separation  from  some  given  portion  of 
the  upper  surface  when  the  airfoil  Is  rotated  at  a constant  angular  rate.  Unlike  the 
static  separation  angle,  this  angle  of  attack  Is  not  constant  but  Is  affected  by  factors 
such  as  angular  pitching  rate,  velocity,  and  airfoil  shape. 

One  measure  of  the  performance  of  an  airfoil  can  be  represented  by  Its  lift  curve, 
which  shows  lift  coefficient  vs.  angle  of  attack,  or  CL  vs.  a . As  shown  In  Figure  1, 
the  lift  coefficient  Increases  up  to  a certain  maximum  value  and  then  drops  as  the  stall 
occurs. 
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FIGURE  1.  Lift  Curve  for  NACA  0015 
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There  Is  an  Inherent  problem  In  studying  the  effects  of  separation  because  the  flow 
does  not  separate  at  a single  point;  Instead,  the  flow  begins  separating  from  the 
trailing  edge  of  the  wing  and  progresses  forward  along  the  upper  surface  as  angle  of 
attack  Is  Increased.  Point  A In  Figure  1 Is  the  approximate  angle  of  attack  at  which 
separation  begins  near  the  trailing  edge.  Point  B represents  the  maximum  lift  angle 
of  attack  (the  stalling  angle)  after  which  the  lift  coefficient  decreases  with  an 
increase  in  a . Point  C closely  represents  separation  across  the  entire  wing  section. 


Research  has  shown  that  the  CL  for  a wing  can  be  Increased  by  exploiting  dynamic 

max 

stall  (Ref  1).  This  has  been  done  by  oscillating  the  airfoil  through  a range  of  angles 
of  attack  that  Includes  the  static  stall  angle.  Although  the  reasons  for  the  Increase 
In  lift  are  still  unclear,  aerodynamlcists  have  developed  several  theories  that  account 
for  the  Increased  performance.  They  claim  that  no  observable  change  occurs  In  the  lift 
curve  slope  as  the  wing  section  Initially  rotates  past  the  static  stall  angle  of  attack. 
The  lift  curve  simply  extends  to  higher  angles  of  attack.  If  conditions  are  correct, 
at  a certain  angle  of  attack  a separation  bubble,  followed  by  a strong  vortex,  appears 
near  the  leading  edge.  As  this  vortex  moves  rearward  It  alters  the  upper  surface 
pressure  distribution  enough  (suction)  to  produce  high  lift  values  and  a rapid  nose-down 
pitching  moment  (Ref  1).  A graph  of  normal  force  coefficient,  C^,  vs.  a taken  from 
Reference  1 Is  shown  in  Figure  2.  Since  the  normal  force  coefficient  corresponds 
roughly  to  the  lift  coefficient,  a plot  of  the  CL  vs.  a curve  associated  with  dynamic 
stall  might  appear  as  in  Figure  3. 


FIGURE  2.  Normal  Force  Coefficient  vs.  Angle  of  Attack 
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Notice  that  C Increases  with  an  Increase  In  the  sinusoidal  variation  In  angular  rate. 
L 

max 

The  shape  of  the  CL  vs. a curve  effectively  changes  from  the  normal  static  stall  shape. 
Since  the  maximum  lift  coefficient  Is  roughly  proportional  to  the  separation  angle  of 
attack,  this  experiment  Involved  analyzing  the  changes  In  separation  angle  of  attack. 

It  ha  <>een  speculated  that  the  dynamic  stall  can  best  he  studied  by  using  a fixed 
angular  rate  Instead  of  a sinusoidal  variation  In  angular  rate. 


To  study  the  results  of  various  fixed  angular  rates,  one  can  plot  a graph  of  the 
separation  angle  of  attack  vs.  the  angular  rotation  rate  of  the  airfoil  section.  The 
angular  rate  was  non-dlmenslonallzed  using  Eqn  (1)  In  an  attempt  to  collapse  data. 
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where  ^ d Is  the  non-dlmenalonal  angular  rate 
a Is  the  angular  rate  In  rad/sec 

)jc  Is  the  distance  from  the  pivot  point  to  the  leading  edge 

U Is  the  free  stream  velocity 
m 


(i) 


i 

i 


l 


A 


-i 


i 


i 


5 


USAFA-TR-79-1 


Since  the  airfoil  was  pivoted  about  the  mid-chord  point,  the  numerator  In  the  above 
equation  corresponds  to  the  velocity  of  the  leading  edge  as  It  makes  Its  circular  path. 

The  half-chord  corresponds  to  the  radius  of  the  circle  etched  by  the  leading  edge.  Then 

a ^c  - ru>  » V,  (2) 

where  r Is  the  radius 

Ui  Is  the  angular  velocity 

V.  is  the  velocity  of  the  leading  edge 

ii6t  * 

The  non-dlmenslonal  angular  rate  Is  then  equivalent  to  the  ratio  of  the  velocity  of  the 
leading  edge  to  the  velocity  of  the  air  flowing  past  the  airfoil  section. 

Accordingly,  the  purpose  of  this  experiment  was  to  examine  the  relationship  between 
the  separation  angle  of  attack  and  the  non-dimensional  angular  rate  determined  by  varying 
the  wing  section  angle  of  attack  at  a constant  rate  as  opposed  to  a sinusoidal  variation 
In  angular  rate.  This  was  accomplished  by  studying  films  of  smoke  tunnel  tests  of  this 
phenomenon  and  plotting  data  In  an  attempt  to  discover  the  relationship  between  the 
separation  angle  of  attack  and  angular  rate. 

III.  Experimental  Approach 

To  analyze  the  relationship  between  separation  and  a, an  NACA  0015  airfoil  section 
was  placed  In  a low  Reynolds  number  Kenny  Engineering  smoke  tunnel  at  the  United  States 
Air  Force  Academy  Aeronautics  Laboratory.  During  a series  of  experiments,  the  rate  of 
increase  In  angle  of  attack  was  varied  parametrically,  and  motion  pictures  were  taken  of 
thn  flow  about  the  airfoil  section. 

The  tunnel  was  run  for  three  airspeeds  over  a range  of  a's  . Figure  4 shows  a 
schematic  of  the  smoke  tunnel  arrangement  and  the  placement  of  the  airfoil.  The  airspeed 
was  varied  by  turning  a knob  on  the  front  of  the  siuc.  e tunnel,  and  the  three  airspeeds 
were  later  determined  by  film  analysis. 
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Figure  4.  Schematic  of  Smoke  Tunnel  as  Modified  for  Experiment 


Figure  5.  Method  Uaed  to  Rotate  Wing  Section  Figure  6.  Digital  Clock  Mechanism 

’ *1 ' 1 
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The  airfoil  section  was  mounted  on  a pivot  rod  at  the  mid-chord.  The  pivot  rod, in 
turn,  was  connected  to  a high-torque  electric  gearmotor  at  the  back  of  the  smoke  tunnel. 

A variable  D.C.  voltage  supply  was  connected  to  the  gearmotor,  which  could  be  actuated 
by  a toggle  switch.  The  high  torque  of  the  gearmotor  (approximately  30  foot-pounds) 
resulted  in  a nearly  constant  angular  rate  with  minimal  acceleration  time.  The, voltage 
supply,  electric  gearmotor,  and  the  connection  with  the  pivot  rod  are  shown  in  Figure  5. 
By  varying  the  voltage  to  the  gearmotor  from  3 to  12  volts  in  increments  of  3 volts,  the 
experiments  were  conducted  at  four  fixed  angular  rates.  The  glass  plate  covering  the 
smoke  tunnel  test  section  was  marked  in  increments  of  10  degrees,  from  0 to  90  degrees, 
to  indicate  the  angle  of  attack  of  the  airfoil  section  at  any  given  Instant. 

To  determine  the  correlation  of  all  the  events  in  each  test  run.  a high-speed, 
digital  display  clock  was  constructed  and  placed  adjacent  to  the  smoke  tunnel  test 
section.  A movie  camera  simultaneously  recorded  the  behavior  of  the  smoke  traces  and  the 
clock  display.  The  clock  was  made  up  of  a 5-volt  voltage  supply  connected  to  a pulse 
generator,  which,  in  turn,  was  connected  to  a cumulative  digital  display  as  shown  in 
Figure  6.  The  digital  display  was  equipped  with  a calibrating  mechanism  that,  when 
activated  by  a switch,  would  allow  the  display  to  run  for  up  to  ten  seconds  (accuracy  to 
within  .0001  second).  By  slightly  varying  the  capacitance  of  the  pulse  generator  with 
control  knobs,  the  entire  clock  display  could  be  calibrated  to  accumulate  in  increments  of 
.01  second.  The  clock  could  thus  be  used  as  a reference  for  determining  such  parameters 
as  airspeed  and  rate  of  change  in  angle  of  attack. 


The  movie  camera  used  to  film  the  tests  was  a variable  speed  Aeroflex  Model  16S, 
which  was  provided  by  the  USAF  Academy  Department  of  Instructional  Technology.  The  film 
was  ASA  400  Video  News  Film  shot  at  approximately  48  frames  per  second.  Although  the 
film  speed  gave  high  resolution  with  respect  to  the  smoke  tunnel  tests,  each  frame  of 
film  recorded  approximately  two  digits  (one  Increment)  in  the  hundredths  column  of  the 
digital  clock. 


At  the  beginning  of  each  experiment,  the  static  separation  angle  of  attack  was 
visually  determined,  and  an  airspeed  calibration  run  was  made  by  filming  a puff  of  smoke 
traversing  a given  distance  across  the  test  section  with  the  airfoil  at  zero  angle  of 
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attack.  The  first  voltage  setting  was  selected  and  the  airfoil  was  rotated  from  a ■ 0° 
to  a given  angle  of  attack.  This  process  was  repeated  for  each  of  several  higher  angles 
of  attack  until  some  angle  well  past  the  dynamic  separation  angle  was  reached.  Each 
rotation  was  filmed  with  the  digital  clock  running  to  determine  the  angular  rate.  Similar 
rotations  were  made  with  the  electric  gearmotor  set  at  the  successively  higher  voltage 
settings  previously  mentioned.  This  entire  process  was  repeated  at  each  of  the  three 
airspeeds. 

At  the  completion  of  the  experiments,  the  film  was  developed  by  the  Photographies 
Division  of  the  Department  of  Instructional  Technology.  The  film  was  then  reviewed  using 
a photo-optical  data  analyzer  (L-W  Photo,  model  224-A-l)  which  could  operate  over  a 
range  of  0 to  24  frames  per  second. 

IV.  Observations  and  Discussion 


Using  the  digital  clock  as  a reference,  the  three  airspeeds  which  varied  slightly 
over  the  experiment  were  determined  to  be  approximately  2.71,  3.38,  and  6.07  feet  per 
second.  These  yielded  Reynolds  numbers  based  on  the  chord  length  of  14,520,  18,110  and 
32,520,  Indicating  laminar  flow. 


The  first  test  runs  were  made  with  the  electric  gearmotor  set  at  3 volts.  The  airspeed 
for  this  run  was  the  lowest  possible  value,  2.71  feet  per  second.  The  3-volt  setting 
gave  an  angular  rate  of  25.38  degrees  per  second.  At  this  airspeed  the  static  separation 
angle  (separation  at  approximately  quarter  chord)  was  observed  to  be  12.7  degrees; 
however,  the  dynamic  separation  at  this  angular  rate  was  delayed  until  an  angle  of  18 
degrees.  Due  to  possible  error  In  determining  the  point  of  separation,  we  judge  that 
these  angles  were  accurate  to  within  3 degrees.  The  flow  remained  completely  attached 
until  this  angle  was  achieved.  However,  If  the  airfoil  was  rotated  at  this  rate  and 
stopped  at  any  angle  between  static  and  dynamic  separation,  the  flow  separated  within 
0.2  seconds  after  the  airfoil  motion  was  stopped.  The  film  also  revealed  that  once  the 
airfoil  rotated  beyond  the  dynamic  stall  angle,  separation  occurred  In  the  same  fashion 
as  In  the  static  case,  whether  the  airfoil  continued  to  rotate  or  not.  Excerpts  fron  the 
film  o£  this  run  are  shown  In  Figure  7. 


30 


a - 40 
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Figure  7.  Film  Excerpts  at  Low  Airspeed  with  25.38  Deg  sec 
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The  next  case  examined  was  at  the  same  airspeed  but  at  6 volts,  yielding  a pitching 
rate  of  44.18  degrees  per  second.  The  dynamic  separation  angle  of  attack  was  observed  at 
28  degrees.  When  the  airfoil  section  was  rotated  to  30  degrees,  a small  vortex  formed 
Just  behind  the  leading  edge  (approximately  ^ » .10)  as  soon  as  rotation  was  stopped. 

The  vortex  Induced  some  reverse  flow  along  the  surface  and  shed  within  .5  seconds  after 
formation.  At  angles  of  attack  above  40  degrees,  there  was  some  leakage  between  the 
airfoil  section  and  the  glass  plate  covering  the  test  section;  however,  this  appeared 
to  have  negligible  effect  on  the  airflow  about  the  airfoil.  When  the  airfoil  was  rotated 
to  higher  angles  of  attack,  another  small  vortex  formed  behind  the  leading  edge  at  an 
angle  of  41  degrees  and  shed  at  48  degrees.  A large  vortex  formed  near  the  trailing  edge 
on  the  bottom  surface  whenever  the  airfoil  section  rotated  to  an  angle  of  50  degrees  or 


more.  The  fluid  rotated  In  a clockwise  sense,  as  opposed  to  the  counterclockwise  rotation 
of  the  upper  surface  vortices;  however,  since  It  shed  Immediately  upon  formation,  there 
was  no  noticeable  reverse  flow  on  the  airfoil  Itself. 


The  angular  rate  at  the  9 volt  setting  was  68.07  degrees  per  second,  and  dynamic 
separation  occurred  at  an  angle  of  36  degrees.  A small  vortex  formed  just  behind  the 
leading  edge  when  the  airfoil  section  rotated  to  40  degrees  at  this  rate.  The  vortex 
Induced  reverse  flow,  but  shed  Immediately  after  the  pitching  motion  stopped.  When  the 
airfoil  rotated  to  higher  angles  of  attack,  a similar  vortex  formed  at  50  degrees  and  was 
shed  at  58  degrees.  At  still  higher  angles  of  attack,  a large  vortex  formed  over  the 
entire  upper  surface  and  shed  at  an  angle  of  75  degrees. 


At  a rotation  rate  of  65.88  degrees  per  second  (12  volts),  the  dynamic  separation 
angle  of  attack  was  observed  at  48  degrees.  Whenever  the  wing  section  was  rotated  beyond 
48  degrees,  either  a small  vortex,  or  a separation  bubble  formed  near  the  leading  edge. 

At  an  angle  of  90  degrees,  a large  vortex  formed  similar  to  the  one  at  the  9 volt 
setting;  this  vortex  Induced  considerable  reverse  flow  along  the  entire  upper  surface. 
Excerpts  from  the  film  of  this  experimental  condition  are  shown  In  Figure  8.  (Note  the 
significant  Increase  In  the  dynamic  separation  angle  of  attack  compared  to  Figure  7.) 

The  effect  of  angular  rate  on  separation  may  also  be  seen  In  a plot  in  Figure  9. 


When  the  airspeed  was  increased  to  3.38  feet  per  second,  the  wing  section  was  tested 


j 

v 


11 


USAFA-TR-79-I 


at  the  same  four  angular  rates  as  during  the  low  speed  tests.  At  these  angular  rates, 
dynamic  separation  occurred  at  angles  of  2l,  32,  38,  and  42  degrees  respectively.  Static 
separation  angle  of  attack  at  this  airspeed  was  measured  at  12  degrees.  Large  and  small 
vortices  were  observed  to  form  In  a manner  Identical  to  that  found  at  the  low  airspeed; 
however,  the  large  vortex  separated  Immediately  after  formation  If  the  airfoil  was  rotated 
to  80  degrees  or  more.  If  the  pitching  motlon’was  stopped  below  this  angle  of  attack, 
the  large  vortex  remained  attached  for  up  to  one  second.  A plot  of  the  results  obtained 
at  this  airspeed  Is  shown  In  Figure  10. 

Finally,  the  airspeed  was  Increased  to  6.07  feet  per  second.  The  four  dynamic 
separation  angles  were  18,  27,  31  and  39  degrees,  Indicating  a decrease  In  separation 
Incidence  with  an  Increase  In  airspeed.  Vortices  formed  In  Identical  fashion  as  at  the 
lower  airspeeds, but  they  appeared  to  shed  sooner  after  formation.  The  plot  for  these 
results  are  shown  In  Figure  11. 

There  were  a few  common  phenomena  In  all  the  test  runs.  The  higher  the  angular  rate, 
the  further  beyond  the  static  stall  angle  the  airfoil  section  rotated  before  separation 
occurred.  The  high  airspeed  reduced  the  dynamic  separation  angle  of  attack  only  5 to  10 
degrees  from  the  angles  at  the  low  airspeed.  A separation  bubble  formed  over  the  upper 
surface  whenever  the  airfoil  section  rotated  well  beyond  the  static  separation  angle  of 
attack. 


One  possible  explanation  for  dynamic  separation  occurring  at  high  angles  of  attack 
Is  the  downward  relative  velocity  component  of  the  free  stream  (relative  to  the  airfoil 
section)  near  the  leading  edge.  This  explanation,  along  with  the  suction  Induced  near 
the  trailing  edge  could  have  caused  the  flow  to  remain  attached  as  the  wing  section 
rotated.  The  large  vortex  which  appeared  over  the  upper  surface  probably  did  not  affect 
the  dynamic  stall  angle  of  attack  (as  defined  In  this  paper),  since  dynamic  separation 


began  well  before  this  vortex  formed  In  tests  tritere  the  airfoil  rotated  to  high  angles 


of  attack. 

Although  many  of  these  phenomena  are  not  yet  fully  understood,  there  appears  to  be 

a predictable  relationship  (as  shown  by  Figure  12)  between  a constant  angular  rate 

(non-dimens Iona llzed)  and  dynamic  separation, 
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V.  Conclusions  and  Recommendations 

Rotating  an  airfoil  section  through  the  static  separation  angle  of  attack  at  high 
rates  delays  separation.  The  dynamic  separation  angle  appears  to  be  directly  related  to 
angular  rate  and  Inversely  related  to  airspeed;  In  fact,  by  non-dimensional lzlng  the 
angular  rate,  the  plot  In  Figure  12  can  be  used  to  predict  the  angle  of  separation  at 
any  airspeed  In  the  experimental  range. 

Since  little  is  understood  In  the  area  of  unsteady  aerodynamics,  the  following  Is 
suggested; 

1.  Further  study  could  reveal  more  exact  relationships  between  vortex  formation, 
Induced  velocity  components,  and  dynamic  separation. 

2.  This  report  could  be  used  as  a data  base  for  future  analysis  In  this  area. 

3.  The  performance  of  a wing  at  high  angles  of  attack  could  be  greatly  improved 
with  further  analysis  and  applications  of  these  results. 
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EXPERIMENTAL  MEASUREMENT  OF  SAILPLANE  TOWROPE  LOADS 
J.  P.  Retelle,  Jr.  and  F.  H.  Zeitz 


Abstract 


Sailplane  towrope  loads  were  measured  in  flight  using  an  Instrumented  towplane  tow 
hook  and  two  trainer  sailplanes.  The  measured  loads  reflect  both  the  good  and  bad 
performance  of  student  pilots  on  calm  soaring  days,  but  do  not  Include  the  effects  of 
strong  turbulence  or  rotor  activity.  Additional  data  on  the  Influence  of  towrope  age, 
as  well  as  current  practices  of  placing  knots  In  ropes,  should  be  of  interest  to  soaring 
clubs  and  other  flight  instruction  activities.  It  was  found  that,  even  for  the  most 
violent  maneuvers  tested,  the  towrope  loads  were  less  than  15  percent  of  the  measured 
towrope  breaking  strength,  and  less  than  34  percent  of  the  gross  weight  of  the  Schweitzer 
2-33  glider.  The  results  are  meant  to  provide  a baseline  of  experimental  data  as  back- 
ground for  possible  alterations  of  the  Federal  Aviation  Administration  regulations  on 
required  towrope  strengths. 


I.  Introduction 


Current  Federal  Aviation  Administration  (FAA)  regulations  regarding  sailplane  tow- 
rope  strengths  are  quite  specific:  The  towrope  must  have  a breaking  strength  no  greater 
then  200  percent  and  no  less  than  80  percent  of  the  maximum  gross  weight  of  the  sailplane 
(F.A.R.  91.17(3))(Pef.  l).  Both  a literature  search  and  contact  with  the  FAA  reveal  that 
this  regulation  is  based  upon  a very  meager  amount  of  experimental  technical  data.  The 
current  towline  strength  regulation  appears  to  be  based  upon  early  experiments  which 
measured  towllne  loads  from  a single  gilder  during  a winch  tow.  The  experiment  described 
below  measured  towrope  loads  during  aerotows  in  order  to  provide  additional  data  to 
serve  as  background  for  possible  updating  of  the  current  regulations. 


The  requirement  on  towrope  loads  is  actually  a reflection  of  concern  for  the  struc- 
tural loads  sustained  by  the  sailplane  tow  hook  and  towplane  tow  hook.  The  sailplane  tow 
hook  and  surrounding  structure,  although  designed  for  low  weight,  is  one  of  the  strongest 
parts  of  the  aircraft  and  does  not  represent  a structural  load  problem.  The  tail  struc- 
ture of  the  towplane,  however,  Is  usually  designed  to  save  weight  and  usually  represents 
the  principle  structural  constraint  when  towing  gliders.  To  satisfy  the  requirements 
cheaply,  and  also  to  allow  commercially  available  light  aircraft  to  serve  as  towplanes,  the 
FAA  placed  the  requirements  on  the  towrope  . Weak  links  of  smaller  rope,  which  themselves 
satisfy  the  towrope  strength  requirements'*',  are  permitted  to  be  used  when  sailplanes  of 
widely  different  gross  weights  are  towed  by  the  same  towplane  and  towrope  . This  experi- 
ment measured  the  towrope  loads  during  aerotows  of  gilders  during  normal  operations. 


Captain,  USAF,  Assistant  Professor  of  Aeronautics,  DFAN 

P 
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"a  25  percent  greater  strength  link  Is  allowed  at  the  towplane  end  but  not  more  than 
200  percent  of  the  gilder  gross  weight 
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If  the  loads  collected  during  aerotows  were  lower  than  those  measured  during  winch  tows, 
perhaps  the  requirements  of  the  regulation  could  be  relaxed  somewhat,  providing  weight 
and  cos t savings  for  aerotow  operators,  while  still  providing  adequate  safety.  An 
extension  of  this  analysis  could  be  that  towplanes  currently  certified  to  tow  only  light- 
weight gliders  could  tow  heavier  gliders  if  appropriate  weak  links  were  used.  Additional 
experiments  described  below  provided  some  experimental  data  on  the  influence  of  towrope 
breaking  strength  due  to  rope  age  and  the  presence  of  knots  in  the  rope. 

II.  Description  of  Experiment 

The  tests  were  conducted  using  Schweitzer  2-32  and  2-33  gliders  at  the  United  States 
Air  Force  Academy.  These  gliders,  normally  used  for  cadet  flight  instruction,  were 
flown  during  the  test  by  a certified  flight  Instructor  and  towed  by  a Piper  Super-Cub 
aircraft  with  a 180  horsepower  engine.  Weights  and  other  data  for  all  aircraft  are 
shown  in  the  Appendix. 

Towing  loads  were  measured  in  flight  using  a strain-gage  glued  to  the  tovplane  tow 
hook.  The  strain-gage  calibration  shown  in  Figure  1,  was  obtained  on  a tensile  strength 
test  stand  in  the  Mechanical  Engineering  Laboratory  at  the  Academy.  The  strain  gage 
was  wired  as  one  leg  of  a Wheatstone  bridge, and  the  bridge  output,  measured  on  a volt- 
meter, was  calibrated  directly  in  pounds.  The  voltmeter  output  was  hand-recorded  in 
flight  by  an  observer  in  the  rear  seat  of  the  towplane. 

Flight  profiles  were  planned  to  include  the  towline  loads  produced  by  both  tfca 
typical  performance,  and  the  typical  mistakes,  of  student  pilots.  A baseline  flight 
generated  data  on  standard  loads  with  no  perturbations.  Later  flights,  described  in 
Table  1,  included  normal  training  maneuvers,  such  as  boxing  the  propwash  and  slack  line  i 

recovery,  as  well  as  maneuvers  involving  poor  tow  position  and  coordination.  Our  test 
plan  was  developed  by  dividing  the  actual  aerotow  into  four  segments:  the  ground  roll, 
the  takeoff  climbout,  the  main  tow,  and  the  release.  Each  segment  of  the  tow  was 
tested  Individually  to  determine  what  kind  of  loads  would  be  encountered  in  that  par- 
ticular portion. 


/ 

v 


The  test  for  high  loads  on  takeoff  roll  was  made  by  starting  the  takeoff  with  a 
slight  amount  of  slack  in  the  rope  rather  than  by  starting  with  the  rope  taut.  A 
second  test  for  high  loads  on  takeoff  Involved  dragging  the  skid  of  the  sailplane  during 
the  takeoff  roll.  The  takeoff  cllinbout  waa  the  portion  of  flight  lmDadlately  after 
leaving  the  ground  before  normal  towing  conditions  were  established  (somewhere  before 
the  first  turn  in  the  cadet  soaring  operation).  The  two  sources  of  high  loads  which 
were  unique  to  this  segment  of  tow  are  sideslipping  immediately  after  takeoff  as  in  a 
cross-wind,  and  a rapid  pitch-up  of  the  sailplane  to  catch  up  to  the  towplane  as  it 
begins  its  climbout.  The  high  loads  during  normal  tow  wars  caused  mainly  by  being  out 
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Figure  1.  Calibration  of  Tow  Hook  Strain  Gage  Transducer 


Table  1 

FLIGHT  PROFILES 


Flight  1 

Baseline  Flight  - Gather  data  on  standard  loads  (no  perturbations). 

Flight  2 

Slack  Rope  Takeoff 

Box  the  Wash  (Improper  Tow  Position) 

Accelerated  Release 

Flight  3 

Drag  the  skid  on  takeoff. 

Side  slip  after  takeoff. 

Slack  rope  on  tow  (level  flight). 

Accelerated  release. 

Flight  4 

Zooming  after  takeoff. 

Out  of  position  In  turns. 

Slack  line  In  turns. 

Flight  3 

Level  flight  loads 

Turbulence  loads 
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of  position  or  by  having  slack  In  the  line.  Both  of  these  situations  were  tested  In 
level  flight  and  in  turns.  The  only  high  loads  encountered  during  release  were  caused 
by  an  "accelerated"  release,  one  where  the  sailplane  was  in  a climbing  right  turn 
before  release.  The  last  flight  was  used  to  generate  towline  loads  in  level  flight  In 
an  attempt  to  provide  a measure  of  aerodynamic  drag. 


III.  Results 


The  breaking  strengths  of  various  towrope  configurations  were  determined  experi- 
mentally using  the  tensile  strength  test  stand  In  the  Academy's  Engineering  Mechanics 
Laboratory.  The  results  of  these  tests  are  shown  In  Figure  2. 


ROPE  A:  INCH  POLYPROPYLENE 
ROPE  B:  INCH  POLYPROPYLENE 


OLD  ROPE  (APPROX.  2000  FLIGHTS) 


EUROPEAN"  WEAK  LINK  (USED 


WITH  2-INCH  RING 


WITH  ONE 
OVERHAND  KNOT 


ROPE  A 
NO  KNOT 


ROPE  A 

ONE  OVERHAND 
KNOT 


Figure  2.  Breaking  Strengths  for  Various  Rope  Configurations 


The  first  set  of  data  on  rope  A,  the  3/8  inch  polypropylene  rope  normally  used  for 
Academy  towing, shows  a 27  percent  decrease  In  breaking  strength  after  the  rope  has  been 
used  for  approximately  2000  flights.  (It  should  be  noted  that  the  measured  breaking 
strength  was  40  pounds  more  than  the  manufacturer's  Indicated  breaking  strength  of  2150 
pounds.)  The  second  set  of  data  shews  the  affect  of  putting  one  overhand  knot  in  the 
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line  which  was,  until  recently,  a standard  practice  at  the  Academy,  since  the  prescribed 
rope  breaking  strength  was  more  than  200  percent  of  the  glider  gross  weight.  For  a new 
rope,  the  breaking  strength  with  one  knot  was  reduced  by  43  percent.  After  about  2000 
flights,  the  knotted  rope  had  a 42.5  percent  decrease  in  breaking  strength.  The  third 
set  of  data  are  for  rope  B,  a k,  inch  polypropylene  rope  of  1100  pound  (prescribed) 
breaking  strength  used  for  weak  links.  The  figure  shows  the  degradation  of  breaking 
strength  using  a three-inch  ring  (the  standard  configuration)  and  then  an  additional 
overhand  knot.  The  addition  of  the  knot  in  the  weak  link  produced  a 79.9  percent  de- 
crease in  breaking  strength. 

Flight  test  data,  presented  in  Figure  3 for  the  Schweitzer  2-33  glider,  compare  the 
two  rope  loads  measured  during  each  segment  of  the  baseline  flight  to  those  measured 
for  the  unusual  maneuvers  performed  during  later  flights.  The  numbers  shown  represent 
average  values,  except  where  the  data  indicated  a peak  or  impulse  force  in  which  case 
the  peak  loads  are  presented.  The  baseline  data  show  a fairly  constant  load  of  128 
pounds  average.  The  snapping  of  the  towrope  on  takeoff  roll  produced  a peak  of  a 23 
percent  increased  load, and  the  snapping  of  the  line  during  an  intentionally  unskilled 
slack-rope  recovery  produced  a peak  of  a 182  percent.  Pull-up  maneuvers  durinR  the 
cllmbout  and  release  produced  load  increases  of  47  percent  and  78  percent  respectively. 
The  load  changes  produced  by  Intentionally  poor  technique  during  the  turns  produr  1 
very  small  maximum  load  changes  of  25  percent.  Figure  4 presents  the  towrope  loads 
during  the  "boxing-the-prop-wash"  training  maneuver.  The  strongest  loads  were  noted 
when  the  sailplane  was  ascending  outside  the  propwash,  and  were  50  percent  higher  then 
in  straight  flight.  By  ascending  through  the  prop  wash,  the  sailplane  pilot  increased 
the  towrope  loads  to  80  percent  above  those  measured  in  straight  flight. 

The  towrope  loads  obtained  while  towing  a heavier  glider,  the  Schweitzer  2-32,  are 
shown  in  Figure  5.  T^e  peak  load  during  the  takeoff  roll  is  higher  than  for  the  2-33, 
as  might  be  expected.  It  is  interesting  to  note  that  the  loads  experienced  while 
maneuvering  are  a fairly  constant  14  percent  higher  than  the  loads  measured  in  straight 
flight.  The  high  drag  produced  by  the  spoilers  on  this  glider  are  also  noted. 

Data  were  collected  for  both  gliders  on  tow  to  determine  the  glider  drag  In  straight 
and  level  flight  by  towing  at  a constant  altitude  and  airspeed  for  approximately  three 
minutes.  At  a tow  speed  of  60  miles  per  hour  and  at  8500  feet  altitude,  the  Schweitzer 
2-33  produced  an  average  rope  load  of  81  pounds.  At  a tow  speed  of  70  miles  per  hour 
at  the  same  altitude,  the  Schweitzer  2-32  produced  an  average  load  of  110  pounds.  These 
loads  represent  drag  coefficients  of  • .031  and  ■ .062  for  the  2-33  and  2-32 
gliders,  respectively.  While  somewhat  Inaccurate,  these  numbers  are  not  unreasonable. 

The  rope  alone  was  found  to  produce  an  average  load  of  15  pounds  after  sailplane  release. 


< 
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TAKEOFF  ROL  LIFT-OFF 


STRAIGHT  FLIGHT 


RELEASE 


BASELINE  FLIGHT 


Figure  3.  Mae cured  Tow  Load a for  SCS  2-33 
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The  maximum  towrope  loads  measured  during  this  flight  test  were  found  to  be  less 
than  15  percent  of  the  measured  towrope  breaking  strength  and  less  than  34  percent  of 
the  2-33  glider  gross  weight.  The  maneuvers  were  performed  in  calm  air  and  were  meant 
to  represent  the  performance  and  mistakes  that  were  common  for  student  glider  pilots. 
These  results  do  not  show  the  strong  loads  that  are  possible  during  periods  of  strong 
turbulence  or  rotor  activity,  the  strong  gusts  which  occasionally  break  towropes  here 
at  the  Academy.  The  data  presented  herein  should  represent  a baseline  for  normal  tow- 
rope  loads,  and  further  testing  should  be  done  to  measure  towrope  loads  under  more 
extraordinary  flight  conditions. 
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APPENDIX 
AIRCRAFT  DATA 


Schweitzer  2-33A,  N 7760S 

gross  weight  (with  crew)  - W * 980  pounds  (Max:  1040) 

eg  position  - 94.10  Inch 

tow  speed  - 60  mph 

wing  surface  area  - S - 219.5  ft2 

Schweitzer  2-32,  N 2698H 

gross  weight  (with  crew)  - W ■ 1265  pounds  (Max:  1430) 

eg  position  - 112.34  inch 

tow  speed  - 70  mph 

wing  surface  area  - S ■ 180  ft8 

Piper  PA  18-180  Super  Cub,  180  HP 
N 3594P 

gross  weight  (with  crew)  - W • 1240  lb 

TOWROPE  Data  3/8  Inch  polypropylene,  2150  pound  breaking  strength 
\ inch  polypropylene,  1100  pound  breaking  strength 
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SECTION  II 
FLUID  MECHANICS 
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AN  INTUITIVE  LOOK  AT  VISCOUS  FLOW 


Eric  J,  Jumper 


Abstract 


This  paper  discusses  viscous  flow  based  on  some  simple  observations  of  a fictitious 
experiment.  The  experiment  Is  a contrived  apparatus  which  simulates  Impulsive  start 
motion  of  a flat  plate  In  an  at-rest  fluid  bath.  Based  on  the  observations,  predictions 
about  the  nature  of  boundary  layer  phenomena  are  made  and  compared  to  the  Exact  Blaslus 
solution  of  a boundary  layer  on  a flat  plate. 


Introduction 


There  are  many  who  have  derived  the  boundary  layer  equations  and  even  solved  and  used 
the  results  to  make  predictions;  but,  there  are  surprisingly  few  who  feel  that  they  have 
a good  Intuitive  feel  for  Just  what  Is  happening  In  a boundary  layer  and,  further,  Just 
what  to  expect  If  a given  set  of  conditions  for  which  the  behavior  Is  known  are  changed. 
In  what  follows  I will  describe  a set  of  simple  observations  upon  which  I feel  the  entire 
understanding  of  boundary  layer  theory  rests  (or  at  least  a large  part  of  it). 


In  order  to  Interpret  the  observations  we  will  make.  It  will  be  necessary  to  first 


understand  some  preliminaries 


A . Backgound 


If  you  want  to  limit  the  required  understanding  of  preliminaries  to  a bare  minimum, 
yqu  would  have  to  deal  with  the  Ideas  of  the  no-sllp  boundary  condition  and  of  shear 
stress  In  a fluid.  Of  course,  these  two  Ideas  open  up  an  entire  array  of  topics  which 
are  already  documented  far  more  completely  than  I could  ever  do  here  (for  example,  see 
references  1 and  2).  Let  us  start  by  realizing  that  we  will  be  treating  a fluid  as  a 
continuum  (and  a Newtonian  fluid).  In  doing  so,  we  are  able  to  refer  to  a rather  compll 
cated  two-phase  molecular  Interaction  taking  place  on  a solid/fluid  Interface  as  simply 
a no-sllp  condition.  In  the  usual  sense  of  a fluid  moving  past  a flat  plate,  this  would 
mean  the  portion  of  the  fluid  continuum  right  In  contact  with  the  surface  Is  not  moving 
with  respect  to  the  plate.  Of  course,  since  the  fluid  Is,  on  the  whole,  moving  with 


Captain,  USAF,  Associate  Professor  of  Aeronautics,  DFAN 


USAFA-TR-79-1 


I 

J 


respect  to  the  plate,  at  some  free  stream  velocity,  there  must  be  present  a region  where 
the  fluid  undergoes  a velocity  gradient  from  zero  at  the  wall  to  free  stream  velocity  at  the 
edge  of  the  region. 


This  region  In  which  the  effect  of  the  no-sllp  boundary  condition  Is  contained  is 
commonly  referred  to  as  the  boundary  layer.  It  is  defined  as  the  region  next  to  a solid 
boundary  where  the  fluid  goes  from  zero  velocity  at  the  wall  out  to  997.  of  the  freestream 
velocity.  The  distance  from  the  wall  to  the  99%  point  is  called  the  boundary  layer 
thickness,  4. 

Since  the  fluid  Is  "losing"  velocity  as  one  scans  from  the  boundary  layer  thickness 
toward  the  wall,  one  might  think  that  the  fluid  has,  in  effect,  a region  which  is  ex- 
periencing a momentum  defect,  momentum  being  the  mass  times  the  velocity.  Since  there  is 
a change  in  momentum  throughout  the  boundary  layer,  it  is  not  surprising  to  find  the 
presense  of  forces  and,  in  particular,  the  appearance  of  significant  shear  stresses 
(force  per  unit  area).  The  shear  stress,  t , In  the  fluid  is  defined  as 


(1) 


where  ^ is  the  coefficient  of  viscosity  of  the  fluid,  u the  velocity  and  y the  distance 
from  the  wall.  It  is  indeed  reassuring  to  find  that  because  the  shear  stress  is  directly 
proportional  to  the  velocity  gradient,  the  equation  would  predict  a shear  stress  in  the 
boundary  layer. 


Since  we  have  already  assumed  the  fluid  is  a continuum,  any  shear  stress  in  the 
bottom  of  a layer  of  fluid  will  be  balanced  by  the  top  of  the  layer  just  below  it.  This 
must  also  hold  right  at  the  wall.  The  shear  stress  in  the  fluid  right  at  the  wall  shows 
up  as  a tangential  force.  Thus,  the  shear  stress  at  the  wall 


Twall 


0 


(2) 
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hung  over  the  plate  so  that  It  can  be  repositioned  between  experiments,  but  once  positioned, 
remains  fixed  throughout  a given  run.  The  sensor  Is  able  to  detect  a velocity  disturbance, 
and  once  detected,  shut  off  a digital  timer  which  Is  Initially  activated  at  the  start  of 
an  experiment.  An  experiment  begins  when  the  plate  Is  set  Into  motion. 


B.  Taking  the  Data 


We  now  begin  to  take  the  data.  The  procedure  Is  to  select  a sensor  distance.  /. 

above  the  plate,  the  desired  plate  speed,  V , assure  the  fluid  Is  still  and,  flnsllv, 
throw  the  start  switch  which  jerks  the  plate  and  starts  the  timer.  The  timer  shuts  off 

ii  * 

when  the  velocity  disturbance  ratio,  /V^  , reaches  0.01,  where  u'  Is  the  velocity 
measured  by  the  sensor.  In  order  to  assure  that  Important  trends  are  noticed,  we  leave 
l untouched  over  a wide  range  of  Vp's  before  moving  the  sensor.  The  same  procedure, 
then,  Is  repeated  over  a correspondingly  wide  range  of  jt's.  Then  the  entire  procedure  Is 
repeated  with  a different  gas.  We  obtain  then,  an  abundance  of  data  as  shown  In  Table  1. 
Before  trying  to  Identify  trends  In  the  data,  it  Is  Important  to  realize  what  has  happened 
In  our  experiment.  Initially  the  fluid  was  at  rest,  but  because  of  our  no-sllp  boundary 
condition, we  caused  a sudden  velocity  disturbance  to  take  place  when  the  plate  was  Jerked, 
This  disturbance  propagated  through  the  fluid  and  eventually  triggered  the  sensor.  The 
disturbance’,  then,  was  In  effect  being  transported  away  from  the  plate.  Let  us  now  look 
at  how  this  transport  phenomena  was  effected  by  our  parameters. 


C.  What  the  Data  Tells  Us 


Probably  a good  place  to  start  Is  to  simply  look  at  the  time  It  takes  for  the 
disturbance  to  reach  the  sensor  at  a given  distance,  i,  for  gas  A.  For  the  purpose  of 
discussion  let's  assume  that  jtj  Is  our  closest  sampling  distance.  First  we  look  at  the 
time  tjj  which  Is  the  time  for  the  disturbance  to  travel  to  1 j for  a plate  speed  of 
V . Again  let's  assume  V Is  our  slowest  speed.  The  time  has  little  or  no  meaning  to 

Pi  Pi 

us  at  this  point  except  as  a reference.  One  thing  we  might  wonder  about  Is  how  this  time 
might  be  affected  by  Increasing  the  speed,  ; after  all.  It  do*.i  seem  that  If  were 
Increased,  a more  "violent"  disturbance  might  be  formed  and  our  logic  might  tell  us  that 
Che  time  to  travel  a distance  j ^ would  be  decreased.  In  examining  tla.  however,  we  find 
this  is  not  the  case.  In  fact,  we  find  that  all  the  times  connected  with  the  distance 


30 


USAFA-TR-79-1 


for  a given  gas  are  exactly  the  same. 


In  order  to  assure  ourselves  that  this  is  a rule,  we  examine  the  times  for  j tg  , £3  , 
etc.  and  find  that,  indeed,  the  time  doesn't  depend  on  the  speed  Vp  . We  do  find, 
however,  that  the  time  Increases  as  £ Increases;  the  relationship  is  a square,  that  is, 


t a £a 


(3) 


We  also  find  that  at  a given  £ , t changes  depending  on  the  type  gas  we  use.  Thus,  t is 
not  only  a property  of  the  distance,  £ , but  also  of  the  particular  fluid  medium  we  are 
using. 

In  examining  the  properties  of  the  fluids  and  comparing  the  times,  we  find  that  the 
times  for  a given  fluid  as  compared  to  another  fluid  differ  by  the  ratio  of  kinematic 
viscosities  , y , (more  often  given  as  y/p  where  p is  the  density)  so  that 


(4) 


where  the  subscripts  compare  gases  A and  B,  for  example,  and  the  relationship  holds  as 
long  as  “ ijj  • The  time,  then,  to  travel  a distance  £ now  beccxnes 


t a (5) 

V 

It  Is  Interesting  to  examine  the  meaning  of  the  kinematic  viscosity  , v » In  light 

of  the  above.  It  appears  that  y Is,  In  effect,  a transport  property,  transport  In  the 

sense  that  It  relates  to  the  transport  efficiency  of  a velocity  (or  momentum)  disturbance. 

The  units  of  y are  In  SI  units.  Comparing  the  units  of  the  kinematic  viscosity  with 

other  transport  properties  such  as  the  dlffuslvlty, in  mass  transport  and  the  thermal 

vP 

dlffuslvlty,  a , In  heat  transport,  we  see  that  they  all  have  the  same  units,  — . 

Perhaps  a better  name  for  1/  might  be  the  momentum  dlffuslvlty. 
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Returning  to  the  lesson  our  experiment  tells  us,  the  time  for  a velocity  (or  momentum) 
disturbance  to  travel  a distance  l Is  given  by 


(6) 


where  C Is  the  constant  of  proportionality.  We  must  now  ask  the  question  "how  does  this 
result  tell  us  anything  about  boundary  layers?" 

Ill,  What  About  Boundary  Lavers 

It  seems  to  me  that  we  now  have  the  essential  key  to  describing  what  goes  on  In  a 

boundary  layer.  In  order  to  use  It,  however,  we  should  first  formulate  a picture  of 

what  might  be  going  on.  Let' s look  at  the  formation  of  a boundary  layer  on  a flat  plate. 
A hypothetical  picture  of  what  might  be  going  on  Is  shown  schematically  In  Figure  2.  The 

key  is  that  we  are  looking  on  the  boundary  layer  as  simply  the  manifestation  of  a sudden 

velocity  (or  momentum) disturbance  (at  B)  being  propagated  out  away  from  the  plate  Into  an 
otherwise  undisturbed  fluid  (at  A). 


DISTURBANCE  BEING  CONVECTED 

■DOWNSTREAM  AT  U 

A B 

C 

^ . a.  ! 





- . 

^ 

^ 

(*'mS=L 

VELOCITY 

DISTURBANCE 

! DISTURBANCE 

PROPAGATING  AWAY 

FROM  PLATE 

Figure  2.  Boundary  Layer  Being  Looked  at  aa  tha  Propagation  of  a 
Disturbance 


33 


USAFA-TR-79-1 


The  only  difference  between  what  we  see  here  and  our  earlier  experiment  is  that  the  fluid 
is  being  convected  along  at  a velocity  U^.  If  this  picture  is  correct,  then  1)  we  must  be 
able  to  make  predictions  and  compare  them  favorably  against  the  exact  solution,  and  2)  if 
the  predictions  check  we  must  accept  the  fact  that  boundary  layers  are  not  as  complicated 
as  they  have  been  cracked  up  to  be.  First  the  predictions. 

A.  Predictions 

If  our  model  is  correct,  we  should  easily  be  able  to  predict  the  boundary  layer 

thickness  as  a function  of  distance  down  the  plate,  x . From  our  experiment  we  already 

know  that  the  rate  at  which  the  disturbance  propagates  away  from  the  plate  is  independent 

of  the  plate  velocity  or  In  our  new  case  (a  change  of  reference  frame)  the  fluid 

velocity,  U . Therefore,  we  should  be  able  to  expect  that  the  disturbance  will  arrive 
00 

at  a distance,  d , above  the  plate  In  time 


being  only  dependent  on  the  fluid  property,  v • Recall  that  at  this  point,  d , th^ 
disturbance  Is  a 1%  velocity  disturbance.  Recall,  also,  that  we  defined  the  boundary 
layer  thickness  as  the  distance  from  the  plate  to  the  1%  disturbance  point  (l.e.,  the 
velocity,  u , equals  99%  of  the  free  stream  velocity,  ).  The  starting  point  for  the 
disturbance  is  the  leading  edge  of  the  plate  (x  “ 0),  In  the  same  time,  t,  that  the 
disturbance  has  propagated  out  a distance, d.  It  has  also  been  convected  downstream  a 
distance,  £ , equal  to  U^t  (see  Figure  3). 

♦ 

I 

d 

I 

I 


Figure  3.  Disturbance  Being  Propagated  Out  From  the  Plate  and 
Convected  Downstream 
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By  definition, this  distance  Is  equal  to  the  boundary  layer  thickness,  fi,  corresponding 
to  the  d In  Eqn  (7).  So  that  In  general, 


1 


I 


CO 


(8) 


and 


or 


(9) 


j 


% 

(10) 


making  the  substitution  for  the  Reynolds  number  based  on  x , R , Eqn  (10)  can  be 

x 

written 


(11) 


This  is  our  first  prediction.  Notice  that  we  have  predicted  the  functional  form  of  the 
boundary  layer  thickness  Including  free  stream  velocity , distance  down  the  plate,  and  the 
appropriate  fluid  property;  It  would  be  more  than  coincidence  If  this  prediction  turned 
out  to  be  correct. 

Now  that  we  have  a formula  for  predicting  the  boundary  layer  thickness,  we  can  make 
some  predictions  about  what  happens  to  the  boundary  layer  as  we  change  the  free  stream 
velocity.  From  Eqn  (10)  we  can  see  that  as  the  velocity  Increases,  the  boundary  layer 
thickness  at  the  same  location  decreases.  Intuitively,  then,  the  rate  at  which  the 
disturbance  propagates  away  from  the  plate  Is  unchanged,  but  In  the  same  time  It  has  been 
convected  farther  downstream.  The  slower  the  freestream  velocity,  then,  the  thicker  the 
boundary  layer. 

Even  though  we  don't  know  the  exact  shape  of  the  velocity  profile  In  the  boundary 
layer,  we  can  also  make  predictions  about  the  shear  stress  at  the  wall  (drag  per  unit 
area)  using  the  following  argument.  We  know  the  shear  stress  at  the  wall  Is  given  by 
Eqn  (2);  we  are,  therefore,  looking  for  the  slope  of  the  velocity  profile  right  at  the 
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wall.  We  also  know  that  the  velocity  must  go  from  zero  to  U over  the  distance  A . As 

00  u 

a quick  approximation,  let ' s assume  the  change  Is  linear,  thus 


U 


ay'y-o  6 


(12) 


at  a minimum  we  know  this  will  give  us  a ball  park  estimate.  Thus,  the  gradient  of  the 
velocity  Is  lnversly  proportional  to  the  boundary  layer  thickness.  As  the  velocity  In 
the  free  stream  Is  Increased,  the  boundary  layer  thickness  decreases  and  the  gradient  must 
Increase.  So,  as  the  velocity  Increases,  the  shear  stress  (thus,  the  drag)  Increases. 

Our  predictions  are  summarized  In  Table  2. 

Table  2 

Summary  of  Predictions  Based  on  Intuitive  Model 


Change 

In 

U 

30 

Effect  on 

6 

Tu 

Increase 

decrease 

Increase 

decrease 

Increase 

decrease 

functional 
form  of  6 

e 

x 


V 


B.  Checking  our  Predictions 


Fortunately,  the  exact  solution  for  the  boundary  layer  on  a flat  plate  Is  well  known. 
The  first  solution  of  this  problem  (Incompressible)  was  obtained  by  H.  Blaslus  In  1908 
(Ref  4)  and  Is  thus  known  as  the  Blaslus  Solution.  The  governing  equations  are,  of 
course,  the  continuity  and  momentum  equations, 


continuity: 


momentum: 


3x  9y 


au  . au  a 8 u 

U rr  + V “r  - i;  * 


dx 


ay 


ay8 


(13) 


(14) 
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where  u Is  the  velocity  In  the  x direction  and  v In  the  y.  The  solution  Is  given  In 
Figure  4 (Ref  5).  We  can  see  from  Figure  4 that  the  IX  disturbance  point  (U/U^  “ .99) 
occurs  at  7J  ■ 5.0  . At  this  point  y must  equal  0 . Thus, 


5.0 


(15) 


or 


0 " 5.0 


5.0  x 


(16) 


Figure  4.  Blaslus  Solution  of  the  Boundary  Layer  for  Flow  Over 
a Flat  Plate 


We  see  that  our  prediction  of  the  functional  form  of  0 as  given  by  Eqn  (10)  and  (11)  is 
exactly  correct.  Since  the  equations  are  the  same,  our  predictions  of  the  effect  of 
boundary  layer  thickness  on  velocity  are  also  exactly  correct.  Note,  also,  that  included 
In  Figure  4 la  our  linear  approximation  to  the  velocity  profile  which  does  give  a ball 
park  feel  for  the  shear  stress^  So,  we  find  that  all  our  predictions  ere  exactly 
correct. 
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IV.  Conclusions 

We  have,  over  the  course  of  these  few  pages,  built  what  I feel  Is  a very  good 
Intuitive  model  of  what  Is  happening  In  a boundary  layer.  The  Important  Idea  of  dis- 
turbance transport  Is  the  foundation  of  this  model  and  should  be  the  key  on  which  you 
continue  to  build  your  Intuition. 

I should  also  state  that  we  have  dealt  here  only  with  laminar  boundary  layers  In 
which  disturbances  are  propagated  via  molecular  Interactions.  The  Ideas  here,  however, 
can  be  extended  to  turbulent  boundary  layers  by  simply  realizing  that  our  disturbance  Is 
still  being  propagated  or  transported  away  from  the  surface  but  by  different  mechanisms, 
!•*.,  eddy  motion  of  the  flow.  The  rate  of  transport,  then,  Is  no  longer  simply  a 
property  of  the  flow  condition.  While  things  are  somewhat  more  complicated,  the  same 
simple  Intuitive  approach  Is  valid;  don't  take  my  word  for  It,  check  some  solutions. 

My  last  remark  Is  this  : Intuition  In  engineering  la  essential.  No  matter  how 
complicated  analysis  becomes,  there  are  always  a few  basic  physics  principles  which 
hold  the  key  to  understanding  the  phenomena.  Your  job  as  a good  engineer  Is  to  seek 
these  principles  and  develop  good  Intuition  In  your  field.  Remember,  a computer  will 
always  let  you  down  when  you  need  It  most,  but  your  Intuition  Is  a constant  reliable 
friend. 
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FLOW  FIELD  CHARACTERIZATION  OF  FLOW  DIRECTION  AND  MACH 
NUMBER  PROBES  USING  SLENDER  BODY 
THEORY 

* 

G.  David  Huffman 

Abstract 

Slender  body  theory  has  been  used  to  model  the  flow  field  past  a flow  direction  and 
Mach  number  probe  at  an  angle  of  attack  and  a sideslip  angle.  An  expression  has  been 
derived  for  the  body  surface  pressure  coefficient.  This  has  been  converted  Into  ang- 
ular and  static  pressure  coefficients  which  are  evaluated  for  a specific  sensor  geom- 
etry. These  expressions  are  employed  In  the  quantification  of  misalignment  effects. 
Finally,  the  velocity  expressions  are  used  to  visualize  the  three-dimensional  flow  past 
a cone  probe  at  an  angle  of  attack  and  a sideslip  angle. 

I.  Background 

Flow  fields  encountered  in  engineering  applications  are  often  three-dimensional  and, 
thus,  have  both  a varying  flow  direction  and  velocity /Mach  number.  Turbomachine  and 
aircraft  environments  are  often  harsh,  necessitating  a rugged,  durable  sensor.  This 
requirement  makes  multi-ported  pressure  probes  particularity  attractive  for  these 
applications. 

Pressure  probes  frequently  consist  of  an  aerodynamic  body  with  a symmetrical 
arrangement  of  sensing  holes.  A number  of  different  geometries  have  been  investigated 
and  some  typical  cases  are  reviewed  In  references  1 through  15.  Probes  are  treated  in 
more  general  terms  in  references  16  and  17. 

Pressure  probes  are  normally  employed  in  either  the  stationary  or  nulling  moda.  In 
the  nulling  or  equal  pressure  approach,  the  probe  Is  oriented  such  that  each  of  the 
side  ports,  see  Figure  1,  reads  the  same  pressure.  The  probe  position  is  noted  and  the 
flow  direction  determined.  in  the  stationary  mode , the  probe  Is  fixed  and  the  top- 
to-bottom  and  side-to-slde  pressure  differences  measured.  Calibration  functions  are 
then  used  to  find  the  flow  angles  a and  fi.  The  static  and  total  prassuras  and,  thus, 
Mach  number  can  be  determined  In  a similar  manner.  Both  these  methods  have  merit; 
however,  turbomachine  and  aircraft  applications  normally  require  an  ionovable  sensor. 

As  a consequence,  only  stationary  probes  will  be  considered. 

Since  the  probes  are  to  be  used  without  rotation,  tha  sensitivity  to  flow  angularity 
Is  axtremely  Important.  While  values  of  this  parameter  are  available  for  a number  of 
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Figure  1.  Typical  Flow-Direction  Probe  (Each  Unit  on  Scale  Represents  0.5  mm) 


calibrated  probes,  this  data  is  of  United  value  in  synthesizing  a sensor  geometry  for 
a specific  application.  Therefore,  a general  analytic  model  has  been  developed  using 
slender  body  theory  and  this  is  discussed  in  Section  II. 

In  addition  to  the  probe  sensitivity  which  is  largely  dictated  by  aerodynamic  con- 
siderations, alignment  and  manufacturing  defects  also  Influence  the  accuracy  with  which 
flew  angles  and  static  and  total  presaurea  can  be  determined.  Effects  of  this  nature 
are  discussed  In  Section  III.  The  mathematical  model  is  used  to  visualize  the  flow 
field  In  Section  IV.  The  results  of  this  investigation  are  summarized  in  Section  V. 

II.  A Mathematical  Model  of  Probe  Aerodynamic  Behavior 

A.  Introduction 

The  objective  of  an  aerodynamic  probe  is  to  determine  the  magnitude  and  direction 
of  the  velocity  vector.  This  translates  into  a measurement  of  pressures  which-by  means 
of  calibration  functlons-are  then  converted  Into  flow  angles  and  total  and  static 
preasures.  The  flow  field  parameters  should  be  accurately  measured  with  the  probe 
itself  creating  a minimal  disturbance. 
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A typical  flow-direction  probe  having  an  aerodynamic  shape  can  be  thought  of  as  a 
slender  body  (the  body  radius  is  much  less  than  the  body  length.)  The  use  of  slender 
body  theory  in  the  analysis  of  probes,  thus.  Immediately  comes  to  mind.  There  are, 
however,  shortcomings  associated  with  this  approach,  the  major  one  being  that  the  rate 
of  change  of  body  radius  with  respect  to  body  length  must  be  small.  This  feature  pre- 
cludes stagnation  regions  in  a slender  body  approach. 

Of  course,  a complete  description  of  the  flow  field  around  an  aerodynamic  probe 
can  be  generated  by  numerically  solving  the  three-dimensional  potential  flow  equations, 
e.g.,  references  18  and  19.  While  this  approach  is  more  accurate  than  slender  body 
theory  under  some  conditions,  it  does  not  lend  itself  to  synthesis  of  probe  shapes. 
Furthermore,  analytic  calibration  relations  are  not  obtained,  and  considerable  insight 
into  the  physical  process  is  lost.  Consequently,  a slender  body  approach  will  be 
employed  in  the  ensuing  analysis. 

B.  Basic  Formulation  of  Slender  Body  Theory 

A slender  body  of  revolution  in  a cross-flow  is  shown  in  Figure  2.  The  body  Itself 
can  be  described  in  terms  of  an  axial,  z,  and  radial,  r,  coordinate  with  r - R(z)  on 
the  body  surface. 
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Following  Shapiro,  reference  20,  and  Llepmann  and  Roshko,  reference  21,  Che  perturbation 
potential  can  be  written  as 


«Vr  + r + r*  % + (1  * M»*>  " 0 


(1) 


where  the  subscripts  r,  0 and  z denote  partial  differentiation.  Equation  (1)  Is  de- 
rived by  assuming  that  the  potential  function  $ can  be  written  as  the  sum  of  the  per- 
turbation due  to  the  body,  <p,  and  the  potential  due  to  the  external  flow,  Vm.  It  is 
further  assumed  that  the  perturbation  velocities  are  much  smaller  than  the  free-stream 
velocity.  The  perturbation  velocities  are  related  to  the  potentials  via 


Vr  m v<r>,  ^ ■ rv(S),  <fiT  - v(r) 

The  boundary  condition  on  the  body  surface  can  be  written  as 

grad  $ • grad  F ■ 0 on  F(r,  fi,  *)  ■ 0 

where  F(r,  0,  *)  - r - R(s)  and  grad  4>  - grad  <p  • \T,  * lT. 

Evaluation  of  equation  (3)  yields 

[v0B(r)  + v(r>]  - [vj*>  + v(r)j  R'  ■ 0 on  F(r,  0,  *)  - 0 

with  R ' ■ . The  perturbation  must  vanish  as  r approaches  infinity  and 

v(r>,  v(0),  v(8)  - 0 as  r - » 


(2) 


(3) 


(4) 


(5) 


Equations  (1),  (4)  and  (3)  comprise  the  basic  relations  for  the  perturbation  flow 
around  a slender  body.  The  solution  of  equation  (1),  subject  to  Che  appropriate 
boundary  conditions,  Is  discussed  In  the  following  section. 

C.  Solution  of  the  Partial  Differential  Equations 

The  system  of  equations  of  Section  II. B Is  solved  by  dividing  the  flow  field  Into 
two  components:  an  axlsynsatrlc  flow  past  the  sensor  and  a lateral  flow  past  the 
same  body.  The  two  solutions  are  then  summed  and  the  complete  flow  field  Is,  thus, 
modelled.  The  two  conditions  can  be  expressed  mathematically  as 


*rr(->  + 7 V>r<*)  + <l  " **»*  > " 0 

V>r(4)  - V.(,)  R on  F(r,  0,  «)  - 0 

- 0 as  r - * 


(6) 
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<pJL)  - - V„(r'  on  F(r,fl,s)  - 0 


*,< *>.  . - 0 , 


where  Che  superscript  a denotes  the  axlsynmetric  and  l the  lateral  component.  It  then 
follows  that  <p  » + V>(i>. 

Equation  (6)  can  be  written  In  terms  of  a source  distribution,  q(z),  following 
Sears,  reference  22,  as 


9>(a)  - - kvf* 


/(z  - C)“  + 62 


\ 

h 


where  6 » o/l  - M^2  . The  velocity  components  are  then  given  as 


<«>  . — L_  fl  («  - C)  1 (O  « C 3,2 

' [(.-O' 


,<r>  - <p  (“>  _ Se  /“ 

4IT  / 
to 


[(*  - O*  + 6*  r*  ]' 


At  this  point  the  source  distribution  function  Is  unspecified.  In  essence,  the 
specification  of  q produces  an  arbitrary  body.  Since  the  body  Is  generally  known,  a 
direct  determination  of  q (()  In  terms  of  R (z)  Is  required.  An  approximate  solution 
of  this  type  can  be  obtained  by  writing 


q(C)  - q<z  + 6 r ?J>  “ E q(n)  (*) 

n-0  nI 


where  qv;  (s)  denotes  the  n-th  derivative  with  respect  to  s.  Substitution  of  equation 
(11)  into  equations  (9)  and  (10)  and  term-by-term  integration  yields 

*.w  - • Th  1^  [if f- ] - -'<*>  • •]  * ...  I <«> 


P 
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<P 


(a)  „ 1_  I 2 q(z) 


4 n I 6 r 


for  terms  of  order  less  than  fir  with  6r<  <1  -z.  Application  of  the  boundary  conditions 
on  the  body  show  that 

q(z)  - 6 vjz)  s' 

where  S denotes  the  body  cross-sectional  area,  ff[R(z)]a  . 


The  lateral  potential  function  can  be  determined  in  a similar  manner  by  noting  that 
if  (r,  z)  is  a solution  of  equation  (6),  then  both  sin  6 <p  and  cos  9 are  solutions 
of  equation  (7).  Again  using  a series  expansion  for  the  source  distribution,  term-by- 
term  Integration  and  application  of  the  boundary  conditions  yields 


<J> 


V_(x)  COB  9 + V-(y)  sin  8 


■ 


(y) 


cos  8 - V„(x)  sin  9 


e.) 


(!) 


U)  V„(x)  cos  0 + V.(y)  sin  0 . 

! ir  lr  / 


(15) 


(16) 


(17) 


The  potential  functions  can  be  combined  and  the  velocity  components  and  pressure 
coefficient  determined. 

D.  The  Velocity  Components  and  Pressure  Coefficient 

Equations  (12)  and  (17)  and  (13)  and  (15)  can  be  combined  to  yield  the  pertur- 
bation velocities  and 


(r) 


(r !') 


\ cos  a cos  P - (sin  P cos  0 + Bin  a cos  P sin  0r^- 


„(8)  „ 

— » (sin  a cos  fi  cos  0 - sin  fi  sin  0)  ~ 


(*) 


~ % cos  a cos  fi  f + (sin  fi  cos  0 + sin  a cos  fi  sin  0) 


(18) 


(19) 


(20) 
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where  f denotes  the  body-geometry  function  and  is 


- - h (r8)'(  / • " I ' ) + \ (R*)"  f j — 

z®  + 6*  r®  V (-2-z)®  + 6®  r®  ' L J (Jl-z)®  + 6®  rs 


z \1  - t a*)"'/  -** 

+ 6®  r®  ' A-z  + J" (&- z)®  F 6®  r®  ’ ^ ' '/  z*  + 6® 


«-z)® 


(f-z)®  + 6®  r® 


1-2  V (i-z)®  + 6*  r®  - 2 V z®  + 6*  r 


2 * A*  -2 


| + • • • 


Equation  (21)  is  not  subject  to  the  constraint  that  fir  < < f-z. 

In  probe  analysis,  the  pressure  coefficient  is  of  major  importance.  The  pressure 
itself  can  be  related  to  the  velocity  by  means  of 


P “ Poo  - ^ PV« 


V„  . v 


where  v is  evaluated  on  the  body  surface.  The  pressure  coefficient  is  defined  as 

P-P®  V®  • 7 a 


This  can  be  written  as 

Cp  * cos*  a cos*/J  [ -f  -(r')*  ] + sin®  S [ 1 - 4 sir£  0]  + 
sin*  a cos*  f)  [ 1 - 4 cos®  0 j + cos  a sin  2 [ - 2 R1  cos  0]  + 

sin  2 a cos*  jS  [ - 2 R*  sin  0]  + sin  2 a cos  fi  [ 4 sin  0 cos  0]  (24 

where  f is  to  be  evaluated  on  the  body  surface,  l.e.,  r ■ R.  Equation  (24)  represents 
the  major  result  of  the  analysis  and  can  be  used  to  determine  probe  angular  and  static 
pressure  sensitivity  for  arbitrary  geometries. 
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III.  Aerodynamic  Probe  Angular  and  Static  Pressure  Sensitivity 

A.  Differential  and  Averaged  Pressure  Coefficients 

The  flow  angularity  Is  determined  by  differencing  the  measured  pressures  between 
the  side  and  top  and  bottom  ports,  see  Figure  1, for  instance.  Since  these  ports  are 
at  the  same  z and  R locations,  the  pressure  difference  Is  generated  by  subtracting  Cp 
values  at  different  8 locations  and 

Ac  - - 4 sin*  JJ  (sin*  9 > - 4 sin*  a cos*  B (cos*  9 cos*  9'  1 - 
P l Si 

2 R'  cos  a sin  2 B (cos  8 - cos  8 ) - 2 R1  sin  2 a cos*>9  (sln8  - sin  8 ) + 

2 1 2 l 


4 sin  2 a cos  B (sin  8 cos  8 - sin  9 cos  8 ) 

Si  S ii 


(25) 


The  sensitivity  to  changes  In  flow  angle  is  defined  as  6 ACp/0 a or  b Acp/b0.  The 
former  can  be  written  as 

b Ac 


b a 


^ ■ 8 R1  cos  2 a cos*  B 


(26) 


where  8 and  8 have  been  taken  as  90  and  270  respectively.  Note  that  the  sensitivity 

it 

to  flow  direction  in  one  plane  depends  to  some  extent  on  the  flow  angle  in  the  other 
direction. 

Equation  (26)  is  plotted  In  Figure  3.  It  is  quite  apparent  from  the  figure  that  the 
probe  angular  sensitivity  — regardless  of  shape  — depends  on  a and  /9.  Approximately  a 
107.  reduction  In  sensitivity  occurs  for  O and  B values  of  10°.  Note  that  b A Cp/b  a Is 
the  same  for  both  positive  and  negative  values  of  a and  B- 

The  averaged  pressure  coefficient  can  be  used  to  determine  the  flow>field  static  pres- 
•ur«»  pro*  cp  avaluated  at  a number  of  theta  values,  the  results  summed  and  divided 

by  the  number  of  theta  values.  Mathematically  this  process  yields 

4 

< Cp>  - cos*  a cos*  fi  [-f  • (R')*  ] + sin*  fi  (1  - £ sin*  0 + sin  a cos */3  (1  - 

4 4 4 

£ cos2  0 ) - h R'  cos  a sin  2 fi  £ cos  0.  - h R'  sin  2 a cos*  fi  £ ,ln  0<  + 

1-1  \ 1-1  1-1 


sin  2 a cos  0E  sin  0.  cos  0. 

1-1  1 1 


(27) 


With  0t  as  0,  90°,  180°  and  270°,  equation  (27)  becosMS 

46 


k 

I * 


1 


USAFA-TR-79-1 


<C  > 

P 


a cos*  0 [-f  - (R')«  ] - sin*  a cos  * 0 


which  is  approximately  parabolic  with  respect  to  q and  0. 

The  analysis  is  restricted  to  small  perturbations  velocities  and,  thus,  a total 
pressure  coefficient,  i.e.,  velocity  stagnated,  cannot  be  directly  derived.  A quasi- 
total pressure  can  be  formulated  by  integrating  the  pressure  over  the  body  surface  and 
resolving  this  force  into  an  axial  component.  Dividing  by  the  cross-sectional  area 
then  yields  an  axial  pressure  which  is  of  the  form  of  a total  pressure.  The  functional 
form  of  this  relationship  is 


<C  > 

p z 


J (R*)'  [-f  - (R')*]dz  - sin*j3  - sin*o?  cos*  0 


where  l denotes  the  integration  length. 


B.  Pressure  Coefficients  for  a Typical  Probe 


Equations  (25)  and  (27)  have  been  evaluated  for  a typical  conical  pressure  probe. 

The  probe  geometry  is  shown  in  Figure  1 with  an  included  conic  angle  of  30°,  a maximum 
diameter  of  3.175mm  (0.125  in),  a length  of  4.740mm  (0.1866  in)  and  the  probe  side 
ports  located  3.073mm  (0.121  in)  from  the  probe  tip.  The  side  port  holes  are  0.452mm 
(0.0176  In)  in  diameter. 

The  pressure  coefficients  are  plotted  in  Figure  4a  in  terms  of  (A  C ) and  ( &C  ) 

P a P j3’ 

with  q and  A chosen  as  either  90°  and  270°  or  0°  and  160°.  The  flow  angles  <y  and  8 

1 8 

along  with  the  pressure  port  locations  in  terms  of  are  specified  on  the  figure.  The 
averaged  pressure  coefficient,  <C^>,  is  shown  In  Figure  4b  for  the  same  probe. 

These  calculations  were  carried  out  at  a value  of  0.2.  The  theory  has  been  veri- 
fied for  Mach  numbers  up  to  0.980  in  Ref.  23. 

C.  Effect  of  Probe  and/or  Side  Port  Alignment  on  the  Pressure  Coefficients 

The  Influence  of  probe  and/or  side  part  alignment  on  the  angular  and  averaged 
pressure  coefficients  can  be  calculated  from  equations  (25)  and  (27).  The  port  locations 
denoted  by  ff  can  be  rotated  relative  to  the  design  locations  of  0°,  90°,  180°  and  270°, 
and  the  resulting  coefficients  calculated.  Coefficients  generated  in  this  manner  are 
plotted  in  Figures  5 and  6.  Note  that  the  cone  of  Section  III.B  has  been  used  in  these 
calculations. 
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Figure  4a.  (AC  ) Versus  (AC  for  a Typical  Flow  Direction  Probe  (6.  • 0°,  90°, 
opa  o pP  1 

180  , and  270  ) 
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Figure  6a.  (AC  ^ Versus  (AC  )o  for  a Typical  Flow  Direction  Probe 
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Misalignments  are  of  two  types.  The  first  consists  of  a rotation  of  the  probe  body 
Itself  relative  to  a fixed  or  pre -determined  coordinate  system.  This  could  occur  in  a 
turbomachine  application  during  installation  of  the  sensor.  In  this  case  all  f)  values 
are  displaced  by  the  same  amount.  The  angular  and  averaged  coefficients  are  plotted 

in  Figure  5 for  a port  displacement  of  2.5°.  Note  that  (£C  ) and  (A  C ) are 

P a P 0 


strongly  influenced  by  the  rotation  but  <£p>  1®  unchanged. 


The  second  misalignment  consists  of  a manufacturing  defect  with  only  one  port 
displaced.  This  case  is  plotted  in  Figure  6 for  one  side  port  located  at  2.5°  rather 
than  0.  Now  both  the  angular  and  averaged  pressure  coefficients  are  Influenced. 


IV.  Flow  Field  Visualization  Using  Slender  Body  Theory 

The  flow  field  past  a slender  body  at  both  an  angle  of  attack  and  sideslip  angle 
can  be  visualized  by  generating  a series  of  particle  trajectories.  The  velocity  com- 
ponents can  be  Integrated  yielding  the  Instantaneous  particle  position,  i.e., 

r - ro  + f V(r)  dt  , * flo  +/  r * >dt  (30) 

0 

t 

z - ZQ  + y V(r)  dt 
0 


where  V^  ■ + v^r\  etc.  The  flow  field  past  the  cone  probe  of  section  III.B  at 

a and  8 values  of  7°  and  11°  is  shown  in  Figure  7. 

V.  Sutanary 

Slender  body  theory  has  been  used  to  model  the  flow  field  past  an  arbitrary  body 

of  revolution  at  a sideslip  angle  and  angle  of  attack.  Expressions  for  the  pressure 
coefficient  on  the  sensor  surface  have  been  derived.  From  these  relations,  angular 
and  averaged  (static)  pressure  coefficients  hsve  been  determined.  Expression 
for  the  angular  probe  sensitivity  have  also  been  evaluated.  Finally,  the  effects  of 
probe  and  probe  side  port  misalignments  have  been  quantified  for  a given  sensor 
geometry. 
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Symbols 


Pressure  Coefficient 


<cp> 


Angular  Pressure  Coefficient 

Averaged  Pressure  Coefficient 

Body  Geometry  Function 
length 

Free-stream  Mach  Number 


Pressure 


Source  Distribution  per  Unit  Length 
Cylindrical  Coordinates 
Body  Radius 

Body  Cross-Sectional  Area 
Perturbation  Velocity 
Velocity 

Cartesian  Coordinates 
Angle  of  Attack 
Sideslip  Angle 

J 1 - M„2 

Density 

Perturbation  Potential 
Potential  Function 


Subscripts 

Partial  Differentiation  with  Respect  to  the  Cylindrical  Coordinates 

Initial  Condition 

Free-stream 


Superscripts 


Transverse  or  lateral 
Cylindrical  Coordinates 
Cartesian  Coordinates 


57 


USAFA-TR-79-1 


1. 

2. 
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WAVE  DRAG  PREDICTIONS  ON  SLENDER  BODIES  OF 
REVOLUTION  USING  THE  SUPERSONIC  AREA  RULE 

: k 

Glen  R.  Schlotterbeck 
Abstract 


This  report  discusses  the  mechanics  and  results  of  a computer  program  designed  to 
compute  the  wave  drag  of  slender  bodies  of  revolution  In  supersonic  flow.  The  program  Is 
based  on  the  supersonic  area  rule  and  Is  Intended  for  use  as  a tool  In  the  first  cut  of 
the  design  process.  Hopefully,  the  program  will  Indicate  the  most  advantageous  arrange- 
ment of  vehicle  components  for  a minimization  of  wave  drag  over  a range  of  Mach  numbers. 


I.  Introduction 


A.  Background 


Von  Karman,  In  working  with  source  distributions,  first  made  the  observation  that  the 
wave  drag  Integral  was  analogous  to  the  classical,  Incompressible  equation  for  the  Induced 


drag  of  a lifting  wing  of  length  l (Ref.  1).  Specifically,  von  Karman  noted  that 


D 

w 


d £ V (x)  f' (£)  log  j x 
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(1) 


where  f (x)  was  the  Intensity  of  the  source  distribution  and  p was  the  density  of  air. 
If  the  function  f (x)  were  changed  to  denote  the  spanwlse  circulation  distribution  on  a 
loaded  line,  the  induced  drag  would  be  given  by 


D 


1 


dx  d C f'(x)  f'<C>  log 
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(2) 


which  Is  Identical  to  the  wave  drag  equation. 

Based  on  this  observation  and  slender  body  theory,  Whitcomb  developed  what  Is  now 
known  as  the  "area  rule"  (Ref.  2).  Using  stream  tube  observations  and  choking  character- 
istics, Whitcomb  realized  that  three-dimensional  flow  fields  within  each  Individual 
stream  tube  must  obey  one-dlmenslonal  flow  laws.  His  "area  rule,"  good  In  the  transonic 

* 
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region,  stated  that  a slender  body  contained  within  a flow  near  the  speed  of  sound  had 
the  same  drag  as  a body  of  revolution  with  the  Identical  cross-sectional  areas  of  planes 
cut  normal  to  the  longitudinal  axis. 

Using  the  relationship  between  three-dimensional  flow  and  one-dimensional  flow,  along 
with  supersonic  influence  characteristics,  Jones  expanded  upon  Whitcomb's  earlier  work  to 
Include  the  supersonic  regime  (Ref.  3).  His  expansion  was  based  on  the  Idea  that  the 

i 

three-dimensional  disturbance  field  could  be  modeled  by  the  superposition  of  one- 
dimensional disturbances  taking  the  form  of  plane  waves.  Instead  of  using  cross-sectional 
areas  cut  by  planes  normal  to  the  longitudinal  axis,  Jones'  method  used  areas  formed  by 
planes  inclined  at  the  Mach  angle,  called  ''Mach  planes"  (see  Figure  1). 


Figure  1.  Mach  Plane  Areas  (Ref.  4) 


Experimental  observation  led  Nelson  and  Welsh  to  conclude  that  greater  accuracy  could 


be  obtained  by  using  the  frontal  projection  of  these  obliquely  cut  Mach  plane  areas 


rather  than  the  actual  cut  areas  (Ref.  5) 


Thus,  by  redefining  the  area  distribution  necessary  for  the  von  Karman  wave  drag 


equation,  the  area  rule,  previously  valid  only  In  the  transonic  region,  was  expanded  to 


Include  the  supersonic  regime 
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B.  Scope 


This  report  describes  a program  which  uses  the  von  Karman  formula  coupled  with  the 
Mach  plane  area  distribution  concept  advanced  by  Jones,  together  with  the  frontal  projec- 
tion technique  developed  by  Nelson  and  Welsh,  to  compute  the  wave  drag  of  a given  slender 
body  of  revolution.  The  program's  approach  Is  to  1)  calculate  the  frontal  projection  of 
the  areas  formed  by  Mach  plane  cuts  at  regular  Intervals  along  the  body  axis  for  a given 
Mach  number,  2)  calculate  the  second  derivative  of  this  area  distribution,  and  3)  compute 
the  wave  drag  coefficient  based  on  the  wave  drag  formula  presented  by  von  Karman. 


II.  Theory 


Based  on  slender  body  and  linear  theory,  the  formulation  of  von  Karman' s wave  drag 
equation  begins  with  the  Idea  that  a distribution  of  sources  located  along  an  axial  line 
describes  a body  of  revolution.  It  uses  power  series  expansions  to  arrive  at  the  form 
given  In  Eqn  (1).  Reference  1 offers  the  most  logical  and  comprehensive  development  of 
the  wave  drag  equation. 

In  order  to  apply  Eqn  (1)  to  predict  the  wave  drag  in  the  supersonic  region,  the 
distribution  of  the  frontally  projected  Mach  plane  areas  for  the  body  must  be  calculated. 
The  corresponding  distribution  of  second  derivative  values  of  these  areas  must  then  be 
determined  for  direct  substitution  Into  the  wave  drag  equation. 

For  practical  applications  the  non-dimensional  coefficient  form  of  the  wave  drag 
equation  Is  most  often  used.  However,  as  presented  by  von  Karman,  Eqn  (1),  transformed 
Into  coefficient  form,  Is  still  not  non-dimens  Iona  Used.  Specifically,  units  of  length 
remain  In  the  argument  of  the  logarithm  term.  To  non-dlmenslonallze  this  term,  therefore, 
both  x and  g need  to  be  divided  by  a length  value,  namely,  the  total  length  of  the  body 
(Ref.  6).  This  transforms  Eqn  (1)  Into  the  truly  non-dimensional  coefficient  form  of  the 
wave  drag  equation,  given  by 
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where  SRef  Is  normally  taken  as  the  maximum  cross-sectional  area  of  the  body.  Eqn  (3) 
appears  to  be  the  only  working  form  of  the  wave  drag  equation  In  truly  non-dimensional 
terms  and  has  not  been  found  In  any  other  presentation  of  the  formula. 

Numerical  Integration  of  Eqn  (3)  Is  the  logical  solution  approach  due  to  the  finite 
number  of  second  derivative  values  calculated.  However,  In  the  region  where  £ approaches 
x the  value  of  the  natural  logarithm  term  goes  to  Infinity  and  Is  actually  undefined 
where  g equals  x . This  necessitated  the  use  of  analytical  methods  employing  a limiting 
technique  In  the  region  where  g and  x were  one  unit  apart  (Ref.  7).  For  this  reason, 

Eqn  (3)  was  separated  Into  three  regions,  two  of  which  were  numerically  evaluated  and  the 
third  region,  where  g approaches  x , which  was  analytically  evaluated.  Eqn  (4)  gives 
this  expanded  form  of  the  wave  drag  equation  In  coefficient  form  as  used  In  this  computer 
program 
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III.  Discussion 

The  program  divides  the  calculation  of  the  wave  drag  coefficient  Into  three  logical 
phases.  The  first  phase  deals  with  the  generation  of  a frontally  projected  Mach  plane 
area  distribution  describing  the  body  In  queatlon.  A maximum  of  one  hundred  axial  loca-s 
tlons  were  Investigated  with  respect  to  Mach  plane/body  Intersection  although  more  could 
be  added  If  required.  One  hundred  eighty  radial  distances,  from  body  axis  to  Mach  plane/ 
body  surface  Intersection,  were  calculated  In  order  to  describe  the  Mach  plane  area  (see 
Figure  2). 
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Figure  2.  Radial  Distances  of  Mach  Plane  Area 

The  bodies  Investigated  to  date  have  been  tapered  and  since  the  Mach  planes  ate 
Inclined,  the  Mach  plane  areas  were  egg  shaped  rather  than  circular.  This  required  the 
calculation  of  more  than  one  radius  value  for  Mach' plane  area  sizing.  Due  to  symmetry 
about  the  Y-Y  axis,  as  seen  In  Figure  2,180  radius  values,  corresponding  to  180  degrees 
revolution  about  the  body  axis,  adequately  described  each  Mach  plane  area.  From  the 
radius  values  obtained,  the  program  calculates  the  total  Mach  plane  area  for  each  axial 
location,  generating  an  area  distribution  of  the  body. 

The  second  phase  deals  with  generating  the  second  derivative  distribution  from  the 
area  distribution  of  the  first  phase.  Initially,  a cubic  spline  fit  was  attempted  for 
the  first  and  second  derivation  calculations,  resulting  In  extreme  fluctuations  In  the 
second  derivative  distribution.  For  this  reason  a five-point  nested  averaging  technique 
Is  now  used.  Starting  with  the  area  distribution,  a range  of  four  slopes  surrounding  a 
given  axis  location  are  averaged  to  arrive  at  a first  derivative  value  for  that  point. 

The  same  procedure  Is  again  used  to  calculate  the  second  derivative.  Each  axial  location 
along  the  body  Is  handled  In  this  manner  except  the  extreme  ends  of  the  body  where  a flve- 
polnt  range  Is  not  available.  In  these  two  regions  approaching  either  end  of  the  body, 
progressively  smaller  slope  averaging  ranges  are  used.  The  end  result  Is  the  generation 
of  the  second  derivative  of  the  area  distribution  for  the  body  under  Investigation. 
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Figure  3.  Slender  Body  (not  to  scale) 
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Th®  third  and  final  phase  of  the  program  deals  with  actually  evaluating  the  wave  drag 
equation  as  presented  In  Eqn  (4).  Knowing  the  reference  area,  S , of  the  body  and  the 

R6l 

necessary  area  second  derivative  values,  the  program  calculates  the  body  wave  drag 
coefficient  for  the  particular  Mach  nusiber  specified. 

IV.  Test  Case  Results 

Several  body  configurations  were  Investigated, but  three  offer 
capabilities  and  will  alone  be  discussed  here.  The  first  body  Is 
revolution  specified  lr.  Figure  3. 

For  this  smooth  body,  the  program  generated  the  vs.  M curve 

w 


Figure  4.  Smooth  Body  Wave  Drag 

Values  near  the  speed  of  sound  are  invalid  due  to  the  limitations  In  the  supersonic 
area  rula  theory.  Also,  the  wave  drag  coefficient  decreases  with  Increasing  Mach  number 
within  the  velocity  range  Investigated. 

The  second  body  used  the  first  body  as  a baseline  but  had  an  axially  symmetric  pro- 
tuberance added  to  a forward  section  of  the  body.  The  additional  volume  of  the  protuber- 
ance, given  In  Figure  S,  was  a amall  football-shaped  volume  addition  distributed  equally 
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«bout  the  slender  baseline  body  axis. 

For  this  body  with  forward  protuberance, the  program  generated  the  CQ  vs.  M curve 

w 

depicted  In  Figure  6. 


Figure  6.  Body  With  Forward  Protuberance  Wave  Drag 


As  with  the  smooth -body  wave  drag,  the  same  trends  were  evidenced  for  the  body  with  the 
forward  protuberance. 

The  third  body  also  used  the  same  smooth  baseline  body  and  tha  Identical  football- 
shaped  volume  addition  but  located  It  In  an  aft  region  of  the  body.  Figure  7 depicts 
this  third  body.  ^ 

For  this  body  with  aft  protuberance, the  program  generated  the  vs.  M curve  depicted 

w 

In  Figure  8. 
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Figure  7.  Body  With  Aft  Protuberance  (not  to  scale) 
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Figure  8.  Body  With  Aft  Protuberance  Wave  Dreg 

The  shape  of  this  curve  was  radically  different  from  the  previous  two  cases.  The  general 
trends  still  remain;  however,  the  slope  of  this  curve  Is  noticeably  different. 

The  three  curves  are  shown  together  In  Figure  9 along  with  smooth-body  curves  from 
reference  5. 
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Figure  9.  Wave  Drag  Comparison 
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This  figure  clearly  shows  that  the  curve  trends  Identified  for  the  three  test  cases  are 
consistent  with  both  NASA  experimental  data  and  NASA  area  rule  theory.  If  the  NASA 
curves  are  extrapolated  toward  higher  Mach  numbers,  the  difference  between  results  becomes 
less  pronounced.  The  forward  protuberance  does  not  appear  to  Increase  the  wave  drag 
appreciably  over  the  velocity  range  studied.  Conversely,  the  aft  protuberance  does, 
especially  In  the  region  between  Mach  1.2  and  Mach  2,20.  The  qualitative  Judgment  can 
be  made  that  the  forward  location  Is  desirable  over  the  aft  location  If  the  specified 
protuberance  Is  to  be  added.  Thus,  this  program  offers  an  additional  tool  to  the 
designer  with  respect  to  cosiponent  location  for  minimum  wave  drag, 

V.  Conclusions 

The  sensitivity  of  this  program  allows  even  the  slightest  protuberance  movement  to 
be  reflected  In  curve  movement  and  shape  changes.  Ideal  protuberance  locations,  over  a 
range  of  Mach  numbers,  could  possibly  be  theoretically  determined  for  least  wave  drag. 

So  far,  values  at  higher  Mach  numbers  more  nearly  approximate  experimental  data. 

VI.  Recommendations 
— 

Further  research  Is  needed  In  several  areas.  The  models  Investigated  In  this  report 

should  be  wind  -tunnel  tested  for  comparative  purposes  with  the  values  obtained  here. 

Additional  baseline  bodies  and  protuberance  shapes  need  to  be  Investigated  to  verify  the 

capabilities  of  this  program.  The  accuracy  of  both  the  area  determination  and  second 

derivative  calculation  phases  should  be  Increased  to  determine  the  worth  of  such  accuracy. 

A graphics  capability  to  generate  C vs.  M curves  for  calculated  values  would  drastically 

w 

Improve  the  program.  Finally,  a searching  procedure  to  locate  Ideal  protuberance  location 
In  order  to  achieve  minimum  wave  drag  over  a desired  Mach  number  range,  should  be  added 
to  the  program  to  make  It  virtually  self-contained. 
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APPENDIX 

COMPUTER  LISTING 
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A METHOD  FOR  MEASURING  ATHLETIC  POTENTIAL 
John  P.  Jackson*  and  James  E.  Scott** 

Abstract 

This  paper  explores  the  use  of  the  Keller  runner  model  to  obtain  physiological  para- 
meters for  use  as  an  athletic  aptitude  measure.  The  Keller  model  Is  presented  In  detail 
and  related  to  sprints  and  distance  runs.  In  this  way  four  physiological  parameters 
representing  strength,  muscular  endurance,  flexibility  and  cardiorespiratory  endurance 
can  be  assigned  numerical  values.  The  technique  of  assigning  values  Is  demonstrated  for 
fourteen  runners  on  which  the  appropriate  data  was  taken.  The  consequences  of  this  work 
towark  athletic  aptitude  measurement  is  discussed. 

I . Introduction 

There  are  a variety  of  techniques  to  measure  athletic  aptitudes  of  persons  enrolled 
in  physical  education  programs.  Often  these  are  expressed  in  terms  of  accomplishment  of 
some  feat:  running  so  fast  or  jumping  so  far.  Sometimes,  however,  athletic  aptitudes 
are  related  to  physiological  parameters  of  the  human  body  but  often  involve  expensive 
test  equipment  such  as  the  treadmill.  The  first  category  of  measurements  are  usually 
easier  to  make  while  the  second  correlates  with  physiological  processes  of  the  human 
body.  In  this  paper,  we  present  a technique  which  combines  the  best  of  each  category: 
convenience  of  measurement  and  tie-in  with  human  physiology.  Only  stop  watches  are 
required,  and  yet, the  final  output  of  the  measurement  process  is  four  physiological  para- 
meters characteristic  of  the  human  body. 

These  parameters  are  closely  associated  with  four  of  the  generally  accepted  com- 
ponents of  physical  fitness  (Ref  1).  They  are  Strength — the  capacity  of  a muscle 
to  exert  force  against  a resistance,  Muscular  Endurance--the  capacity  of  a muscle  to 
exert  a force  repeatedly  over  a period  of  time,  Flexibility — the  ability  to  use  a 
muscle  throughout  its  maximum  range  of  motion,  and  Cardiorespiratory  Endurance  — 
the  capacity  of  the  heart,  blood  vessels,  and  lungs  to  deliver  nutrients  and  oxygen  to 
the  tissues  and  to  remove  wastes. 


♦Captain,  USAF,  Associate  Professor  of  Physics,  DFP 
♦♦Captain,  USAF,  Director,  Human  Performance  Lab,  AHPR 
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II.  Theory 

A.  Basic  Equations ' 

The  athletic  measurement  process  discussed  herein  depends  upon  the  requirement  that 
a runner  must  conserve  energy  and  momentum.  We  use  the  analysis  expressing  this  conser- 
vation as  formulated  by  Joseph  B.  Keller  (Ref  2 and  3). 

A runner  moving  with  a speed  v(t)  will  cover  a distance  x in  time  t according  to  the 


relation 


r " 


The  speed  v(t),  however,  is  not  arbitrary,  but  must  satisfy  Newton's  Law,  that  is,  con- 
servation of  momentum: 


y 

■where  f (t)  is  the  propulsive  force  per  unit  mass  (which  the  runner  exerts)  and  — the 
friction  force  per  unit  mass.  Strictly  speaking,  f(t)  is  the  forward  contact  force  of 
the  running  surface  upon  the  runner;  but  by  Newton's  Third  Law,  f(t)  is  equal  in  magni- 
tude to  the  actual  force  the  runner  exerts  against  the  running  surface  (for  which  we 
assume  no  slip).  Like  Keller,  we  assume  that  body  friction  force  per  unit  mass  is  pro- 
portional to  speed  with  a constant  of  proportionality  1/x.  Thus,  the  friction  coeffi- 
cient r is  the  first  physiological  parameter  characteristic  of  the  runner  that  will  be 
measured. 

Following  Keller,  we  next  make  several  conditions  upon  v(t)  and  f(t).  First,  we 
shall  require  that  the  runner  starts  from  rest;  hence, 

v(0)  - 0 (3) 

«••  <U.  though  f(t)  is  under  the  control  of  the  runner,  we  assume  there  exists  an  upper 
m»<>  >a  ft  determined  by  the  runner's  physiology.  ThuB,  for  some  F, 

ft)  « F for  all  t.  (4) 

« IM  awed  stars  to  logical  parameter. 

•a.  e>»  Sill  rseeel  that  equations  (1)  through  (4)  are  insufficient  to 
««..  • . ' lining  situations.  During  a race,  the  runner  does  work 

w m#  h statrti  serose  to  deplete  hla  body's  energy  reserves.  In 
.....  | - i.  ability  to  acquire  energy  by  aerobic  means  if 
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f » F.  For  this  reason,  a runner  can  sprint  only  so  far  at  full  force  before  he  exhausts 
all  available  energy  resources.  Thus,  additional  equations  must  be  written  to  account 
for  the  energy  balance. 

If  E(t)  is  the  energy  per  unit  mass  available  for  running  at  some  time  t,  then, 
according  to  Keller's  simple  formula, 

S - 0 - fv  <5) 

where  o is  the  rate  energy  (per  unit  mass)  can  be  made  available  to  the  runner  by  the 

aerobic  process.  Obviously,  the  runner  always  has  an  initial  reserve  of  available  energy 

E , for  otherwise  he  would  be  unable  to  start  running  at  all.  Hence, 
o 

E(0)  - E , (6) 

o 

o and  Eq  are  the  final  two  physiological  parameters. 

In  order  to  couple  Equation  (5)  with  the  running  dynamics  expressed  by  Equations 
(1)  to  (4),  one  more  requirement  must  be  stated:  the  available  energy  must  always  be 
positive  during  a run.  If  E becomes  zero  at  any  time,  then  the  race  must  terminate.  That 
is, 

E(t)  >_  0 for  t <_  T (7) 

where  T is  the  duration  of  the  race. 

Equations  (1)  through  (7)  are  the  basis  for  interpreting  the  scientifically  con- 
trolled race  described  below.  We  obtained  volunteer  running  subjects  who  we  asked  to 
run  a set  of  four  distances,  D,  in  a minimal  time,  T.  For  short  distances,  this  implied 
that  f(t)  equaled  F.  For  long  distances,  f (t)  was  of  necessity  less  than  F;  but  in 
order  to  achieve  a run  in  minimal  time,  all  energy  should  have  been  used  up  at  time,T, 
the  time  for  the  entire  race.  That  is,  E(T)  • 0.  As  described  below,  these  two  condi- 
tions allowed  a unique  solution  of  Equations  (1)  through  (7),  and  from  these  solutions 
the  four  athletic  parameters  were  calculated  from  measurement  data. 

B.  Sprints 

First,  consider  the  short  distance  runs  (sprints).  With  f(t)  « F,  Equation  (2) 
can  be  combined  with  Equation  (3)  to  give 

v(t)  "Ft  (1  - e C^T).  (8) 

Placing  Equation  (8)  into  Equation  (1)  gives 

x(t)  - Fr2(t/T  + e~t/T-  1).  (9) 
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For  times  t «t,  the  Initial  acceleration  of  the  runner  is  governed  by  the  exponential 
term;  but  for  times  t >>x,  the  maximum  force  is  balanced  by  friction  at  a constant  run- 
ning speed  and,  as  can  be  seen  from  Equation  (9),  x is  linear  with  t. 

The  energy  Equation  (5)  can  be  combined  with  Equation  (2)  and  the  result  integrated 
to  give 


E(t) 


2 i f 2 

+ at  - t v - Mv  ds> 


(10) 


For  the  sprint,  v(t)  is  given  by  Equation  (8)  so  that  Equation  (10)  can  be  integrated  to 

give  E(t).  However,  as  Already  indicated,  there  will  be  a maximum  time.T  , for  which  a 

c 

runner  can  spring  with  f(t)  « F.  Using  the  relation 
E(Tc)  = 0, 

we  calculate  from  Equation  (10) 


(ID 


F r - a 


(12) 


(where  it  is  assumed  that  Tc  >>x).  Then,  if  we  insert  Tc  into  Equation  (9),  we  compute 
the  maximum  distance,  Dc  **  x (Tc),  for  which  a runner  can  sprint  with  f(t)  = F: 


, . (13) 

1 - cr/(F  t) 

Though  this  distance  obviously  changes  from  runner  to  runner,  we  interpret  it  as  the 
natural  dividing  line  between  sprint  and  distance  runs. 

C.  Distance  Runs 

Now,  consider  the  distance  run  for  which  f (t)  < F.  Using  optimal  control  theory, 
Keller  showed  that,  with  the  exception  of  the  start  and  finish  (which  may  be  ignored 
here) , the  speed  v of  the  runner  must  be  constant  in  order  to  minimize  the  time  T for  the 
given  distance  D > Dc>  Perhaps  we  can  make  this  plausible  by  the  following  argument: 
suppose  the  runner  runs  "slow"  for  part  of  the  race  to  conserve  energy.  It  is  obvious 
he  cannot  run  "slow"  for  the  entire  race,  for  then  he  would  be  unable  to  totally  deplete 
his  available  running  energy  by  the  end  of  the  race.  This  means  he  would  then  have 
energy  left  over  which  could  have  been  used  to  give  him  a better  time.  Thus,  he  must 
speed  up  during  some  other  interval  of  the  race.  However,  according  to  Equation  (10), 
the  runner  consumes  his  available  energy, E(t),  as  the  speed  squared;  whereas,  according  to 
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Equation  (1),  he  only  buys  an  advantage  in  a shorter  time  according  to  the  first  power 
of  his  velocity.  This  strategy  of  variable  speed  wastes  energy  during  times  when  he 
must  run  fast.  Hence,  he  would  make  optimal  use  of  his  energy  (and  guarantee-  the  shor- 
test time)  by  running  at  a constant  speed. 

Substituting  v - constant  into  Equation  (10),  and  using  E(T)  - 0,  we  get 

i 2 v^ 

Eq  + aT  - jv ^ — T = 0 (14) 


from  which  we  calculate  v to  be  (using  T/t  >>  1). 


v 


+ aT)  ' 
T 


Since  v - constant  « D/T, 


D 


I 1/2 

(E  + aT)  T ' 
o 


(15) 


(16) 


in  agreement  with  Keller. 

III.  Measurement  of  the  Physiological  Parameters 

Equations  (9)  for  sprints  and  (16)  for  distance  runs  show  how  running  distance  is 
related  to  time.  Since  these  relations  involve  the  four  physiological  parameters 
F,  t,  a,  and  Eq,  it  should  be  possible  to  compute  these  parameters  for  any  given  runner. 
A.  Sprint 8 

In  the  case  of  sprints,  the  exponential  term  of.  Equation  (9)  can  be  neglected  for 
times  t >>t.  Under  this  approximation,  Equation  (9)  becomes 

t - ( j^-)  x + t (17) 

Thus,  if  x(t)  were  plotted  for  a given  runner,  the  slope  of  the  line  would  be  (Ft)-1  and 
the  t-intercept  would  be  x.  Following  this  reasoning,  we  timed  volunteer  runners  using 
stopwatches  at  intervals  of  10  meters  as  they  sprinted  a 100  meter  dash.  We  analyzed  the 
resulting  data  according  to  the  method  of  least  squares  and  calculated  the  best  linear 
fit  to  the  data.  Since  the  method  of  least  squares  gives  the  error  in  both  the  slope 
and  t-intercept,  we  were  able  to  compute  the  statistical  errors  in  F and  x.  Figure  1 
shows  a typical  graph  of  the  least  squares  fit  to  a set  of  collected  data  on  one  individ- 
ual. 
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x (meters) 


Figure  1.  Typical  time  versus  distance  curve  for  volunteer  subject.  Line  represents 

Least  Squares  fit  to  data  measured  by  timers  stationed  at  ten  meter  Intervals. 


B.  Distance  Runs 

For  distance  runs,  we  had  to  collect  and  analyse  the  data  quite  differently.  It  was 
not  possible  to  compute  o and  Eq  from  just  one  race  as  for  F and  t.  The  reason  was  that 
different  running  speeds  must  be  used  for  different  distances  and  hence  we  could  not 
derive  one  equation  analogous  to  Equation  (17). 

Using  the  same  volunteers  as  in  the  sprints,  we  measured  the  time  for  run  distances 
of  400,  800  meters,  and  the  mile  and  a half  (on  separate  days,  of  course).  Actually, 
only  two  races  are  required  to  calculate  the  two  parameters  o and  Eq,  but  we  used  a 
third  to  estimate  the  error  involved.  We  insured  that  the  condition  E(T)  - 0 (Equation 
(14)  ) was  well  approximated  by  asking  all  subjects  how  much  faster  they  thought  they 
could  have  run  the  distance  race;  in  all  cases  the  estimated  time  decrease  was  much 
smaller  than  the  total  running  time  and  so  we  assumed  that  Equation  (14)  was  adequately 
satisfied. 
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Figure  2.  Distance  versus  time  graphs  of  Equation  16  drawn  through  three  measured  data 
points.  Corresponding  to  each  curve  (only  two  of  three  possible  are  shown) 
is  a unique  set  of  O and  E0.  Note  that  each  curve  passes  through  only  two  of 
the  three  points  (since  each  curve  is  uniquely  determined  by  two  parameters 
<7  and  E0.) 

Through  any  two  points,  it  is  possible  to  plot  a curve  according  to  Equation  (16) 
which  corresponds  to  a unique  set  of  a and  Eq.  Thus,  as  shown  in  Figure  2,  we  con- 
structed three  unique  curves  from  the  three  data  points,  each  corresponding  to  a 
specific  o and  Eq.  Since  we  could  not  establish  which  curve  was  the  best  characteriza- 
tion of  the  given  runner,  we  treated  each  curve  as  having  an  equal  a priori  probability 
of  being  correct  and  simply  computed  the  average  a and  Eq  as  well  as  their  standard 
deviations  as  an  estimate  of  the  error.  It  should  be  pointed  out  that  we  also  took  into 
consideration  the  error  involved  in  our  previous  measurements  of  t.  This,  of  course, 
was  necessary  since  Equation  (16),  from  which  a and  Eq  were  calculated,  also  involves  t. 
Hence,  our  reported  error  was  increased  accordingly. 

IV.  Results 

Figure  3 shows  the  results  of  our  measurements  for  F,  x,  a,  and  Eq.  Though  small, 
the  data  base  of  fourteen  individuals  (eight  male,  six  female)  reveals  some  noteworthy 
results.  First,  the  distribution  of  mean  values  for  all  parameters  fluctuate  signifi- 
cantly more  than  the  associated  error  limits.  This  means  that  this  method  of  measuring 
athletic  performance  based  on  simple  physiological  parameters  is  capable  of  discrimi- 
nating between  individuals.  Second,  the  data  shows  that  Equations  (9)  and  (16)  predict 

the  general  features  of  how  a runner  covers  a specific  distance  with  time.  If  this  were 
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WOMEN  (6) 


WOMEN  (6) 


W0?1EN  (?) 


WOMEN  (6) 


Measured  values  of  physiological  parameters  for  each  subject. 
Shaded  area  Indicates  uncertainty  In  each  parameter;  mean 
value  is  at  center.  Superman  values  arer-  .89  s,  F ■ 12.2  m/s 
E ■ 575  Cal/Kg,  <r  - 9.9  Cal/Kg  s (Ref.  3).  Point  values  for 
athletic  evaluation  scheme  are  shown  at  right  of  each  chart. 


not  so,  then  our  first  observation  would  not  be  possible  to  make.  It  would  seem,  there- 
fore, that  other  physiological  factors  (which  certainly  must  play  a role  in  a runner's 
performance  but  were  not  treated  within  the  mathematical  model)  probably  are  second  - 
order  effects.  Third,  according  to  Figure  4,  there  appears  to  be  little  or  no  correla- 
tion between  F and  t,  and  a and  E , that  is,  between  the  parameters  associated  with 


anaerobic  and  aerobic  running. 

Relating  the  four  parameters  to  the  components  of  physical  fitness  mentioned  pre- 
viously, we  can  make  the  following  associations:  the  parameter  F is  most  closely  related 
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to  the  strength  component;  the  parameter  t Is  probably  greatly  affected  by  the  flexibil- 
ity component;  the  parameter  a is  the  aerobic,  or  cardiorespiratory  endurance;  and  Eq 
is  anaerobic,  or  muscular  endurance.  However,  it  appears  as  though  there  may  be  some 
correlation  between  F and  o.  That  is,  it  appears  as  though  a strong  runner  tends  to 
have  a high  aerobic  conversion  efficiency.  This  is  reasonable  since  runners  who  develop 
a by  long  distance  running  also  probably  develop  F as  well.  But  this  does  not  mean, 
however,  that  a good  sprinter  is  a good  distance  runner  because  the  effect  of  body 
friction  must  be  taken  into  account.  Fourth,  on  the  basis  of  our  limited  data  set,  it 
appears  as  though  men  and  women  differ  significantly  only  in  two  of  the  four  parameters. 
As  indicated  in  Table  1,  these  parameters  are  a and  t;  women  tend  to  have  a smaller  t 
(greater  frictional  resistance)  and  a lower  o (reduced  aerobic  conversion  efficiency). 
Women's  strength  per  unit  mass  F and  stored  energy  per  unit  mass,  E^,  however,  seem  com- 
parable to  men's.  Table  1 also  indicates  that  the  calculated  conversion  distance,  Dc , 
between  sprints  and  distance  runs  is  about  the  same  for  men  (233  + 44  meters)  and  women 
(220  + 33  meters). 


Table  1 

Comparison  of  Men  and  Women's 
Physiological  Parameters 


Young 
Women : 


Young 

Men: 


.956+  .066s 


1.057+  .081s 


7.26+  .87m/s 


7.76+  .46m/s 


286.2+  55.8 


cal 

kg 


309.9+  40.9 


cal 

kg 


2.95+  .68 


cal 

kgs 


4.18+  .77 


cal 

kgs 


220+  33m 


233+  44m 


Fifth,  the  ratio  of  aerodynamic  friction  per  unit  mass,  — A (p,  density  of  air; 
v,  runner's  sprint  speed;  A, cross-sectional  area  of  sprinter;  M,mass  of  runner),  to  the 
total  body  friction  v/t  Is  no  greater  than  about  10Z;  thus,  most  of  the  frictional 
resistance  of  a runner  must  be  internally  generated.  Sixth,  Figure  3 also  shows  the 
values  calculated  by  Keller  for  a hypothetical  athlete  who  runs  every  race  in  world 
record  time  regardless  of  distance.  Appropriately  enough,  we  refer  to  him  as  Superman; 
his  parameters  represent  the  "ideal"  of  human  physiology. 

V.  Recommendations 

Like  measurements  of  academic  potential,  this  method  of  determining  athletic 
potential  does  not  produce  highly  discriminating  results.  Nevertheless,  it  does  seem 
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feasible  to  construct  an  evaluation  system  which  arranges  athletic  potential  in  broad 
categories.  To  illustrate  how  this  might  be  done,  we  have  constructed  the  following 
simple  evaluation  scheme.  We  are  not  suggesting  that  this  is  the  most  suitable  method, 
but  it  does  show  how  athletic  aptitude  might  be  displayed.  In  our  scheme,  we  assign  the 
measurement  of  each  physiological  parameter  an  Integer  between  one  and  ten.  We  chose 
the  Interval  between  each  integer  to  correspond  to  the  approximate  measurement  error  of 
the  parameter,  for  it  is  of  no  use  to  be  concerned  about  fractions  of  a point  since,  as 
we  have  shown,  the  measurements  are  not  that  accurate.  We  arbitrarily  selected  the 
number  five  to  be  the  approximate  mean  value  of  all  subjects,  regardless  of  sex.  Figure 
3 shows  the  way  we  assigned  the  integers  to  the  four  physiological  parameters.  Then,  we 
summed  the  integer  values  for  each  parameter  and  plotted  the  resulting  distribution  as 
shown  in  Figure  5.  Assuming  that  better  athletic  potential  corresponds  to  higher  values 
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Figure  5.  Distribution  of  total  physiological  score 

according  to  sex.  Two  overlaplng  distributions 
are  suggested,  one  for  women  (14-20)  and  the 
other  for  men  (19-24). 

of  each  parameter,  then,  it  would  seem  that  the  higher  the  summed  score,  the  better  the 
overall  athletic  aptitude.  We  note  further  that  the  use  of  a summed  score  allows  defi- 
ciencies in  one  physiological  parameter  to  be  compensated  by  a stronger  value  in  another. 

Figure  5 shows  the  kind  of  distribution  of  parameters  which  can  be  expected  using 
this  running  model.  Differences  of  one  point  value  on  the  summed  scores  correspond  to 
significant  differences  in  at  least  one  parameter  because  of  the  way  we  assigned  the 
integers  to  the  measured  data.  It  is  noteworthy  that  two  distributions  seem  to  be 
j apparent:  one  centered  at  15  (an  average  of  3.75/parameter)  and  the  other  at  20  (an 
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average  of  5/parameter).  The  lower  distribution  consists  entirely  of  women;  the  upper 
consists  of  eight  men  and  two  women.  Thid  indicates  that  men  and  women  should  be  evaluated 
separately,  a fact  not  readily  obvious  since  the  physiological  parameters  do  not  depend 
on  body  mass. 

An  obvious  question  is  what  does  one  do  with  these  data.  We  can  only  provide  a 
partial  answer  at  this  time.  Suppose  a subject  has  a point  rating  of  18  which  is  dis- 
tributed as  follows:  t - 3,  F = 5,  E =>  3,  o “ 7 (Recall  from  Figure  5,  3.75  is  average 

o 

for  women,  5 for  men).  If  the  subject  is  a woman,  her  athletic  aptitude  should  be  con- 
sidered high  and,  therefore,  it  would  probably  be  difficult  to  raise  her  value  appre- 
ciably (note  that  she  is  already  almost  average  or  above  in  all  parameters.)  But  if  the 
subject  is  a male,  then  his  summed  score  might  be  raised  considerably  if  certain  para- 
meters could  be  Increased.  These  probably  would  be  x and  Eq  since  F and  a are  already 
average  or  above.  Thus,  an  athletic  program  for  him  should  involve  those  things  which 

tend  to  increase  his  x and  E . More  research  is  needed  to  determine  what  these  things 

o 

might  be,  but  it  probably  isn't  too  difficult  to  make  an  intelligent  guess  knowing  the 
role  each  parameter  plays  In  the  dynamic  running  equations  (see  Table  2). 


Table  2 

The  Role  and  Suggested  Method  of 
Increasing  each  of  the  Physiological  Parameters 


Parameter 

Role 

Suggested  Method  of 
Increasing 
Parameter 

T 

Internal  body  friction, 
flexibility 

Gymnastic  exercise, 
limbering  up. 

F 

Body  strength  per  unit 
mas  8 

Weight  training. 

E 

o 

Energy  stored  in  body 
available  for  physical 
exercise,  muscular 
endurance . 

Diet,  rest. 

0 

Aerobic  efficiency, 

cardiorespiratory 

endurance. 

Distance  running, 
distance  swimming. 
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Finally,  it  should  be  pointed  out  that  it  is  possible  to  determine  via  Equations 
(9)  and  (16)  the  effect  on  running  time  of  a change  in  the  physiological  parameters. 

Thus,  an  indication  how  a given  exercise  program  will  translate  into  athletic  performance 
is  possible.  This  could  be  of  value  to  track  athletes  in  particular. 

In  conclusion,  we  present  this  study  of  a way  to  measure  athletic  potential  as  a 
feasibility  study,  not  necessarily  the  way  in  which  such  a method  would  actually  be 
implemented.  Questions  still  remain.  Could  more  accurate  measurements  be  made, thereby 
reducing  the  error  of  each  parameter?  How  would  the  data  distribute  for  a much  larger 
sample?  Do  the  suggested  exercises  of  Table  2 really  affect  the  parameter  significantly? 
Are  the  measured  values  of  a physiological  parameter  repeatable  for  a given  individual? 
Can  the  progress  of  training  be  followed  by  measuring  athletic  potential?  How  applicable 
is  the  friction  coefficient  t (as  measured  for  sprints)  to  distance  runs?  In  spite  of 
these  and  other  unanswered  questions,  it  does  appear  to  us  that  the  method  described 
herein  could  be  used  as  an  inexpensive  method  of  measuring  athletic  aptitude  and  obtain- 
ing useful  physiological  Information. 
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DESIGN  PROCEDURE  FOR  COUPLING 
SAVONIUS  AND  DARRIEUS  WIND  TURBINES 

. . * ** 

R.  C.  Oliver  and  P.  R.  Nixon 


Abstract 

A design  procedure  for  a combined  Savonius  and  Darrieus  Wind  Turbine  System  is 
presented  and  demonstrated.  Individual  system  characteristics  are  described  and  dis- 
cussed. System  characteristics  are  combined  and  procedures  prescribed  to  design  either 
maximum  power  output  or  minimum  starting  time  systems.  Generalizations  are  provided 
for  system  design  when  specific  test  data  does  not  exist  and  a possible  application  of 
the  resulting  system  is  described. 

I.  Introduction 

The  recent  upsurge  of  interest  in  alternate  energy  sources  has  included  wind  power. 
Historical  uses  included,  among  others,  grinding  grain  and  pumping  water.  This  use  was 
discontinued  to  a large  degree  as  thermodynamics  advanced  and  engines  of  high  power 
capacities  were  developed.  Wind  power  could  not  compete  with  low-cost  fuels  and  the 
high-capacity  engines.  Recent  developments  increasing  wind  turbine  efficiency  and 
the  continued  rise  of  conventional  fuel  costs  have  made  wind-power  systems  a poten- 
tially viable  power  source  for  certain  situations.  There  is  little  doubt  as  to  the 
technical  attainability  of  wind  energy  today,  but  the  question  of  economic  compet- 
itiveness for  general  use  remains  (Ref.  1).  The  suitability  of  wind  power  to  supple- 
ment Air  Force  energy  requirements  has  not  been  generally  determined;  however,  several 
current  applications  exist  where  wind  power  may  be  a viable  alternative  or  supplement 
for  conventional  fuels.  Perhaps  the  most  obvious  use  is  a remote-site  power  source. 
Wind  energy  is  quiet,  unobtrusive  and  has  no  identifying  heat  or  noise  emissions. 
Depending  on  the  location,  the  wind  could  supply  all  or  a major  portion  of  the  power 
required.  Even  used  as  a supplemental  power  source,  a wind  system  can  result  in  a 
significant  reduction  of  site  logistical  support.  Supply  and  logistical  support  of 
many  such  sites  is  costly,  and  often  fuel  must  be  transported  over  large  distances  or 
difficult  terrain.  The  amount  of  power  needed  is  generally  not  large  and  within  the 
capabilities  of  moderate  to  small  wind  systems.  In  such  a situation  a wind  system 
may  be  the  optimal  choice.  Used  alone,  some  type  of  energy  storage  system  sized  for 
windless  periods  would  be  required. 
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II.  Background 

With  the  remote  site  application  in  mind,  a variety  of  wind-power  systems  were  con- 
sidered. The  Darrieus  Vertical  Axis  Wind  Turbine  (DVAWT)  seems  to  be  the  most  promising. 
Although  invented  in  1926,  new  interest  has  triggered  significant  study  and  advancement 
since  1971.  Large  programs  by  Sandia  Laboratories  and  the  Canadian  National  Research 
Council  have  been  supplemented  by  NASA  Langley,  Oregon  State  University  and  other 
organizations  (Ref.  2).  The  DVAWT,  recognized  by  its  eggbeater  shape  (see  Figure  1), 
has  several  important  advantages  over  other  systems.  First,  due  to  shape,  the  wind  can 
be  accepted  from  any  direction, and  no  yaw  control  system  is  required  to  orient  the 
turbine  into  the  wind.  Additionally,  the  relatively  simple  tower  construction,  delivery 
of  the  mechanical  energy  naturally  to  the  ground,  and  the  lack  of  a pitch  control 
system  for  synchronous  applications  contribute  to  the  advantages  of  the  DVAWT  (Ref.  3). 

As  a result  of  these  factors, the  power  output  per  rotor  weight  and  cost  are  among  the 
best  (Ref.  4),  resulting  in  a cost  which  may  approach  one  sixth  that  of  horizontal 
axis  systems  (Ref.  5,  6 and  7).  One  item  not  yet  addressed,  however,  is  the  inability 
of  a Darrieus  system  to  self  start,  even  in  gale  force  winds.  The  task  reported  here 
began  with  selection  of  a starter  system  for  a remote-site  Darrieus  wind -power  system. 

III.  Theory 

A variety  of  concepts  from  electric  motors  to  aerodynamic  systems  was  considered 
as  possible  starting  systems.  The  specific  application  greatly  influenced  the  choice 
as  did  the  previous  selection  of  the  primary  system.  An  electric  starter  with  its 
control  and  energy  supply  systems  would  add  complexity  and  Increase  system  weight  while 
decreasing  portability.  An  aerodynamic  starter  would  have  to  be  compatable  with  the 
vertical  axis  primary  system  (a  horizontal  turbine  would  require  a yaw-control  system 
and  a transfer  system  to  start  the  Darrieus.)  A self-starting  vertical  axis  turbine, 
however,  avoids  these  complications  and, therefore, a Savonius  vertical  axis  turbine  was 
selected.  Invented  in  1929, and  finding  a multitude  of  uses  since  then  (Ref.  8), the 
Savonius  Vertical  Axis  Wind  Turbine  (SVAWT)  is  capable  of  self-starting  although  it 
operates  at  a drastically  different  tip  speed  range  than  the  DVAWT.  The  tip  speed 
parameter  is  a ratio  of  the  device  rotational  velocity  to  the  wind  speed.  If  used, 
the  combined  system  must  permit  each  device  to  operate  in  its  own  regime,  or  drag 
rather  than  thrust  will  result. 

Luckily,  the  performance  characteristics  of  both  the  Darrieus  (Ref.  9,  10)  and  the 
Savonius  (Ref.  11,  12)  are  known  and  documented  for  positive  power  conditions.  Figure  2 
illustrates  the  performance  of  characteristic  Darrieus  and  Savonius  turbines.  The 
optimum  tip-speed  ratio  for  the  best  performance  (maximum  power)  occurs  at  about  0.85 
for  the  Savonius  and  about  6.0  for  the  Darrieus. 
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Figure  1.  Schematic  of  Combined  Darrieus  and  Savonlua  Self -Starting  Wind  Turbine 
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Maximum  power  Cor  Che  combined  system  would  be  produced  when  both  rotors  operate  at 
or  near  their  respective  tip-speed  ratios  for  maximum  power.  The  minimum  starting 
time  condition  for  the  combined  system,  however,  would  exist  when  maximum  Savonlus 
power  is  produced  during  the  Darrleus  starting  period.  Both  of  these  cases  are  addressed 
below.  Since  both  the  DVAWT  and  the  SVAWT  are  on  the  same  axis,  the  Individual  radii 
must  be  adjusted  to  provide  different  tip-speed  ratios  for  the  common  angular  velocity. 
Consider  the  maximum  power  case,  the  tip-speed  ratio,  TSR,  Is 


TSR  - 


R U 
V 


(1) 


where  R Is  the  device  (either  Savonlus  or  Darrleus)  Radius,  a)  the  angular  velocity 
(common  for  both  systems)  and  V the  wind  velocity  (also  common).  In  order  for  both  the 
SVAWT  and  the  DVAWT  to  operate  at  their  maximum  power  conditions,  the  ratio  of  the  SVAWT 
radius  to  the  DVAWT  radius  must  be  the  same  as  the  ratio  of  the  tip-speed  ratios  corre- 
sponding to  their  respective  maxlum  conditions.  In  equation  form 

SVAWT  m SVAWT.  Max  Pwr  ... 

d 

DVAWT  DVAWT,  Max  Pwr 

For  the  devices  cited, 


R 

r" 


SVAWT 

DVAWT 


0.85 

6.0 


0.142 


(3) 


Therefore,  for  a Darrleus  rotor  radius  of  10  meters  the  Savonlus  rotor  radius  for  the 
maximum  power  condition  should  be  1.42  meters.  This  radius  combination  permits  both 
devices  to  operate  at  maximum  power  output  conditions  simultaneously. 


The  minimum  starting  time  condition  exists  when  maximum  excess  torque  occurs  during 
the  period  when  the  Darrleus  is  accelerating  to  starting  or  self-sustaining  speed.  The 
condition  where  this  occurs  will  be  analysed  for  a specific  case  Illustrating  the 
general  procedure  for  known  device  performance  characteristics.  Additionally,  a simple 
design  technique  Is  presented  which  provides  a reasonable  minimum  starting  time  design. 

For  the  combined  system  to  start  and  accelerate  ,the  torque  produced  by  the  Sevonlus 
end  Darrleus  must  exceed  that  required  to  overcome  friction.  The  Darrleus  produces  no 
positive  torque  at  low  tip -speed  retlos.  Although  the  Darrleus  torque  may  even  be 
slightly  negative,  by  far  the  most  Important  factor  during  starting  Is  the  system 
friction  (Ref.  13).  Knowing  the  performance  characteristics  of  the  Darrleus  and  the 
Savonlus  and  the  system  friction,  we  can  estimate  the  time  required  to  eccelerate  the 
system  to  self-sustaining  speed,  where  the  Darrleus  produces  sufficient  torque  to 
overcome  friction  and  continue  accelerating  to  design  conditions.  The  difference 
between  the  combined  or  total  torque  produced the  friction  torque,  TprlctloB* 


t 
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Is  a measure  of  the  systems  ability  to  accelerate,  this  can  be  stated  mathematically  as 


Total 


'Friction 


tftere  L is  the  system  polar  moment  of  Inertia  and  a the  angular  acceleration 


where  u is  the  angular  velocity ,and  t is  time 
the  acceleration. 


The  net  torque,  T^^,  then,  provides 


Total  Friction 


A constant  acceleration  or  a short  time  period  is  assumed,  to  permit  the  difference 
form  for  the  angular  acceleration. 


For  a given  system  the  performance  characteristics  can  be  determined  based  on  the 
selected  radius  ratio  and  the  starting  time  estimated  using  equation  5.  In  order  to 
better  understand  the  process, a sample  problem  for  a representative  DVAWT  system  is 
provided. 


IV.  Example  Calculation 


The  characteristics  of  the  system  are  presented  in  Table  1 


TABLE  I 


CHARACTERISTICS  OF  A DARRIEUS  WIND 
TURBINE  FOR  REMOTE-SITE  APPLICATION 


Height 


Radius: 

Polar  Moment  of  Inertia 


Starting  Hind  Velocity: 

Saif -Sustaining  Spaed 
Frictional  torque,  including  electrical 
generation  ays tarn: 

Operational  Altitude 
Tip  Spaed  Ratio: 


Botii  Darriaue  and  Savon 1 us  characteristics  are  known  and  presented  as  a plot  of  power 
coefficient,  C versus  the  tip-speed  ratio  and  expressed  as: 


idiere  T is  the  device  torque,  y the  angular  velocity,  V the  wind  speed,  p the  air 
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density,  R the  device  radius,  and  H the  height.  The  data  are  experimentally  obtained. 
For  any  given  system  we  can  calculate  the  combined  torque.  The  torque  required  to  over- 
come friction  can  be  estimated,  but  Is  more  accurate  if  measured,  when  possible.  Having 
the  net  torque  and  polar  moment  of  Inertia,  we  can  use  Equation  5 to  solve  for  the 
required  time. 

Table  2 is  an  example  calculation  for  the  Darrieus  data  given  in  Table  1 combined 
with  a Savonlus  having  a radius  of  1.79  feet  and  a height  of  6 feet.  Using  a Savonlus 
height  of  half  the  height  of  the  Darrieus  prevents  possible  interference  with  the 
highly  loaded  Darrieus  center  section.  The  Savonlus  starter  would  be  split  and  mounted 
at  the  top  and  bottom  portions  of  the  Darrieus,  (see  Figure  3).  The  common  angular 
velocity  is  co.and  X,  Cp,  T represent  the  tip-speed  ratio,  power  coefficient  and  torque 
of  the  Savonlus  and  Darrieus.  The  total  torque  Is  reduced  by  the  frictional  torque 
yielding  the  net  torque.  Equation  5 Is  then  used  to  solve  for  the  time  required.  The 
starting  time  Is  the  sum  of  all  times  before  the  starting  point.  The  Darrieus  Is 
considered  started  at  a Darrieus  tip-speed  ratio  of  about  4.7.  Note  (Table  2)  that  at 
this  point  the  Darrieus  torque  has  begun  to  exceed  friction  with  sufficient  excess  to 
provide  continued  acceleration. 


TABLE  2 


Sample  Calculations  for  Combined  System 


TRad/Sec] 

SAVONIUS 

DARRIEUS 

TOR 

QUE 

Time 

At 

Sec 

* 

X 

• 

* 

C 

p. 

T. 

** 

*D 

irk 

CPD 

Ft\f 

Total 

TX 

Ft  lbf 

Net 

Tn 

Ft  lbf 

1.64 

0.2 

0.085 

3.386 

0.67 

0 

3.386 

1.586 

6.78 

3.28 

0.4 

0.16 

3.1864 

1.34 

0 

3.1864 

1.386 

7.76 

4.92 

0.6 

0.22 

2.92 

2.0 

0 

2.92 

1.12 

9.61 

6.56 

0.8 

0.24 

2.39 

2.68 

0 

0 

2.39 

0.59 

L8.24 

8.2 

1.0 

0.23 

1.83 

3.35 

0.84 

1.2 

3.032 

1.232 

8.72 

9.83 

1.2 

0.21 

1.46 

4.02 

0.125 

3.78 

5.38 

3.58 

3.0 

11.4 

1.4 

0.16 

0.9 

4.69 

0.22 

6.67 

7.57 

5.77 

1.86 

*Data  Reference  12,  pg  86 
**Data  Reference  13,  pg  37 


i 


i 


\ 
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The  total  time  until  the  system  starts  for  this  case  Is  estimated  at  about  56  seconds. 
Any  other  case  could  be  Investigated  In  the  same  way  by  Inputting  the  new  radius  and 
height.  The  specific  performance  characteristics  of  the  devices  also  affect  the  results; 
thus.good  data  is  a necessity. 


In  the  general  case,  or  for  the  designer  who  lacks  comparative  data,  there  exists 
a simpler  technique  for  estimating  the  minimum  starting  time  and  maximum  power  conditions. 
Approximate  data  exist  with  respect  to  the  tip-speed  ratio  for  maximum  power.  The 
Savonius  ratio  for  maximum  power  varies  only  fractionally  with  size,  so  a ratio  of  0.85 
will  always  be  close.  For  a Darrleus  of  less  than  5 meters  radius  a ratio  of  about 
5.5  to  6.0  will  result  in  maximum  power.  In  general, a Darrleus  will  start  at  a ratio 
greater  than  three  for  a large  rotor  (Ref.  14),  nearing  4.5  as  the  rotor  size  decreases. 
Based  on  these  guidelines,  we  could  choose  the  following  conditions. 


Maximum  Power: 


SVAWT 
R DVAWT 


X 

X 


XVAWT. 

DVAWT, 


max  pwr  “ 
max  pwr 


0.85 

5.75 


0.148 


(7) 


So  for  a small  turbine  a Savonius  with  15%  the  radius  of  the  Darrleus  provides  maximum 
combined  power. 


Minimum 
Starting  Time: 


SVAWT  X SVAWT. 

” r ’positive  pwr 

DVAWT  DVAWT, g tart 


.311  (8) 


For  minimum  starting  time  a Savonius  radius  of  about  30%  of  the  Darrleus  rotor  should 
be  used.  XSVAWT  for  positive  power  Is  the  highest  value  of  tip-speed  ratio  for  which 
the  Savonius  produces  useable  positive  power  and  la  generally  a constant. 


Although  these  results  are  based  on  general  conditions , an  error  of  a few  percent 
does  not  have  a significant  Input  on  the  results.  Recall,  however,  that  this  case  Is 
based  on  a Savonius  height  of  half  of  the  Darrleus  height. 

Having  described  the  dynamic  starting  problem  for  a specific  case.  Including  sev- 
eral general  guidelines,  it  Is  now  necessary  to  assure  a static  start  for  the  design 
chosen. 

A ststic  starting  check  of  the  system  can  be  made  using  the  functional  relationship 
proposed  by  Blackwell  (Ref.  15)  similar  to  Equation  6. 

CO  - ^ (») 

09  % p V*  R Ag 

Where  CQ  Is  an  axperimently  determined  static  torque  coefficient,  Q Is  the  net  static 
torque  produced,  A(  Is  the  turbine  swept  area,  and  the  other  terms  have  been  previously 
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defined.  The  swept  area  is  twice  (2  buckets)  the  height  tines  the  radius 


Using  the  minimum  torque  coefficient  suggested,  0.25,  and  the  device  parameters  as 
specified, including  a height  of  6 feet,  we  find  that  the  static  torque  is  2 ft-lbf 
•which  is  sufficient  to  overcome  the  frictional  drag  of  1.8  ft-lbf.  In  «ll  cases 
examined, the  device  designed  for  dynamic  starting  provided  sufficient  static  torque 
for  starting.  An  actual  device  was  constructed  based  on  the  described  design  example 
Details  are  provided  in  a subsequent  section  and  Figure  3 is  a picture  of  the  device. 


One  additional  characteristic  of  the  combined  system  which  must  be  described  is  its 
ability  to  self-regulate  or  prevent  an  overspeed  condition  naturally,  without  controls. 
Both  the  Darrieus  and  Savonlus  contribute  to  this  ability.  Firstly,  the  Darrleus 
blades  are  high  lift-to-drag  ratio  sections  for  irtilch  the  lift  and  drag  forces  are 
strong  functions  of  the  local  blade  angle  of  attack.  The  velocity  of  the  blade, due 
to  rotation  and  the  wind  velocity, combine  to  form  the  resulting  total  flow  which  strikes 
the  blade.  At  low  rotational  speeds  a small  angle  of  attack  results  and  the  lift-to- 
drag  ratio  is  small  (as  in  the  previously  described  case  where  the  torque  is  slightly 
negative.)  As  the  rotational  velocity  approaches  the  design  point, the  angle  of  attack 
increases  and  produces  a large  lift-to-drag  ratio  providing  tha  positive  torque  regime. 
As  rotational  speed  increases  further,  the  angle  of  attack  Increases, causing  a stall 
of  the  blades,  reducing  the  lift-to-drag  ratio  and  preventing  an  overspeed  or  runaway 
condition.  Basic  airfoil  theory  explains  both  the  starting  problem  and  the  natural 
speed  limiting  (Ref.  16). 


The  Savonlus  also  contributes  to  this  Inherent  speed  limiting;  although  not  ex- 
plained aa  clearly  in  tans  of  tha  angle  of  attack,  tha  Savonlus  also  produces  negative 
torque  or  drag  once  its  tip-speed  ratio  exceeds  that  where  positive  power  is  produced, 
bee  Figure  1.  Although  the  tip-speed  ratio  at  which  this  occurs  can  be  estimated,  the 
exact  performance  beyond  this  point  is  in  question.  It  is  postulated  that,  unlike  the 
Darrleus  Which  simply  stops  producing  torque,  the  Savonlus  actually  produces  negative 
torque  or  drag  beyond  its  maximum  tip-speed  ratio  for  positive  power.  No  data  was 
found  to  support  this  contention  because  experimental  results  have  bean  concerned  with 
the  positive  power  output  regime.  Not  only  will  this  additional  drag  help  prevent 
overapeedlng,  but  it  should  also  ba  considered  in  the  maximum  power  output.  For  example, 
tha  starter  (Savonlus)  designed  for  mlnlmimi  starting  time  will  be  in  the  dreg  regime 
during  the  Darrleus  msximmi  power  condition  and,  once  quantified,  should  ba  considered. 
The  ability  to  self-regulate  is  a significant  advantage  relative  to  other  wind  power 
systems  which  often  smploy  active  control  systems  to  prevent  overapeedlng. 


USAFA-TR-79-1 


V.  Construction  Details 


The  Savonlus  starter  designed  by  the  above  procedure  for  minimum  starting  time  was 
constructed  and  will  be  tested  on  a Darrieus  turbine  of  the  design  described;  results 
will  be  reported  when  available.  Two  Savonlus  rotors  of  three- foot  height  were  con- 
structed. Considerable  attention  was  given  to  a design  permitting  easy  assembly  of 
the  entire  apparatus  'a  the  field  by  a minimum  of  personnel.  The  starter  buckets  of 
the  two  Savonlus  rotors  were  mounted  90  degrees  to  each  other  to  provide  a wind- 
catching  surface  in  all  directions. 

Many  experiments  (specifically  Reference  17)  have  shown  that  there  exists  an  optimum 
gap  width  ratio  for  maximum  Savonlus  performance.  This  ratio  is  the  gap  distance 
between  bucket  inner  radii  divided  by  the  rotor  diameter,  see  Figure  4.  The  gap  width 
used  in  our  design  is  the  net  distance,  excluding  the  center  support  pole  which  was 
not  considered  in  the  quoted  reference.  Using  the  net  distance,  then,  the  gap  width 
ratio  should  be  between  0.1  and  0.15  for  maximum  performance  (Ref.  18). 

A minimal  load  analysis  of  the  Savonlus  rotor  buckets  indicated  that  0.063  inch 
aluminum  sheet  buckets  with  0.10  inch  aluminum  end  plates  would  suffice  under  the 
given  conditions  (Ref.  19). 

VI.  Sumnarv 

The  design  procedure  described  provides  for  the  combination  of  Darrieus  and 
Savonlus  Vertical  Axis  Wind  turbines.  The  self-starting,  self-regulation  combined 
system  has  a variety  of  uses  although  this  specific  system  was  designed  for  possible 
remote  site  application.  Specific  design  procedures  were  outlined  for  known  performance 
characteristics.  General  guidelines  were  also  provided  for  designs  when  specific  data 
is  not  available. 

The  period  is  rapidly  approaching  when  various  alternate  energy  schemes  must  be 
considered  to  replace  or  supplement  conventional  sources.  A combined  Darrieus -Savonlus 
system  offers  one  possible  system.  More  concrete  evidence  as  to  exact  performance, 
specifically  the  Savonlus  at  higher  speeds  is  needed  for  a precise  analysis. 


100 


R.  W.  Thrasher  and  R.  W.  Wilson,  "Design  Consideration  for  the  Darrieus  Rotor,* 
Eleventh  Intersocletv  Energy  Conversion  Engineering  Conference  Proceedings.  Vol  II, 


American  institute  of  Chemical  Engineers,  New  York.pp.  1695-1703,  1976. 

E.  G.  Kadlec,  "The  Darrieus  Verticle-Axls  Wind  Turbine  Program,*  Sharing  the  Sun. 
Vol  7,  Winnipeg,  Canada,  p.  351,  1976. 

Department  of  Energy,  "A  Fresh  Interest  in  Wind  Energy,"  The  Sacramento  Bee, 
p.  CL  11-12,  March  4,  1978. 

P.  South  and  R.  S.  Rongi,  "The  Performance  and  Economics  of  the  Verticle-Axis  Wind 
Turbine  Developed  at  the  National  Research  Council,"  Paper  No.  PNW  73-303,  American 
Society  of  Agricultural  Engineers , pp . 14-16,  1974. 

Ibid., Blackwell,  p.  8-9,  1975. 

B.  F.  Blackwell,  The  Vertlcle  Axis  Wind  Turbine  "How  it  Works,1  Sandia  Laboratories 
Energy  Report,  SLA-74-0160,  p.  3,  1974. 

D.  M.  Simmons,  Wind  Power.  Noges  Data  Corp,  Park  Ridge,  NJ.pp.  101-104,  1975. 

R.  J.  Tempi in,  Aero  Dynamic  Performance  for  the  NRC  Verticle-Axls  Wind  Turbine. 
National  Research  Council  of  Canada,  LTR-LA-160,  1974. 

R.  E.  Sheldahl,  B.  F.  Blackwell,  Free-Air  Performance  Tests  of  a 5-Meter  -Diameter 
Darrieus  Turbine.  Sandia  Report,  SAND  77-1063,  1977. 

B.  G.  Newman,  Measurements  on  a Savonius  Rotor  with  Variable  Gap.  Sherbrooke 
University  Symposium  on  Wind  Energy,  1974. 

B.  F.  Blackwell  and  others,  Wind  Tunnel  Performance  Data  for  Two  and  Three 
Bucket  Savonius  Rotors.  Sandia  Energy  Report,  SAND  76-0131,  1977. 

B.  F.  Blackwell,  R.  E.  Sheldahl  and  L.  V.  Feltz,  Wind  Tunnel  Performance  Data  for 
the  Darrieus  Wind  Turbine  with  NACA  0012  Blades.  Sandia  Laboratories  Energy 
Report  Sand  76-0130,  p.  31,  1976. 

Ibid., Blackwell,  p.  8,  1974. 

Ibid., Blackwell,  p.  10-36. 

B.  F.  Blackwell,  and  othera,  "Engineering  Development  States  of  the  Darrieus 
Wind  Turbine,"  Sandia  Energy  Report,  Sand. 

Ibid . .Blackwell,  1977. 


Ibid  , .Blackwell,  p.  31,  1977. 

P.  R.  Nixon,  A Starter  for  th« 
report,  1978. 


Force  Academy  Darrleua  Wind  Turbine,  unpublished 


USAFA-TR-79-1 


References 


101 


B.  F.  Blackwell  and  L.  V.  Feltz,  Wind  Energy  - A Revitalized  Pursuit. Sandia 
Laboratories  Energy  Report  Sand  75-0166,  p.  6,  1975. 


USAFA-TR-79-1 


SECTION  IV 

INSTRUMENTATION  AND  HARDWARE 


102 


USAFA-TR-79-1 


HARMONIC  OPTIMIZATION  OF  A 
PERIODIC  FLOW  WIND  TUNNEL 

•Hr  "fftc 

John  P.  Retelle,  Jr.  and  Donald  A.  Kennedy 


Abstract 


Numerous  methods  have  been  used  to  generate  axial  gustlng  In  wind  tunnels.  Attempts 
to  Increase  the  amplitude  of  velocity  oscillations  about  a mean  velocity  have  usually 
resulted  in  Increasing  the  harmonic  content  of  the  power  spectrum  of  the  free  stream 
velocity.  This  work  describes  a wind  tunnel  modification  which  was  specifically  designed 
to  produce  a sinusoidal  perturbation  to  the  mean  velocity  which  was  free  from  unwanted 
higher  harmonic  effects.  Velocity  perturbations  were  produced  by  rotating  vanes  mounted 
downstream  of  the  tunnel  test  section.  A flow  bypass  mechanism  provided  harmonic 
supresslon  which  could  be  optimized  for  the  various  combinations  of  velocity-oscillation 
amplitude  and  frequency  used  in  unsteady  flow  experiments.  Depending  on  the  flow 
conditions,  results  showed  that  the  second  harmonic  power  amplitude  could  be  greatly 
suppressed  and  that  higher  harmonics  could  not  be  detected. 


I.  Introduction 


Recent  experimental  investigations  of  unsteady  boundary  layer  separation  and  dynamic 
stall  have  used  unsteady  variations  of  airfoil  angle  of  attack  (Ref.  1,  2),  free  stream 
velocity  (Ref.  3,  4,  5), or  both  (Ref.  6).  Of  these  configurations,  experiments  dealing 
with  oscillatory  velocity  perturbations  about  a mean  velocity  have  provided  greater 
ease  of  data  collection  near  the  surface  of  a fixed  airfoil,  while  still  producing  un- 
steady flow  phenomena  resembling  that  noted  during  dynamic  stall  on  an  airfoil  oscil- 

latlns  in  pitch.  These  unsteadv  effects,  such  as  the  presence  of  regions  of  unseparated, 
reversed  flow  in  the  airfoil  boundary  layer,  have  been  observed  (Ref.  5)  to  be  strong 
functions  of  the  frequency  and  amplitude  of  the  oscillatory  free  stream  velocity  per- 
turbations. In  addition,  recent  experimental  data  (Ref.  5)  have  shown  that  the  presence 
of  higher  harmonic  content  in  the  power  spectrum  of  the  oscillating  free  stream  velocity 
produced  large  unsteady  variations  in  the  response  of  the  unseparated  flow  reversal 
regions,  fftus,  free  stream  flow  oscillations  can  be  extremely  useful  in  reproducing 
unsteady  flow  effects.  To  simplify  the  flow  physics,  however,  care  must  be  taken  to 
suppress  or  eliminate  the  second  and  higher  harmonics  (terms  of  factor  2 or  higher  In  a 
Fourier  series  expansion  of  the  free  stream  velocity)  which  are  especially  prevalent  at 
large  oscillation  amplitudes. 
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the  wind  tunnel  modification  described  below  was  specifically  designed  to  produce  a 
large-amplitude  sinusoidal  perturbation  to  the  mean  velocity  which  was  free  from  unwanted 
higher  harmonic  effects.  A set  of  rotating  vanes  at  the  downstream  end  of  the  test 
section  was  designed  to  periodically  block  the  flow  to  the  axial  flow  fan  located  further 
downstream.  A nearly  cons tan  - air  flow  was  provided  to  the  fan  by  a bypass  mechanism 
which  allowed  the  fan  to  operate  at  a nearly  constant  pressure  rise,  improving  the 
harmonic  suppression  in  the  flow.  It  was  found  that  this  harmonic  suppression  could  be 
optimized  for  the  various  combinations  of  velocity  oscillation  amplitude  and  frequency 
used  in  other  experiments. 

Several  experiments  in  wind  tunnels  have  used  the  rotating  vane  principle  to  perturb 
the  flow.  Two  similar  tunnels  made  no  attempt  to  Improve  the  suppression  of  higher 
harmonics.  An  apparatus  on  an  open-circuit  wind  tunnel  described  by  Despard  (Ref.  7) 
used  four  rotating  vanes  downstream  of  the  test  section  to  produce  velocity  variations 
of  up  to  40  percent  at  oscillation  frequencies  of  up  to  240  Hz.  The  second  and  third 
harmonics  contributed  23  percent  and  4 percent  of  the  velocity  power  spectrum,  respec- 
tively. An  almost  identical  apparatus  was  described  by  Morkovin,  Loehrke,  and  Fejer 
(Ref.  8)  which  perturbed  the  flow  in  a closed-circuit  wind  tunnel.  They  found  resonances 
which  corresponded  to  the  tunnel  circuit  length,  the  distance  between  the  rotating 
vanes,  and  the  tunnel  contraction.  They  also  noted  large  departures  from  slnusoidality 
associated  with  large  velocity  oscillation  amplitudes.  Recent  work  in  this  tunnel  by 
Saxena,  Fejer,  and  Morkovin  (Ref.  4)  produced  plots  of  the  velocity  and  pressure 
fluctuations  which  indicated  a higher  harmonic  content  in  the  pressure  fluctuation  than 
in  the  longitudinal  velocity  variation. 


Several  different  schemes  have  been  tried  to  eliminate  higher  harmonic  content  in 
the  velocity  oscillations  of  these  rotatlng-vane  tunnels.  Karlsson  (Ref.  9)  used  four 
rotating  vanes  downstream  of  the  test  section  of  an  open-circuit  wind  tunnel  to  produce 
root-mean-square  velocity  variations  of  34  percent  of  vane  rotation  frequencies  of  up  to 
48  Hz.  At  high  frequencies  the  higher  harmonics  were  responsible  for  about  15  percent 
of  the  total  root-mean-square  amplitude  of  the  fluctuation.  In  an  effort  to  reduce 
these  harmonics,  the  lowest  vane  was  turned  at  a higher  rate,  but  no  mention  is  made 
of  the  success  of  this  device.  Simpson  (Ref.  10)  described  an  apparatus  of  this  type 
which  used  five  rotating  vanes  downstream  of  the  test  section  of  an  open-circuit  tunnel. 
A programmable  motor  was  used  to  adjust  the  angular  velocity  of  the  rotating  vanes, 
thereby  optimizing  the  higher  harmonic  suppression.  The  resulting  velocity  waveforms 
appeared  sinusoidal,  but  no  harmonic  analysis  was  presented. 


An  interesting  variation  to  the  rotating  vane  device  was  described  by  Chamay  and 
Mathieu  (Ref.  11).  Their  eight  rotating  vanes  were  positioned  upstream  of  the  test 
section  in  the  settling  chamber  ahead  of  the  screens  and  honeycomb  sections  of  the  open- 
circuit  wind  tunnel.  They  found  that  the  smallest  harmonic  content  occurred  when  the 
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vane  vortex-shedding  frequencies  matched  the  tunnel  resonant  frequencies.  Turbulent 
velocities  were  measured  to  evaluate  the  effect  of  the  upstream-mounted  vanes.  They 
observed  that  optimum  sinuaoidality  was  indicated  by  maxima  of  the  magnitudes  of  the 
longitudinal  and  vertical  turbulent  fluctuating  velocities,  ui ' and  u^',  and  by  minima 
of  the  skewness  and  flatness  factors  of  ui’.  They  measured  the  turbulence  intensity  in 
the  test  section  downstream  of  the  vanes  to  be  on  the  order  of  0.7  percent. 

The  rotating  vane  apparatus  described  below  was  capable  of  higher  harmonic  suppression 
to  a greater  degree  than  the  devices  described  above,  with  the  possible  exception  of  the 
device  described  by  Simpson  (Ref.  10).  Careful  adjustment  of  the  bypass  mass  flow 
vlded  regulation  of  the  higher  harmonics  in  the  test  section  for  a given  tunnel  mean 
velocity  and  vane  rotation  frequency.  Depending  on  the  vane  size  and  flow  conditions, 
the  second  harmonic  power  amplitude  could  be  suppressed  as  much  as  48  db. 


II.  Experimental  Apparatus 

The  unsteady  flow  device  was  installed  in  the  two  foot  by  two  foot  (0.61  meter  by 
0.61  meter)  cross-section,  low-speed  wind  tunnel  at  the  University  of  Colorado,  Boulder. 
It  was  an  open-return,  suction-type  wind  tunnel  with  the  addition  of  the  unsteady  flow 
apparatus.  Figure  1 shows  a schematic  of  the  wind  tunnel.  The  test  section  was  14  feet 
(4.27  meters)  in  length  with  a 10  foot  (3.05  meter)  long  plexiglass  wall,  % inch  thick. 
Downstream  of  the  test  section  were  the  unsteady  flow  device,  a transition  section  to 
the  fan,  the  fan,  and  the  conical  diffuser. 


A periodic  disturbance  to  the  mean  flow  was  produced  by  a set  of  six  rotating  vanes 
located  just  downstream  of  the  test  section.  Figure  2 shows  a schematic  of  the  vane 
drive  apparatus.  Three  sets  of  vanes  of  varying  width  were  used  during  this  experiment 
to  change  the  amount  of  maximum  tunnel  blockage  and,  thus,  the  velocity  amplitude  variation, 
N.  The  vanes  were  made  in  a modified  diamond  shape,  both  for  ease  of  construction  and 
for  minimal  pressure  and  friction  losses.  Of  the  three  sets  of  vanes,  set  A used  a 
2.000  inch  chord,  set  B used  a 3.495  inch  chord,  and  set  C used  a 3.980  inch  chord. 

The  three  vane  sets  produced  tunnel  area  blockages  of  50,  87.5,  and  99.5  percent, 
respectively.  The  rotating  vanes  were  driven  by  a gear  train  which  was  in  turn  powered 
by  a General  Electric  Stratltrol  II  3/4  horsepower,  full-wave,  direct-current  motor 
controlled  by  a Variac-equlpped, direct-current  power  supply.  An  8-to-l  gear  reduction 
allowed  the  motor  to  run  smoothly  near  the  middle  of  its  0-1759  revolutions-per-minute 
speed  range  preventing  surging  at  low  speeds.  The  gear  drive  was  designed  to  allow 
the  vanes  to  counterrotate  and,  thu%  produced  no  net  local  lift  perturbation  to  the 
flow. 

The  bypass  mechanism  allowed  air  to  enter  the  fan  section  when  the  tost  section  was 
blocked  by  the  front  (free  stream-controlling)  vanes.  Three  rotating  vanes  on  the  top 
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Figure  1.  Schematic  of  Wind  Tunnel 
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Figure  2.  Teet  Section  Schematic 
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of  the  device  and  three  on  the  bottom  were  driven  by  the  same  gear  train  which  drove 
the  front  vanes.  The  vanes  on  the  top  and  bottom  allowed  air  from  the  laboratory  to 
enter  tne  fan  section  at  a varying  rate  which  balanced  the  decrease  flow  in  the  test 
section.  The  adjustment  of  the  bypass  air  flow  into  the  fan  section  was  controlled  by 
moveable  slides  which  are  shown  in  the  test  section  schematic  of  Figure  2.  These  slides 
were  set  to  predetermined  bypass  openings  based  on  the  flow  conditions  desired  in  the 
test  section.  Unlike  the  front  vanes,  the  bypass  vanes  had  a fixed  2-lnch  (5.08  cm) 
chord  (vane  set  A). 


The  free  stream  steady  turbulence  level 


• (A  j=.y- 


for  the  basic  wind 


tunnel  was  measured  by  Francis  (Ref.  12)  at  a nominal  value  of  0.06  percent  for  a free 
stream  velocity  of  60  feet  per  second  (18.29  meters  per  second).  In  this  expression, 
ui  are  the  turbulent  velocity  components , and  is  the  magnitude  of  free  stream  velocity. 
Turbulence  intensity  measurements  with  the  rotating  vane  device  installed  were  made 
using  a Bruel  and  KJaer  type  2417  Random  Noise  Voltmeter  which  had  a bandwidth  of 
2-20,000  Hz.  With  the  large  front  vanes  (set  C)  fixed  at  the  minimum  resistance  (full 
open)  position,  the  steady  flow  turbulence  level  was  measured  to  be  0.08  percent  at  a 
free  stream  velocity  of  60  feet  per  second  (18.29  meters  per  second).  The  vanes  were 
then  rotated  at  3.33  Hz  and  the  turbulence  level  was  measured  to  be  0.54  percent  at  a 
mean  flow  velocity  of  12  feet  per  second  (3.66  meters  per  second).  To  determine  the 
effect  of -unsteadiness  on  the  higher  frequency  turbulence,  this  velocity  data  signal 
was  filtered  through  a 4 Hz  high-pass  filter.  The  higher  frequency  root-mean-square 
turbulence  level  was  then  measured  to  be  0.05  percent. 

III.  Unsteady  Flow  Field  Analysis  and  Optimization 

Results  of  an  analysis  by  McMichael  (Ref.  13)  show  that  a tunnel  of  this  type  could 
be  optimised  to  suppress  flow  harmonics  and  could,  in  certain  configurations,  eliminate 
all  even  harmonics.  The  applications  of  this  analysis  are  described  below  as  well  as 
a description  of  the  mechanical  optimization  using  the  unsteady  flow  device. 

The  analysis  of  McMichael  (detailed  in  the  Appendix)  showed  that  the  smallest 
temporal  change  in  the  pressure  rise  across  the  fan  occurred  when  the  pressure  and 
frictional  losses  across  the  bypass  ports  and  vanes  were  equal  to  the  sum  of  the  pres- 
sure drop  across  the  main  tunnel  vanes  and  the  frictional  loses  in  the  test  section. 
Careful  adjustment  of  the  bypass  port  openings  permitted  a fine  tuning  of  the  bypass 
frictional  losses  and  permitted  an  almost  steady  flow  to  reach  the  fan.  The  resulting 
unsteady  free  stream  velocity  was  of  the  form 


U«.(x,t)  - Ui(x) 


N £ [an 

n“l  L 


cos  (n </£)  + b sin  (nut) 
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where  Ux  (x)  Is  Che  mean  velocity  (as  a function  of  axial  distance  x)  and  N Is  the 
oscillation  amplitude.  The  analysis  showed  that,  once  the  free  stream  velocity  and 
oscillation  frequency  were  set  and  the  resulting  oscillation  amplitude  measured,  the 
bypass  flow  could  be  varied  to  suppress  all  the  coefficients  in  equation  (1)  except  bx. 

In  addition  to  the  harmonic  suppression,  McMlchael  has  shown  that  the  even  harmonics  of 
the  flow  oscillation  can  be  eliminated  when  the  main  flow  vanes  and  bypass  flow  vanes 
are  identical. 

The  goal  of  the  unsteady  flow  analysis  was  to  provide  a longitudinally  oscillating 
flow  at  a desired  amplitude  ratio,  oscillation  frequency,  and  mean  velocity  whose  power 
spectrum  was  as  free  as  possible  of  harmonic  content.  For  the  experiment,  of  the  four 
controllable  tunnel  parameters,  a combination  of  vane  size,  steady  tunnel  velocity,  and 
vane  rotation  frequency  (or  rotation  period)  were  selected  to  produce  the  desired 
frequency  and  oscillation  amplitude  ratio.  The  upper  and  lower  fixed  bypass  openings, 
each  at  a setting  of  L/2  (see  Figure  2),  were  then  adjusted  for  maximum  second  harmonic 
suppression,  according  to  the  parameters  suggested  by  the  analysis.  The  variable  bypass 
provided  a more  uniform  flow  to  the  fan  over  the  period  of,  and  In  correlation  to,  the 
vane  rotation. 

IV.  Results  of  Unsteady  Flow  Optimization 

Mean  and  oscillating  velocities  were  measured  in  the  test  section  using  a single- 
wire hot-wire  anemometer.  The  linearized  output  voltage  was  ensemble-averaged  over 
50  data  entries,  with  the  start  of  each  data-collectlon  cycle  being  triggered  by  a 
magnetic  pickup  mounted  on  the  vane  gear  drive.  A power  spectrum  of  the  resulting 
signal  was  produced  using  a Fourier  analyzer.  A typical  power  spectrum  of  the  linear- 
ized velocity  signal  for  the  small  vanes  (set  A)  is  shown  In  Figure  3 for  a non- 
dimensional  frequency  of  k-0.873  and  a bypass  opening  of  (L/2)*l  inch  (2.54  cm). 

The  second  harmonic  power  was  observed  to  be  approximately  48  db  lower  than  the  power 
of  the  fundamental.  The  absence  of  odd  harmonics  was  difficult  to  confirm  experimentally 
since,  for  the  small  vanes,  usually  only  the  second  harmonic  was  large  enough  to  be  de- 
tected by  the  Fourier  analyzer  on  a logrlthmlc-magnitude  power  spectrum  (maximum  scale  of 
80  db).  The  db  differences  between  the  first  and  second  harmonics,  20  log^g  (V^/V^) , are 
plotted  in  Figure  4 for  a variation  of  L/2  at  this  flow  condition.  For  vane  set  A,  at  a 
mean  velocity  of  51.6  feet  per  second  (15.72  meters  per  second),  the  optimum  or  maximum 
power  difference  Is  observed  to  be  approximately  48  db  and  the  smallest  difference  Is  32 
db.  The  power  differences  between  the  fundamental  and  third  harmonic  and  between  the 
fundamental  and  fourth  harmonic  were  also  noted,  If  measurable.  In  no  case  were  they 
measured  to  be  less  than  60  db  and  70  db  respectively  for  the  small  vanes. 

The  power  spectrum  results  for  the  larger  vane  sets  were  acceptable,  but  not  as 
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free  from  harmonics  as  those  for  the  smallest  vanes.  The  power  differences  between  the 
fundamental  and  second  harmonic  as  functions  of  L/2  are  presented  in  Figure  5 for  both 
the  medium  vanes  (set  B)  at  a mean  velocity  of  31  feet  per  second  (9. 45  meters  per 
second)  and  for  the  large  vanes  (set  C)  at  a mean  velocity  of  12  feet  per  second  (3.65 
meters  per  second).  At  these  flow  conditions,  the  medium  and  large  vane  sets  were 
observed  to  produce  maximum  harmonic  differences  of  15  db  and  11  db,  respectively.  The 
maxima  can  be  clearly  distinguished,  but  they  are  not  as  large  in  magnitude  as  those 
measured  with  the  small  vane  set  (set  A)  Installed.  It  must  be  remembered  that,  for 
mechanical  simplicity,  only  the  main-flow  vanes  were  replaced  when  changing  vane  seta* 
the  bypass  vanes  were  not  replaced. 

The  experimental  results  show  that,  even  with  different  sixes  of  main  and  bypass 
vanes,  a maximum  harmonic  suppression  can  be  achieved  with  the  proper  setting  of  L/2. 

In  order  to  obtain  the  parameters  required  for  the  flow-reversal  experiments,  the 
optimum  setting  of  L/2  could  not  always  be  used.  The  objectives  of  the  flow -reversal 
experiments  required  two  separate  parametric  variations  of  the  flow.  The  first  varied 
the  reduced  frequency  and  held  the  mean  velocity  and  oscillation  amplitude  nearly 
constant.  The  second  varied  the  oscillation  amplitude,  with  the  reduced  frequency  and 
mean  velocity  almost  constant.  It  was  observed  that  the  magnitudes  of  the  mean  velocity 
and  amplitude  ratio  varied  considerably  with  L/2  for  a given  frequency.  Even  though 
off-optimum  settings  of  L/2  were  used,  the  second  harmonic  rejection  was  at  least  10  db 
for  all  cases. 

The  required  combination  of  oscillation  frequency  and  amplitude  are  usually  dictated 
by  experimental  needs  and  not  by  harmonic  optimisation  of  the  flow  field.  The  para- 
meters of  the  unsteady  flows  produced  by  each  vane  set  are  shown  in  Table  1. 


TABLE  1 

FLOW  PARAMETERS  - FLOW  REVERSAL  MEASUREMENTS 


Vane 

L/2 

u0 

T 

Vi 

N 

R 

20  log 

Set 

in 

ft/sec 

•ec 

ft/sec 

7,  _ 

k 

xlO® 

[ViAfc)10 

A 

0.5 

60 

0.5 

51.6 

5.65 

0.101 

2.421 

47 

45 

0,5 

36.6 

5.15 

0.143 

1.717 

49 

B 

7.0 

60 

0.3 

29.5 

37.5 

0.296 

1.383 

18 

0.5 

31.1 

40.0 

0.168 

1.459 

16 

0.7 

31.3 

47.0 

0.119 

1.470 

10 

3.0 

40 

O.f 

12.2 

59.0 

0.429 

0.570 

13 

C 

4.5 

60 

0.5 

11.6 

98.0 

0.453 

0.542 

u — 

As  an  aside,  it  was  found  that  harmonic  optimisation  could  also  be  confirmed  using 
the  method  of  Charnay  and  Mathleu  (Raf.  11)  who  measured  the  skewness  and  kurtosls 
of  the  longitudinal  turbulent  velocity  component,  ui ',  in  an  oscillatory  flow  tunnel. 
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They  found  the  conditions  where  S approached  zero  and  K approached  1.5  represented  an 
optimum.  Figures  6 and  7 show  the  values  of  K and  S versus  L/2  with  the  small  vanes 
(set  A)  for  the  same  flow  conditions  as  considered  above.  It  can  be  seen  that  the 
optimum  values  of  L/2  are  the  same  as  predicted  by  the  spectral  analysis.  Also,  the 
points  of  maximum  S and  K correspond  to  the  points  of  minimum  harmonic  suppression.  The 
zero  value  of  S Indicates  that  the  distribution  of  the  velocity  harmonics  Is  Gaussian 
and  that  the  disturbance  caused  by  the  vanes  only  affects  the  pressure  field.  The  low 
value  of  K shows  that  the  uj ' distribution  Is  broad  for  the  optimized  tunnel  condition. 

V.  Conclusions 

The  higher  harmonic  content  of  the  free  stream  velocity  oscillations  was  reduced 
by  adjustment  of  the  bypass  flow  of  the  rotating  vane  device.  The  best  operating 
conditions,  with  a high  reduced  frequency  and  a small  amplitude  ratio,  provided  a 
second  harmonic  rejection  of  approximately  48  db.  Other  combinations  of  flow  parameters 
used  for  other  experiments  produced  second  harmonic  rejections  of  at  least  10  db.  It 
Is  recommended  that  a feedback  control  system  for  the  vane  motor  drive,  similar  to  that 
used  by  Simpson, be  incorporated  to  improve  the  harmonic  rejection  capability  of  this 
device. 
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Symbols 


A(x) 

Fl 

Ve 

K 

Kmn  ^ 

t/2 

N 

Q(x,t) 

Qi 

Q* 

R 

Ri 

R 

» 

R' 


S 

u. 

u 

o 

Vi 


V 

a 

<v> 

c 

fi(t) 


V#) 

fi(t) 

f#  (t) 

k 

qx(t) 


wind  tunnel  cross-section  area,  as  a function  of  longitudinal  distance 
pressure  drop  coefficient  across  main-flow  rotating  vanes 
pressure  drop  coefficient  across  bypass-flow  rotating  vanes 
kurtosls  factor  of  ui',  i.e.,  (u^*  4 )/(ui' 8 
pressure  changes  between  wind  tunnel  stations  m and  n 

width  of  top  or  bottom  bypass  port  opening  in  the  unstead 

Au  ( 

oscillating  velocity  amplitude  ratio,  l.e.,  N ■ * 

mass  flow  per  unit  volume,  A(x)u(x,t) 
steady-flow  mass  flow  per  unit  volume 


V flow 
^maximum 


rator 

-UM 


mass  flow  per  unit  volume  through  main-flow  rotating  vanes 
mass  flow  through  bypass  rotating  vanes 

Reynolds  number,  ,^^>c 

sum  of  coefficients  of  viscous  and  pressure  losses  for  flow  through  tunnel 
test  section  and  main-flow  vanes 

sum  of  coefficients  of  viscous  and  pressure  drop  losses  for  flow  through 
bypass  ports  and  bypass  vanes 

coefficient  of  viscous  and  pressure  losses  for  Identical  main- flow  and  bypass 
rotating  vanes 

-Try  3/2 

skewness  factor  of  ui',  l.e.,  (u/  3 ) / (u/  x ) 

free  stream  velocity  (steady  or  oscillatory) 

free  stream  velocity  prior  to  operating  unsteady  flow  device 

mean  value  of  oscillating  longitudinal  free  stream  velocity 

voltage  used  to  compute  the  power  contained  in  the  first  harmonic  of  the 
velocity  oscillation  power  spectrum 

voltage  used  to  compute  the  power  contained  in  the  second  harmonic  of  the 
velocity  oscillation  power  spectrum 

average  voltage  resulting  from  the  ensemble  average  process 
airfoil  chord  length 

unsteady  component  of  pressure  drop  coefficient  across  main-flow  rotating 

vanes 

unsteady  component  of  pressure  drop  coefficient  across  bypass  rotating  vanes 

nond lawns lonal  value  of  fi(t),  l.e.,  f1/2R' 
nondlmenalonal  value  of  f^(t),  l.e.,  ^/^R* 

reduced  frequency,  !— 

2Ui 

unsteedy  component  of  Qi(x,t) 
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% (t) 


t 

«i 

us 

x 


P 

y(x) 

Ap„ 


Apr 


p 

p 

Tl 


7b 


01 

*2 

f) 

<(  )> 


unsteady  component  of  (x, t) 

ratio  of  Ri  /R  «*  (Q  /Q  ) ^ 

2 *o  10 

time 

velocity  component  in  longitudinal  direction 
velocity  component  in  vertical  direction 

length  dimension  along  longitudinal  (axial)  coordinate 

position  vector,  |x|  “ J if  + y*  * r* 

mass  flow  passing  between  wind  tunnel  stations  m and  n 

coefficient  of  fan  pressure  drop 
coefficient  of  viscous  losses 
steady  pressure  drop  across  fan 

unsteady  pressure  drop  across'  fan 

coefficient  of  viscosity  of  air 
air  density 

time  constant  associated  with  area  changes  upstream  of  test  section,  l.e., 

ai/RQ0 

time  constant  associated  with  area  changes  downstream  of  test  section,  i.e., 
Ob/RQ0 

vane  rotation  angle  of  main-flow  vane 
vane  rotation  angle  of  bypass-flow  vane 

denotes  a time  average 
denotes  an  ensemble  average 
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APPENDIX 

UNSTEADY  FLOW  ANALYSIS 

The  following  Is  an  analysis  of  an  unsteady  flow  wind  tunnel  which  has  been  adapted 
for  the  tunnel  configuration  described  in  Section  II.  Following  the  discussion  of 
McMichael  (Ref.  5),  the  momentum  equation  for  the  wind  tunnel  flow  may  be  written. 


£>u  bu 

bt  U bx 


(Al) 


where  y(x)  represents  losses  due  to  the  friction  of  the  walls.  The  mass  balance  Is: 


Q(x,t)  - A(x)  u(x,t) 


(A2) 


where  A(x)  is  the  local  cross-sectional  area.  Substituting, 


L + a!  ±t  i \ - - -S-  (£\ 

A(x)  bt  2 bx\W0i))  A*  (x)  bx  \p) 


(A3) 


The  wind  tunnel  stations  and  mass  flows  are  defined  in  Figure  1.  Integrating  between 
arbitrary  stations  m and  n yields: 


(P  - P„)  - a I?  + K 

m n mn  ot  in 


(A  4) 


where: 


Knowing  that  P * P “Pa  (atmospheric  pressure)  and  assuming  quasi-steady,  high 

1 s 

Reynolds  number  flow,  the  pressure  differences  at  each  station  are: 


P-P-O  Q+K  Q * 
x • 1*1  i*  l 


P-P-O  Q+K  Q 8 
3 4 34  1 34  X 


P 

5 


P - 

6 


BC 


(Qx  +S 


) + K (Q  + 


66 


(A5) 
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P„  ' P " (Q,  + 9.  > + K <Q  + f 

1 a 78  1 * 7e  l 2 


P - P = C Q * 

2 3 83  1 


P - P = C Q * 

4 5 48  1 


P - P = C.  Q * 

4 B • B 8 


where  Cmn  Is  the  pressure  drop  coefficient  across  the  indicated  mechanical  devices. 


P - P = - P„  (fan  pressure  rise), 

6 7 I 


Then,  McMichael  writes: 


(P  • P ) « (P  - P ) - (a  + a ) Q + (K  + K + C + C ) Q 8 , 

18  “ 8 18  34  1 18  34  83  4S  1 


Now  define 


a » (a  t-  a ) 

1 12  34 


K - (K  + K + C ) 

1 18  34  *3 


F * C 

1 4B 


F « C 

t *8* 


Combining  these  terms  gives 


F Q * - a Q + (K  + F )Q  8 
8 t 11  1 11 


The  adjustable  bypass  ports  can  be  opened  a distance  L and  essentially  add  a resistance 
(L)  so  that 
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P - P » (K  + F )0  * . 

a b 1 2 11 


(P  - P ) - (K  + F ) Q * = a Q + (K  + F )Q  * . (All) 

A g 2 8*  11  1 11 

Mote  that  O Q was  considered  negligible  because  o£  the  short  distance  between  the  by- 
2 B 

pass  port  and  bypass  vanes.  For  the  fan  duct  exit 


(P  - P ) - (P  - P)  + (P  - P)  + (P  - P) 

s • 6 6 e 7 7 e 


(a  + 0 ) (Q  +Q)  + (K  + K ) (Q  +Q)8-p<: 

B«  7.  1 2 66  76  1 I 1 


define 


K - K + K 

3 66  7e 


0-0  + o . 

3 66  7a 


Collecting  the  previous  relations  yields 

(K2+F2>V 


Pa  Ps 


(K  + F ) Q f + O Q 
1 11  11 


'Ks(Qi  + " * % + V + Pf’ 

A reference  state  for  the  tunnel  flow  way  be  defined  using  the  tine  average  of  the  vane 
pressure  drop  coefficients  and  mass  flows  at  a given  naan  tunnel  velocity.  Then, 


F (t)  - F 

+ f (t) 

1 

10 

l 

F (t)  - F_ 

+ f (t) 

2 

lo 

t 

Q (t)  - Q 

+ q (t) 

i 

10 

1 

V°-So 

+ <^(t) 
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F (t)  - F 
1 10 


F (t)  - F 

S to 


Q (t)  - Q 
l 10 


f (t)  - f (t)  « q (t)  - q (t)  • 0 
12  1* 


P ■ p + Ap 
f f f 

o 


q Q A? 

— 1-  ~ -J — .. £.  < < 1 

Q Q pf 

io  to  £ 


The  reference  state  (tine  average)  of  equation  13  shows 

(K  + F )0  * - (K  + F )Q  * ..... 

* to  l 10  10  (A14) 

(K  + F )Q  8 - P.  - K (Q  + Q P . 

1 lo  io  3 ^o'  (A15) 

Equations  (2),  (3),  and  (4)  are  than  linearised  about  the  reference  state  (second* 
order  terns  dropped): 

2<K.  * F.A.  % * '.V*  ’ 2«, + W * *,«„*  ♦ \ «“> 

2<K.  * V V * t,\‘ m Ltt  ■ V\  * V ‘ V..  * ♦ %>•  (U7> 


For  the  slnpllfled  case  where  the  main- flow  and  bypass  vanes  are  Identical,  equations 
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2 §_  ,*  . * x 2a  a ,*  A 2Ka  .A  . a , 

RQ0  (qi  V RQ„  (qi  ' V " R (qi  V 


(A20) 


Now  define  the  tine  constants 


and  also  let 


T T -fib- 

Ti  RQ0’  T3  RQg 


. -J_  . SflL. 


R * 


Substituting  gives 


T T q + [r  (1-G)  + 2t  ]q  + (l-2G)q 
13  1 l 3 1 \ 

• t 

A A A A 

“G(£1- V-Vfi  V* 


The  damping  ratio  for  this  second  order  system  la 


(A21) 


(l-G)r  + 2t 

{ . — . -1- A 

2 7(1-G)t  T 

1 3 


(A22) 


end  the  resonant  frequency  is 


% 


JUZG  _ ^/^^2G 


1 3 


-J  a a 

1 3 


(A23) 


It  can  also  be  shown  that  these  smell  vanes  of  set  A ideally  generate  only  even  harmonics 
in  the  velocity  field.  Equation  (21)  shows  that  the  forcing  function  Is  (f  - £ ). 
McMlchael  shows  the  angular  dependence  between  the  two  sets  of  vanes  to  be 


20  - wt+  ir/2 
i 


20^  - wt  - W/2. 
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For  a given  vane  resistance  function 


f(0)  “ a cos  (2n0) 

l . n 

n"l 


It  can  be  seen  that 


a a * y 
f (0  ) + f <0. ) - -£  + L 
11  e t 2 


n-0,2,4,... 


an  cos  (n(wt  + ff/2)j 


f (0  ) - £ (0  ' 

11  £ £ 


n* 1|3|5| • • • 


•n  cos  (n(wt  - ff/2)) 


The  forcing  function  (f  - f ) in  the  dynamic  equation  thus  Ideally  only  contains  odd 

l S 

harmonics.  This  was  difficult  to  determine  experimentally  since,  for  the  small  vanes, 

usually  only  the  second  harmonic  was  of  a large  enough  amplitude  to  be  detected  by  the 

Fourier  analyzer  on  a logrlthmlc-magnltude  power  spectrum  (maximum  scale  of  80  db).  The 

equations  cannot  be  so  easily  simplified  for  the  situation  where  the  main-flow  vanes 

and  bypass  vanes  are  different,  because  F it  and  K 4 K . The  first  line  of 

10  fo  i I 

equation  (9)  may  then  be  written  (Including  the  assumption  for  Ap{): 


S.  - A'f  - B'f  + C'i  + D'q 
1 ill 


A'  - Q / (2rR  ) 
10  £ 


B1  - rQ  / (2R  ) 
10  • 


c'  - Of  /(2rR  Q ) 
1 £ lo 


D'  " r. 


Substituting  expression  (B-ll)  Into  the  second  line  of  equation  (B-8)  yields: 

1,  rVi  Mh  (1+  2^u-n>0.„-a  \ . »«.  i 

*•  \ *v..  ) “J 


J 

*4 
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CONFIDENCE  IN  FUNDAMENTALS 
* 

Roger  W.  Gallington 


Abstract 


The  technical  performance  of  recent  engineering  graduates  may  be  limited  by  a lack  of 
confidence  in  fundamentals.  This  contention  is  supported  here  by  three  real  engineering 
problems  requiring  only  undergraduate  knowledge  for  their  solutions,  which  the  author 
doubts  many  recent  graduates  could  solve  with  sufficient  confidence  to  direct  the  activ- 
ities of  others  to  implement  the  solutions.  A means  of  building  confidence  during  edu- 
cational programs  by  systematic  comparisons  between  different  approaches  to  the  same 
problem  is  suggested. 

I,  Introduction 


Many  engineering  educators  comment  that  the  graduates  are  increasingly  unable  to 
solve  the  technical  problems  they  face  after  leaving  the  educational  institution.  I 
believe  there  are  two  reasons  for  this  trend.  The  first  reason  is  that  faculties  are 
becoming  more  isolated  from  Industry  and,  consequently,  the  new  graduate  simply  does  not 
have  an  appropriate  background.  The  second  reason  is  that,  because  of  the  lack  of 
specific  experiences  in  his  education,  the  new  graduate  does  not  have  sufficient  confi- 
dence in  the  fundamentals  of  his  field.  This  paper  addresses  the  second  reason. 

K» 

A lack  of  confidence  in  fundamentals  may  be  limiting  the  technical  performance  of  our 
graduates . 

To  illustrate  the  Importance  of  specific  experiences  in  building  confidence,  I will 
describe  three  practical  engineering  problems  I have  solved  since  my  assignment  to  the 
Academy.  None  of  these  problems  is  especially  difficult  and  none  requires  technical 
information  beyond  that  taught  in  the  typical  modern  undergraduate  Aeronautical  Engineer- 
ing program.  Confidence  in  these  solutions  is  based  on  specific  past  experiences.  As 
you  read  through  these  examples  ask  yourself  what  fraction  of  recent  graduates  could 
correctly  work  these  problems  without  assistance  and  be  confident  enough  of  their  answers 
to  direct  the  activities  of  technicians  in  implementing  the  solutions.  I have  pondered 
this  question  and  think  the  fraction  must  be  very  low.  Herein  lies  a major  challange  in 


a 
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engineering  education.  How  do  we  lncreeee  the  students'  confidence  In  their  fundamental 
technical  education?  After  the  example  problem* , suggestions  for  building  technical 
confidence  are  made. 

II.  Some  Real  Engineering  Problems 

These  three  problems  are  real  In  every  sense  of  the  word.  They  are  facilities 
oriented  simply  because  facilities  happen  to  be  my  present  responsibility.  However, 
similar  practical  problema  occur  In  all  types  of  engineering  Including  the  most  sophisti- 
cated research  and  development.  In  fact.  Inability  to  correctly  and  quickly  solve 
problems  of  this  sort  can  cause  major  and  unnecessary  slippage  In  programs  Intending  to 
address  much  more  challenging  technical  Issues. 


A.  Turning  Vane  Problem 

A bearing  failed  In  the  subsonic  wind  tunnel.  During  the  repair  operation  It  was 
necessary  to  cut  a turning  vane  (several  of  these  are  located  In  each  corner  of  a wind 
tunnel  circuit)  In  two  places.  It  was  obvious  to  the  technicians  that  future  repairs 
would  be  easier  If  during  the  re-assembly  the  turning  vane  was  bolted  together  Instead  of 
being  welded.  Welding  would  return  the  vane  to  Its  original  strength  whereas  the  bolted 


Figure  1.  Hind  Tunnel  Geometry 
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v«ne  would  b«  weaker,  the  question  for  the  engineer  was  whether  to  allow  the  bolted 


repair  or  to  require  welding 


1.  Given: 


The  geometry  of  the  relevant  part  of  the  wind  tunnel  and  needed  dimensions  were  all 
available  from  the  construction  drawings  shown  In  Figure  1.  The  maximum  achievable 
velocity  In  the  test  section  le  400  /sec,  and  the  density  of  air  at  the  Air  Force  Academy 
Is  about  0.0018  8lg/fts  . 


The  air  load  on  the  turning  vane 


Figure  2.  Control  Volume 


3.  Solution 


Draw  the  control  volume  shown  In  Figure  2.  Assume  that  the  row  of  turning  vanes 
approximates  a perfect  two-dimensional  cascade  so  that  the  pressures  on  the  two  curved 
surfaces  of  the  control  volume  balanee  eech  other.  Further,  assume  that  the  flow 
la  Incompressible,  uniform  across  the  Inlet  and  outlet  arees  of  the  control  volume,  and 
attached  to  the  tunnel  wells  and  turning  vanes.  The  continuity  equation  was  applied 
between  station  1 end  2 and  stations  2 and  3 yielding: 
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Bernoull's  equation  was  applied  to  the  flow  through  the  turning  vanes,  showing  that  the 
pressures  on  the  Inlet  and  outlet  areas  were  equal  and, therefore,  balanced.  The  two 
remaining  components  of  the  momentum  equation  became 


Fx  - p V,2  A, 


»/  ‘ o. 

/n  + 1 


0018  x (121. 88)2  x *=  40.50  lbf 


P-pVW  -40. 
y /n  + 1 


50  lbf 


where  n Is  the  number  of  turning  vanes  and  n + 1 Is  the  number  of  channels. 


The  total  force  becomes 


FTOT  “ <V  + Fxa>%  “ 57-27  lbf  • 


Since  the  load  was  assumed  to  be  uniformly  distributed,  the  loading  could  be  described  as 
that  of  a uniformly  loaded  beam  with  rigid  ends  as  shown  in  Figure  3. 


12.91  LBF/ FT 


4 FT  5jIN 


28.64  LBF 


28.64  LBF 


Figure  3.  Load  on  Turning  Vane 
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Having  found  the  load  to  be  small,  approval  was  given  for  the  bolted  repair.  After 
the  turning  vane  was  re-installed  and  before  the  tunnel  was  operated,  the  total  calcu- 
lated load  was  applied  to  the  center  of  the  vane.  (This  loading  condition  is  more  severe 
than  the  calculated  one).  The  vane  deflection  was  negligible. 

The  tunnel  was  carefully  run  up  to  its  top  speed,  verifying  both  the  bearing  repair 
and  the  adequacy  of  the  bolted  turning  vane. 


B.  Wind  Tunnel  Blow-Down  Problem 


The  two  large  (200  horsepower,  600  psig)  air  compressors  which  normally  supply  air 
to  the  supersonic  wind  tunnel  were  being  repaired.  Due  to  contractural  and  funding 
difficulties  it  was  obvious  that  they  would  be  Inoperative  for  some  time.  The  tunnel  is 
of  the  blow-down  type  (meaning  that  the  compressors  charge  large  air  storage  tanks  and 
the  tunnel  is  subsequently  run  by  rapidly  depleting  the  stored  air).  The  question  was: 
Could  the  tunnel  be  operated  to  complete  two  cadet  projects  during  the  current  semester 
by  using  a five-horsepower  compressor  capable  of  delivering  air  at  110  psig?  Suprising- 
ly  enough  the  answer  was  yes.  The  details  follow. 


STORAGE  TANK 


HEAT  SINK 

CONTROL  VALVE, 
■fit  rAND  SCREENS  / 


STILLING  CHAMBER 


TEST  SECTIPV 


V - 5281.3  FT 

M - 2.02 


T , - 530  °R 
si  2 

A - 1.0  FT 


33  PSIA 


Figure  4.  Blow-Down  Tunnel  Schematic 
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1.  Given: 

A schematic  of  the  parts  of  the  system  essential  to  the  solution  is  shown  in  Figure  4. 
The  following  are  fixed  by  the  hardware  and  the  working  gas: 


volume  of  storage  tanks, 


area  of  test  section. 


gas  constant. 


ratio  of  specific  heats, 


l/-  5218.3  ft* 


A - 1 ft3 


ft-lbf 

R - 1715  — ,e>„ 


slug  R 


y - 1.4 


2.  Find: 


The  available  test  time  as  a function  of  test  Mach  number,  initial  storage  pressure 
and  temperature,  and  tunnel  total  pressure  . 

3.  Assumptions: 

A series  of  assumptions  leads  to  an  easy  and  accurate  solution. 

Assume  that  the  air  in  the  tank  expands  lsentroplcally  causing  the  pressure  and  the 
density  to  be  related  by 


4 ’ft)’ 


(5) 


Assume  that  the  heat  sink  cools  from  the  upstream  end,  maintaining  the  discharge  at  the 
initial  storage  temperature,  T «■  T This  assumption  is  well  supported  by  experimental 

O Sj 

data  from  the  tunnel.  Further,  we  will  assume  that  the  stagnation  temperature  is  constant, 
which  really  means  we  sssudm  the  control  valve  is  perfect  and  we  are  neglecting  the 
unsteady  features  associated  with  starting  and  stopping  the  tunnel  flow.  The  control 
valve  on  the  upstream  aide  of  the  stilling  chasfcer  is  choked  during  a run,  and  a run  ends 
when  this  ceases  to  be  true. 
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4.  Solution: 


For  quasi-steady  flow  the  negative  rate  of  change  of  mass  In  the  storage  tanks  must 
equal  the  mass  flow  rate  through  the  test  section  by  conservation  of  mass: 


dp  , 

-v  ~i r~  ■ p VA  - 2-  p0  RT)  m 

o 


y RTq  MA 


*(%,) 


Using  the  perfect  gas  law  and  the  energy  equation, Eqn  (7)  becomes 


■(%,) 


and  Integrating,  yields : 

ps. 

r\  „ 


* t 


or  rearranging, 


0 
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If  we  now  apply  the  lsentroplc  relation  between  pressure  and  density  In  the  storage  tank. 


Eqn  (11)  becomes 


At  the  final  time  (end  of  the  run),  the  flow  between  the  storage  tank  and  the  stilling 
chamber  Is  barely  choked,  which  means  at  t. 


Substituting  Into  Eqn  (12)  yields 


Mv  RT 


I*. L ./Vf UY'\ 

m y1  2 "y  po  I \ps  / V 2 / J* 


By  making  the  appropriate  substitutions  of  H * 2.02,  P - 111  psla,  P - 33  psla,  and 

8j  O 

T « 530°R, the  available  test  time  per  compressor  pumping  was  found  to  be  15.75  seconds. 

sl 

Since  only  force  data  was  required  on  the  cadet  projects,  the  time  was  adequate. 

This  solution  Is  especially  gratifying  because  by  running  the  tunnel  we  were  able  to 
compare  the  actual  run  time  with  that  computed  using  Eqn  (14).  This  comparison  Is  shown 
In  Figure  5.  The  actual  test  time  Is  seen  to  be  within  a fraction  of  a second  of  the 
calculated  time. 

Another  calculation  had  to  be  made  to  determine  the  time  required  to  charge  the 
storage  tanks  to  110  pslg  given  the  substitute  compressor  first-stage  cylinder  displace- 
ment and  the  compressor's  rotational  speed.  Can  your  students  do  thlc  one? 


C.  Stagnation  Pressure  Controller  Design 


The  trlsonlc  tunnel  has  operated  since  195B  with  manual  control  of  the  stagnation 
pressure.  During  the  years,  some  components  of  an  automatic  system  have  been  acquired. 
Completion  of  the  system  has  became  essential  to  reduce  the  total  pressure  fluctuations 
In  the  tunnel  by  an  order  of  magnitude  and  make  the  trlsonlc  tunnel  perform  up  to 
standards. 
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■ 15.3  sec 


NOTE:  START  IDENTIFIED  BY  SUDDEN 

REDUCTION  IN  STATIC  PRESSURE 

STOP  IDENTIFIED  BY  SUDDEN 
INCREASE  IN  STATIC  PRESSURE 


TIME  (sec) 


Figure  5.  Blow-Down  Tunnel  Run  Time 


Several  components  of  the  system  are  on  hand,  and  the  englneera  problem  is  to 
select  the  remaining  components  to  complete  the  system.  Again, the  details  follow.  In 
this  case  It  was  necessary  to  delete  soma  of  the  secondary  calculations  supporting  the 
assumptions  and  present  results  for  only  one  Mach  number  to  fit  the  solution  Into  a 
concise  presentation. 


STATIC  PRESSURE  (PSIG) 
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Figure  6.  Wind  Tunnel  Stagnation  Pressure  Controller 


1.  Given: 


The  available  equipment  Is  shown  by  the  broken  line  In  Figure  6.  Several  pressure 
transducers  are  available.  It  Is  also  known  that  corner  dally  available  servo  amplifiers 
have  gains  from  2.5  to  250  mllll-amps  per  volt.  The  Important  parameters  are: 


gain  of  hydraulic  servo  valve 
piston  area  of  hydraulic  actuator 
area  of  open  control  valve 


stroke  of  valve  and  actuator 

*x* 

gain  of  control  valve  « — 


k 

« r uj 


q 15  Mln-MA 
A - 10.16  loP 

*i*max  ’ 81 
* • 6 In 


volume  of  stilling  chamber 


Ku  - 1.33  In 


V - 27  ft* 
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area  of  throat  at  M ■ 4.5 
nominal  storage  pressure 
gain  of  coeuercial  servo-amp 

2.  Find: 

The  engineering  problem  is  to  determine  the  proper  gain  of  the  pressure  transducer 
to  achieve  the  tightest  possible  control  without  oscillation. 

3.  Solution: 


*e  “ 


16.36 


ft* 


P “ 611  psla 


Ksa  " 2‘5-250  V 


The  first  step  is  to  estimate  the  time  constants  associated  with  the  response  of: 

(1)  the  valve-body/oll -filled  actuator,  (2)  the  servo-valve,  (3)  the  transducer,  (4) 

stilling  chamber  pressure  (to  a change  in  valve  position),  and  (5)  the  change  in  pressure 

in  the  storage  tanks.  Without  proof,  the  conclusion  of  these  estimates  is  that  the  valve 

body,  servo-valve,  and  transducer  are  all  very  fast  compared  to  the  response  of  stilling 

chamber  pressure  to  a change  in  valve  position.  The  response  of  stilling  chamber  pressure 

to  valve  position  is  very  fast  compared  to  the  slow  decay  of  pressure  in  the  tanks. 

Therefore,  it  is  appropriate  to  consider  the  transducer  and  servo-valve  to  be  essentially 

instantaneous  and  to  neglect  mechanical  vibration  of  the  valve  body/ actuator.  It  is 

also  appropriate  to  assume  the  supply  pressure  to  be  constant  for  purposes  of  estimating 

the  remaining  required  gains.  It  is  a fact  that  the  tank  pressure  does  slowly 

decrease  causing  the  actual  stagnation  pressure,  PQ  , to  be  slightly  less  than  the 

CMmanded  pressure,  PQ  . This  slight  difference  between  commanded  and  actual  pressure 
c 

causes  the  servo-amp  to  have  a small  output,  causing  the  servo  valve  to  pass  some  oil; 
this,  then,  causes  the  valve  to  open  slowly  as  the  run  progresses.  After  an  appropriate 
feedback  gain  is  found,  it  is  a straightforward  problem  to  estimate  the  pressure  differ- 
ence required  to  drive  the  valve  at  the  velocity  required  to  maintain  constant  stagnation 
pressura.  This  calculation  uses  the  same  baalc  approach  as  in  the  previous  problem  and 
will  not  be  done  here. 
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This  loog  list  of  results,  then,  brings  us  to  the  problem  an  undergraduate  should  be 
able  to  solve.  Find  the  proper  gain  to  make  the  system  stable  and  non-osclllltory. 

I 

Figure  7 shows  the  Idealized  feedback  loop  to  be  analyzed.  In  that  loop,  note  that  the 
actuator  contributes  to  the  dynamics  of  the  system  because  the  output  of  the  actuator  Is 


proportional  to  the  Input  Integrated  over  time.  Clearly  the  transfer  function 
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If  we  differentiate  the  perfect  gas  law,  P - o RT  - 77  RT  .and 

o (/  o 

substitute,  l.e. 


dPo  dMRTo 
dt  dt  1/  1/ 


★ * 


.y*-i 

£(y-l) 


Rit  the  differential  equation  In  standard  form, 


•px(v^y(^)2<y' 

Take  the  Laplace  transform, 

r *.*(r  «„)*/ 2 \^j 

v»l 

V>  «)*/  2 V(y_1)  * 
i/  U1/  M 


_a±l 

>/  2 \2(V-D 


And  the  sought-after  transfer  function  Is 

p ly  rt  Y*/  VyH (p»i 

».<•>  -4-^ 
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Further,  If  we  let 


-ytl 


1 

pn(.)  A,*  T 


(-*f) 


The  forward  transfer  function  then  becomes 
P , 

v K K — — 1 
sa  q ax  ^ T 


R) 


and  the  feedback  transfer  function  becomes 


H(s)  - Kfp  . 


H(s)  G(s) 


K K r*  A K, 

sa  q ax  Af  r fP 


R) 


R] 


The  closed-loop  transfer  function  becomes 


V*>  _ at 


Ej(s)  1+  G(s)  H(s) 


Therefore,  we  want  to  find  the  valve  of  the  loop  gain  K such  that 
the  roots  of  the  equation. 


1 ♦ G(s)  H(s)  - 0 , 
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Just  11«  on  th«  negative  reel  axis  so  that  any  Increase  In  K will 
imaginary  parts.  Combining  Eqn  (23)  and  (24)  Into  Eqn  (26)  gives 

s»  + - + K - 0 
T 


cause  them  to  have 


(27) 


For  K “ this  becomes  a perfect  square  leading  to  equal,  real  negative  roots . 

Increasing  the  gain  above  K “ leads  to  a pair  of  complex  conjugate  roots  and  an 
oscillating  solution.  Reducing  the  gain  below  K - gives  less  than  achievable 
performance. 


Now  Equation  (24)  is  used  to  find  the  minimum  required  gain  for  the  pressure 
transducer.  If  a transducer  with  greater  than  the  minimum  gain  Is  available,  then  no 
other  amplifiers  other  than  the  servo-ampllfler  are  required. 


K.. K.  K«  i % 

*8 


- K 

*7® 


Pressure  transducers  with 


(29) 


Kfp  ‘ 0.000015  J2i£  . o.oi5 

or  greater  are  readily  available  In  the  laboratory  and,  tharefore,  the  only  new 
equipment  required  la  the  servo-ampllfler. 

III.  What  Experiences  Build  Confidence? 

How  can  we  build  technical  confidence  In  our  students?  To  answer  this  question  I 
turned  a bit  Introspective  and  asked  how  t knew  each  of  the  above  solutions  would  work. 
The  answer  is  experience  - not  years  of  It  - but  specific  experiences  that  built 
confidence  In  each  of  the  analytical  techniques  used. 
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The  turning  vane  problem  reminded  me  o£  a 1959  mechanical  engineering  laboratory 
experiment  that  was  part  of  a required  course  which  I detested.  After  all,  I aspired 
to  be  a theoretician  and  a researcher,  and  knowledge  of  the  torque  of  a Pelton  Wheel 
(an  Impulse  water  turbine)  could  hardly  contribute  to  this  lofty  goal.  I did  remember, 
however,  that  the  calculated  torque  and  the  measured  torque  were  within  a few  percent  of 
each  other.  The  calculation  was  based  on  the  momentum  equation  — the  same  principle 
applied  to  the  turning  vane  problem. 

The  wind  tunnel  blow-down  time  problem  exercised  considerable  knowledge  of  compres- 
sible flow  and  basic  thermodynamics.  Several  physical  observations  supported  confidence 
In  this  solution.  The  detested  mechanical  engineering  laboratory  had  examples  of  air- 
flow measurement  with  choked  and  un-choked  nozzles,  and  I recalled  that  the  correction 
to  the  theoretical  calculations  were  on  the  order  of  2-5%.  Another  feature  of  the  blow- 
down time  problem  Is  that  it  Is  fundamentally  unsteady,  and  an  assumption  of  quasi-steady 
flow  had  to  be  made  to  get  the  solution.  My  confidence  that  this  assixaptlon  was  adequate 
was  based  on  a single  homework  problem  in  Shapiro  which  showed,  when  properly  worked, 
that  the  unsteady  terms  are  only  Important  when  a significant  fraction  of  the  flow  has 
kinetic  energy. 


In  the  stagnation  pressure  control  problem,  confidence  In  the  linearized  analysis  and 
the  use  of  Laplace  Transforms  dates  from  a 1963  course  In  analog  computing  at  the 
University  of  Arkansas  Graduate  Institute  of  Technology  (a  city  extension  campus  located 
In  Little  Rock).  In  that  course  we  modeled  many  different  systems,  both  linear  and  non- 
linear systems.  The  observation  that  the  behavior  of  a non-llneer  system  around  the 
operating  point  Is  very  Closely  predicted  by  the  linear  system  gives  me  absolute  confi- 
dence that  selection  of  equipment  based  on  this  analysis  of  the  trlsonlc  tunnel  will 
yield  an  adequate  controller. 

Since  the  earliest,  confidence-building  experlencea,  many  engineers  have  had  rein- 
forcing experlencea  which  maintained  confidence  in  fundamentals.  However,  It  wee  these 
first  verlflcstions  of  the  fundementals  that  oterted  the  confident  assault  on  problems 
of  lncressing  difficulty. 
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The  critical  first  Impression  of  new  technical  material  Is  always  based  on  at  least 

Ti 

one  comparison.  Possible  candidates  for  comparisons  are:  (1)  analytical  solutions, 

(2)  laboratory  experiments,  (3)  simulations,  and  (4)  numerical  solutions.  Each  of  these 
can  be  exact  or  approximate.  It  seems  clear  that  a large  number  of  successful  comparisons 
would  tend  to  bolster  confidence  In  the  material. 

A personal  note  of  caution  on  making  these  comparisons'  any  comparison  must  be  done 
carefully  and  the  reasons  for  differences  correctly  Identified.  Vague  statements  of 
"experimental  error"  or  "approximations  In  the  theory"  are  not  reasons!  The  exact 
assumptions  or  experimental  measurement  Inadequacies  contributing  to  the  difference  must 
be  Identified.  Our  students  must  constantly  place  relative  values  on  the  technical 
material  we  give  them.  If  we  don't  convincingly  Identify  the  correct  source  of  differ- 
ences In  these  comparisons,  we  fall  to  give  them  the  Information  necessary  to  make  these 
value  judgments;  In  this  way,  we  actually  prevent  them  from  developing  conf ldence  In 
fundamentals. 


Svbol« 


A - area 

G - forward  transfer  function 
H - feedback  transfer  function 
K - sensitivity  factor 


i 


i 


F - force 

A - mass  flow  rate 
M - Mach  number 
P - pressure 
R - gas  constant 

s - Laplace  variable  (frequency  domain) 
t - time 

V - volume 

V - velocity 

x - Cartesian  coordinate 
y - Cartesian  coordinate 
3 - Cartesian  coordinate 


p - density 

y - ratio  of  specific  heats 
x - stroke  of  actuator 


Subscripts  and  superscripts 

f - final 
s - storage  tank 

0 • stagnation  conditions 

* - throat  or  sonic  conditions 

1 - Initial  condition 


142 


£ V 


/ V 


> 

USAFA-TR-79-1 


AERODYNAMIC  CHARACTERISTICS  OF  AN 

★ 

UNSTEADY  SEPARATED  FLOW 

M.  S.  Francis**,  J.  E.  Keesee***,  J.  D.  Lang+,  G.  W.  Sparks"’’!  G.  E.  Sissont"f" 

Extended  Abstract 

During  flight  vehicle  maneuvers  at  high  angles  of  attack,  dynamic  stall  or  other 
complex  flow  phenomena  Involving  time-varying  separation  may  occur  resulting  in  aero- 
dynamic loads  which  differ  significantly  from  those  predicted  by  steady  flow  consider- 
ations. Future  designs  of  advanced  maneuverable  flight  vehicles  may  well  seek  to  in- 
corporate the  advantageous  effects  of  these  highly  energetic  flows  in  a positive  way 
to  enlarge  the  effective  flight  envelope.  The  capability  to  stabilize  and  control  flow 
behavior  necessary  for  productive  improvements  in  performance  must  be  preceded  by  a more 
complete  knowledge  of  the  mechanics  of  unsteady  separation.  The  investigation  discussed 
in  this  paper  is  directed  at  a realistic  modeling  of  these  types  of  flows  through  an 
improved  understanding  of  the  role  of  the  dominant  physical  mechanisms. 


An  experiment  which  provided  a simplified  and  controlled  unsteady  separation  having 
a two-dimensional  flow  geometry  was  devised  using  a NACA  0012  airfoil  with  an  oscill- 
ating fence-type  spoiler  located  at  mid-chord  on  one  surface  lmsersed  in  a subsonic  wind 
tunnel  free  stream.  With  the  capability  of  adjustable  spoiler  mean  height,  oscillation 
amplitude  and  frequency  characteristics,  this  configuration  provided  for  control  of  the 
size  of  the  separated  region  while  fixing  the  separation  point  with  respect  to  the  free- 
3 t re am  coordinate  direction. 

A characteristic  feature  of  the  resultant  flow  field  observable  for  a wide  range  of 
flow  conditions  was  the  generation  and  subsequent  roll  up  of  a large  amount  of  vorticity, 
While  this  variable  represents  a physically  significant  descriptor  of  regions  of  vis- 
cous interaction  in  real  fluids,  it  is  a parameter  whose  experimental  determination  is 
often  difficult.  In  the  experiments  described  below,  the  time-varying  spatial  vorticity 
distribution  associated  with  the  separated  region  was  Inferred  from  both  hot-wire  ane- 
mometer and  two-color  laser  doppler  velocimeter  measurements  using  a spatial  contour 
Integration  technique  (Ref  1).  In  this  procedure,  measurements  of  select  Instantaneous 
velocity  component  values  obtained  over  a mesh  of  spatial  sampling  locations  were  first 

*A1AA  17th  Aerospace  Sciences  Meeting,  AIAA  Paper  No.  79-0283 

**Captaln,  USAF,  Chief,  Mechanics  Division,  FJSRL 

Captain,  USAF,  Research  Associate,  FJSRL 

+ Major,  USAF,  4950/FFMS , W-PAFB 
Captain,  USAF,  Instructor  of  Aeronautics,  DFAN 
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ensemble  averaged  (a  phase-locked  average  correlated  to  spller  position)  using  a digital 

computer.  The  results  were  employed  to  calculate  the  circulation  about  elements  of  a 
planar  array  of  spatial,  rectangular  contours  like  that  shown  in  Figure  1.  The  resultant 
value  of  vorticity  for  any  given  phase  point,  spatially  averaged  over  the  "i"-th  element 
of  the  mesh  was  then  determined  by  the  relationship 


<4  <t)>  - i f U(x,t)  • di(x) 

1 Ct 

p 

where  £ is  the  vorticity  component  perpendicular  to  the  plane  of  the  integration  array  and 
At  is  the  area  of  the  mesh  element.  The  conversion  of  data  time  histories  to  instan- 
taneous spatial  distributions  allowed  for  a more  useful  description  of  results  through 
the  generation  of  iso-vorticity  contour  maps. 

Resultant  vorticity  distributions  were  then  compared  with  airfoil  surface  pressure 
measurements  obtained  from  sixteen  miniature  transducers  located  in  the  wing's  interior 
in  an  effort  to  generate  a correlation  for  modeling  purposes.  A schematic  of  the 
complete  experimental  apparatus  is  provided  in  Figure  2. 


Investigation  of  the  global  ensemble-averaged  vorticity  field  generated  by  simple 
harmonic  spoiler  oscillations  has  revealed  that  the  separation  zone  is  not  merely  a 
growing  and  shrinking  "bubble"  as  previously  thought  (Ref  2),  but  is  a dynamically 
evolving  vortex-like  structure  bearing  a strong  resemblance  to  the  one  encountered  in 
the  dynamic  stall  flow  field  (Ref  3)  (Figure  3).  Iso-vorticity  contour  maps  of  the 
resultant  field  in  the  outer  shear  layer  region  show  how  the  concentrated  vortex  which 
forms  on  the  spoiler  upstroke  begins  to  convect  and  diffuse  as  the  spoiler  retracts. 

An  example  of  this  behavior  is  depicted  in  Figure  4 for  several  phase  points  under 
identical  flow  conditions.  While  the  peak  vorticity  levels  are  observed  to  be  primarily 
a function  of  the  freestream  Reynolds  number,  the  development  of  the  flow  geometry, 
including  a characteristic  lag  in  the  separated  region's  growth  is  primarily  dependent 
on  the  dimensionless  frequency,  k.  The  integral  effect  of  the  variation  of  circulation 
with  time  (phase  angle)  over  the  entire  region  is  useful  in  discussing  temporal  lag 
effects. 

The  behavior  of  the  corresponding  surface  pressure  coefficient  values  (Figure  5) 
indicates  that  the  growth  characteristics  of  a suction  peak  (which  appears  in  the  region 
downstream  of  the  spoiler  on  its  upstroke)  correlate  closely  with  those  of  the  outer 
shear  layer  at  least  for  low  values  of  frequency.  The  earlier  development  of  the  neg- 
ative pressure  peak  at  dimensionless  frequencies  near  and  above  unity  is  attributed  to 
an  energetic  flow  reversal  near  the  surface  which  precedes  the  development  of  the  outer 
vortex  in  these  cases. 
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A a lmpl If lection  made  by  correlating  the  'center  of  vortlclty'  of  the  distribution 
with  the  surface  loading  Is  also  discussed.  A simple  discrete  vortex  model  Is  pre- 
sented to  describe  the  effects  of  vortex  entrapment  and  subsequent  convection  on  the 
pressure  distribution. 
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Figure  1.  Contour  Array  Geometry  for  the  Outer  Shear  Layer 
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Figure  4.  Iso-Vortlclty  Contours  In  the  Separated  Region, 
frequency  * 19.85  Hz,  * 18.4  m/sec  (distance 
between  lines  - 400  per  sec,  unless  noted) 
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Figure  4 (cent).  Iso-Vortlclty  Contours  In  the  Separated 
Region,  frequency  ■ 19.65  Hz,  U ■ 18.4 
m/sec  (distance  between  lines  - 400  per 
sec,  unless  noted) 


